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This  study  was  designed  to  test  the  hypothesis  that  thermal  tempering  by 
compressed  air  is  the  most  effective  strengthening  method  for  layered  feldspathic 
porcelain  structures  and  to  determine  the  influence  of  cooling  method,  initial 
tempering  temperature,  thermal  contraction  mismatch,  and  specimen  thickness  on 
transient  and  residual  stresses  produced  within  dental  ceramic  structures. 

Bilayered  porcelain  disks  with  a  thermal  contraction  mismatch  between  -1.5 
and  +3.2  ppm/°C  and  a  thickness  of  2  or  6  mm  were  prepared  and  thermally 
tempered  in  air,  helium,  and  silicone  oil.  The  transient  temperature  distributions 
were  calculated  by  analytical  and  finite  element  methods  and  were  incorporated  in 
finite  element  and  viscoelastic  thermal  stress  models  to  determine  the  transient  and 
residual  stresses  produced  in  each  disk.  The  finite  element  method  was  also 
employed  to  determine  the  influence  of  ceramic  thickness,  loading  orientation, 
thermal  contraction  mismatch,  initial  tempering  temperature  and  cooling  rate  on  the 
stress  distribution  within  a  dental  ceramic  molar  crown.  Statistical  analysis  revealed 
that  forced  convective  cooling  of  disk  specimens  in  air  resulted  in  a  significantly 
smaller  induced  crack  size  and  increased  the  mean  flexural  strength  on  specimens  by 
a  factor  of  two  compared  with  the  specimens  which  were  conventionally-cooled  by 
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free  convection.  However,  the  thermal  stress  models  predicted  transient  tensile 
stresses  far  above  the  failure  level  of  porcelain  for  the  6-mm  thick  bilayered  disks 
subjected  to  free  convective  cooling  and  for  the  2-mm  thick  bilayered  disks  subjected 
to  quenching  at  high  initial  temperatures  in  silicone  oil.  These  predictions  were 
supported  by  experimental  findings.  The  results  of  this  study  support  the 
hypothesis  that  thermal  tempering  by  compressed  air  is  the  most  effective 
strengthening  method  for  layered  feldspathic  porcelain  structures.  However, 
variables  such  as  specimen  dimensions,  initial  tempering  temperature,  heat  transfer 
coefficient,  and  thermal  contraction  mismatch  have  a  significant  influence  on  the 
effectiveness  of  tempering.  The  most  effective  tempering  treatment,  which  results  in 
maximum  flexure  strength  and  minimum  structural  distortion  by  viscous  flow  and 
air  blasting,  was  achieved  by  forced  convective  cooling  in  air  using  initial  tempering 
temperatures  between  750°C  and  850°C. 
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CHAPTER  I 
INTRODUCTION 

Strengthening  of  Glass 

Influence  of  Flaws  on  Strength 

Microcracks  and  critical  flaws  which  act  as  stress  concentrators,  reduce  the 
strength  of  brittle  materials.  The  interaction  of  these  flaws  with  applied  forces  on 
structures  can  cause  either  premature  fracture  or  crack  growth  which  endangers 
the  integrity  of  these  structures.  The  theoretical  tensile  strength  of  a  flawless 
glass  can  be  as  high  as  1400  MPa,  but  because  of  the  surface  stress  raisers  which 
usually  are  present,  the  strength  may  be  reduced  to  70  MPa  [1].  Defects  and 
cracks  can  be  produced  from  grinding,  machining,  handling,  impact  forces,  and 
thermal  stresses.  Studies  of  the  failure  mechanisms  in  glass  have  been  conducted 
for  several  hundred  years.  In  1678,  Hooke  included  glass  in  his  experiments 
which  led  to  the  famous  stress-strain  law  for  solid  bodies.  A  historical  review  of 
glass  strength  factors  and  failure  mechanisms  is  given  by  Holloway  [2]. 

Inglis  in  1913  showed  that  for  a  flaw  represented  by  a  long  and  narrow 
elliptical  hole  in  a  glass  plate,  the  stress  at  the  flaw  tip  increases  by  a  factor  of 
2  (C/pc)1/2,    where  C  and  pc  are  the  half  crack  length  (semimajor  axis  of  an 

ellipse)  and  radius  of  curvature  at  the  crack  tip,  respectively  [3].    However, 
because  the  radius  (pc)  of  a  crack  tip  is  generally  of  atomic  dimensions  and  is 

difficult  to  measure,  many  workers  have  assumed  that  it  is  constant  for  cracks  of 
different  length  [4-5].  Griffith  [6]  made  a  very  important  theoretical  contribution 


to  the  understanding  of  brittle  fracture  by  demonstrating  that,  for  a  given  value  of 
applied  force,  a  crack  will  propagate  and  cause  fracture  only  if  it  has  a  length 
greater  than  a  certain  value  called  the  critical  crack  length  and  this  is  determined 
by  energy  considerations.  The  strength  of  brittle  materials  is  significantly 
influenced  by  the  number  and  severity  of  the  flaws  and  cracks,  specifically  those 
that  are  located  at  the  surface  regions  [5-12].  Sapphire  crystals  after  normal 
handling  have  a  flexural  strength  in  the  range  of  350-700  MPa,  but  when  the 
critical  surface  flaws  are  reduced  by  flame  polishing,  the  strength  increases  to 
10,000  MPa  [13].  However,  this  technique  is  impractical  in  many  cases  because 
the  material  will  be  subjected  to  surface  damage  in  use. 

Although  it  is  difficult  to  maintain  the  surface  of  brittle  materials  completely 
free  of  flaws  and  microcracks,  it  is  possible  to  significantly  reduce  the  effect  of 
these  surface  stress  raisers  by  producing  compressive  stresses  in  a  thin  layer  near 
and  parallel  to  the  surface.  These  stresses  which  are  present  in  the  unloaded  state 
should  act  over  a  distance  from  the  surface  which  is  greater  than  the  depth  of  the 
deepest  surface  flaw.  If  a  compressive  stress  of  70  MPa  is  produced  at  the 
surface  of  a  glass  plate  which  has  a  tensile  strength  of  70  MPa,  its  strength  can  be 
increased  by  up  to  a  factor  of  two,  because  any  tensile  forces  applied  at  the 
surface  of  this  plate  must  first  overcome  the  existing  compressive  stresses  before 
causing  any  tension  at  the  surface.  Compressive  stresses  can  be  produced  at  the 
surface  by  the  use  of  a  thermal  tempering  treatment,  which  is  a  well-established 
technique  in  the  glass  industry. 

In  general,  strengthening  by  surface  compression  (tempering)  can  be 
classified  into  two  main  categories:  1)  chemical  tempering,  and  2)  physical 
tempering  [1]. 


Chemical  Tempering 

In  chemical  tempering  techniques  such  as  the  Schott  overlay  process, 
surface  crystallization  and  ion-exchange,  surface  compression  is  produced  by 
modifying  the  chemical  composition  of  the  surface  [14,15].  In  the  Schott 
process,  glass  with  a  relatively  larger  thermal  expansion  coefficient  is  sealed 
within  a  smaller  thermal  expansion  glass.  When  the  composite  glass  structure 
cools  from  an  elevated  temperature,  the  interior,  which  contracts  more  than  the 
exterior,  causes  the  development  of  compressive  stresses  at  the  surface.  This 
process  is  currently  being  used  to  produce  glass-ceramic  dinnerware.  Krohn  and 
Copper  [16]  used  this  technique  to  produce  cladded  fibers  by  placing  a  solid  rod 
of  high  thermal  expansion  glass  inside  a  tube  of  low  expansion  glass.  The  mean 
strength  of  fibers  drawn  from  this  composite  increased  to  290  MPa  compared 
with  138  MPa  for  the  untreated  fibers. 

Surface  crystallization  forms  a  layer  with  a  lower  thermal  expansion 
coefficient  compared  with  the  interior  and  develops  compressive  stresses  in  the 
surface  as  glass  is  cooled.  In  most  glasses,  this  technique  may  weaken  the 
sample  because  it  introduces  tensile  stresses  or  cracks.  Also  it  may  result  in 
surface  opacity  and  poor  optical  properties.  However,  certain  lithium- 
aluminosilicate  systems  can  be  treated  in  this  manner  to  achieve  a  strength  of 
about  690  MPa  [15].  For  strengthening  of  a  0.5-cm  diameter  lithium- 
aluminosilicate  rod  the  optimum  thickness  of  the  compressive  surface  layer  was 
found  to  be  0.07  mm.  When  the  crystallization  process  was  continued  too  long, 
it  produced  large  tensile  stresses  which  resulted  in  failure  of  the  specimen  [5]. 

Toughening  of  glass  by  ion-exchange  is  accomplished  when  a  large  ion  is 
exchanged  for  a  smaller  one  at  temperatures  below  the  glass  transition 


temperature.  As  a  result  of  this  ion  transport,  compressive  stresses  build  up  in 
the  exchanged  region.  In  this  technique,  glass  is  immersed  in  a  molten  salt  bath 
containing  potassium  ions.  Sodium  ions  diffusing  from  surface  into  the  fused 
salt  are  replaced  by  an  equal  number  of  larger  potassium  ions.  The  surface 
compressive  stresses  develop  because  of  the  difference  in  size  between  the  alkali 
ions  which  increase  the  local  density  of  surface-near  regions.  The  larger 
potassium  ions  try  to  expand  the  surface  but  this  is  resisted  by  the  chemically 
unchanged  interior  region  which  puts  the  surface  in  compression  [17].  It  is 
possible  to  combine  ion-exchange  and  surface  crystallization  for  some  glasses  to 
enhance  the  effectiveness  of  chemical  tempering.  Above  the  glass  transition 
temperature,  surface  crystallization  which  occurs  as  a  results  of  ion-exchange  can 
reduce  the  thermal  expansion  coefficient  of  the  surface  layer.  This  develops 
larger  compressive  stresses  at  the  surface  [5]. 

Physical  Tempering 

In  physical  tempering  (thermal  tempering),  rapid  cooling  of  glass  at  an 
initial  temperature  above  the  glass  transition  causes  the  glass  surface  to  cool 
before  the  interior  region.  During  the  initial  stage  of  the  tempering  process,  the 
surface  becomes  rigid  and  contracts,  while  the  interior  is  still  soft  and  can  adjust 
its  dimensions  to  the  surface  contraction.  Cooling  is  most  often  accomplished  by 
jets  of  air  or  sometimes  by  an  oil  bath  [18].  At  the  later  stages  of  tempering,  the 
interior  temperature  drops  below  the  glass  transition  temperature  and  it  will 
contract.  However,  the  surface  is  rigid  and  cannot  adjust  itself  to  the  interior 
contraction.  Since  the  volume  elements  at  the  surface  are  constrained  from  any 
expansion  or  contraction,  compressive  stresses  are  produced  at  the  surface, 
simultaneously  balancing  tensile  stresses  of  smaller  magnitude  which  are 


developed  in  the  interior  regions.  This  effect  can  be  enhanced  by  coating  a  glass 
with  another  one  of  lower  thermal  expansion  coefficient  [5].  Precompression  of 
the  glass  surface  is  based  on  the  unique  behavior  of  glass  in  the  transformation 
temperature  range  during  cooling  where  its  properties  change  from  that  of  a 
viscous  liquid  to  that  of  an  elastic  body  [19-34]. 

The  initial  idea  of  tempering  glass  goes  back  to  the  17th  century  when 
experiments  were  conducted  to  produce  the  so-called  "Prince  Rupert's  Drops" 
[20].  To  produce  these  drops,  molten  glass  is  subjected  to  a  water  quench  which 
produces  a  tear-shaped  glass  drop  with  a  long  tail.  Since  the  exterior  surface  of 
this  molten  drop  cools  before  the  interior,  large  compressive  stresses  are  induced 
on  the  surface,  while  the  interior  is  under  tension.  As  a  result,  the  body  of 
"Prince  Rupert's  Drop"  becomes  so  strong  that  it  cannot  be  broken  by  relatively 
large  compressive  forces  applied  by  a  hammer.  However,  the  drop  will  burst  in 
a  quite  rapid  and  violent  manner  when  tensile  stresses  are  applied  by  twisting  its 
tail.  This  is  similar  to  the  less  explosive  fracture  which  occurs  when  a  tempered 
car  windshield  breaks  into  a  large  number  of  fragments  a  few  millimeters  in  size 
[17]. 

One  of  the  earliest  tempering  experiments  for  glass  plates  was  performed  by 
Siemens  [35],  who  introduced  the  method  of  quenching  by  contact  with  solid 
surfaces.  However,  most  of  the  studies  on  controlling  thermal  stresses  in  glass 
prior  to  1 920  were  related  to  annealing  in  which  the  residual  stresses  in  glass 
products  were  removed  by  slow  cooling  in  an  annealing  oven.  In  1920  the 
method  of  quenching  flat  glass  with  impinging  air  jets  was  invented  in  France, 
and  in  early  1930  thermally  tempered  safety  glass  began  to  be  widely  used  in 
European  automobiles.  When  a  toughened  windshield  glass  breaks  in  the 
presence  of  tensile  stresses  at  the  interior,  it  breaks  into  small  pieces  which  are 


relatively  harmless.  In  tempered  automotive  glass,  one  usually  seeks  high 
strength  and,  upon  fracture,  small  fragments  to  minimize  physical  injury  from 
glass  shards.  However,  in  the  tempered  glass  used  in  some  architectural 
applications  and  in  ovenware,  increased  strength  and,  upon  fracture,  large 
fragments  are  desirable.  The  enhancement  in  the  strength  of  tempered  glass  is 
determined  by  the  magnitude  of  compressive  stresses  in  its  surface  while  the 
propensity  to  break  into  many  small  fragments  is  believed  to  be  governed  by  the 
magnitude  of  the  interior  tension.  Perhaps  the  most  extensive  analytical  and 
experimental  studies  on  tempering  of  glass  in  the  United  States  have  been 
conducted  by  Robert  Gardon.  In  a  recent  study  [24],  Gardon  has  shown  that  the 
strength  and  fracture  characteristics  of  tempered  glass  can  be  controlled 
independently  by  programming  the  heat  transfer  coefficient  of  air  flow  during  the 
tempering  process.  In  this  technique,  which  is  called  "modulated  quenching," 
the  ratio  of  surface  compression  to  interior  tension  is  controlled. 

Comparison  Between  Chemical  and  Physical  Tempering 

Doremus  [5]  has  reported  that  glass  can  be  strengthened  by  a  factor  of  17 
for  surface  crystallization,  10  for  ion-exchange,  22  for  combined  ion-exchange 
and  surface  crystallization,  and  6  for  thermal  tempering.  However,  the 
effectiveness  of  a  tempering  technique  depends  on  factors  such  as  its  ability  to 
maintain  the  strength  over  a  long  period  of  time  under  a  variety  of  loading 
conditions,  the  time  and  the  cost  required  for  strengthening,  the  convenience  of 
the  process,  and  the  simplicity  of  the  method. 

Compared  with  thermal  tempering,  chemical  tempering  provides  a  greater 
degree  of  flexibility  in  overcoming  physical  design  problems  such  as  size  and 


shape  restrictions.  Since  the  chemical  exchange  process  is  carried  out  at  lower 
temperatures,  there  is  no  danger  of  distortion  resulting  from  viscous  flow  of 
glass.  However,  compared  with  the  thermal  tempering,  chemically  strengthened 
glass  is  more  susceptible  to  abrasion,  is  less  flexible  in  terms  of  chemical 
composition,  is  more  expensive  and  the  processing  is  longer. 

Shown  in  Fig  1.1  is  the  difference  between  the  residual  stress  profiles 
produced  by  ion-exchange  and  thermal  tempering  [1].  Although  the  chemical 
tempering  produces  greater  compressive  stresses  at  the  surface  and  smaller  tensile 
stresses  at  the  interior,  the  compressive  stresses  are  developed  only  in  a  thin  layer 
near  the  surface  and  the  tensile  stresses  at  the  interior  regions  are  extended 
dangerously  close  to  the  surface  layer  [1]. 

Major  Findings  from  Thermal  Tempering  of  Glass 

Gardon  has  employed  cinematography  to  record  the  transient  photoelastic 
fringe  patterns  corresponding  to  the  transient  stress  distributions  in  a  6.1  mm  thick 
glass  plate  during  the  tempering  process.  The  stress  profiles  for  the  first  20 
seconds  of  quenching  are  shown  schematically  in  Fig.  1.2.  For  a  constant  heat 
transfer  coefficient  of  222  W/m2-K,  stresses  are  influenced  by  the  initial  tempering 
temperature.  At  the  final  stage  of  tempering,  the  stress  distribution  is  parabolic 
with  the  surface  in  compression  and  the  interior  in  tension.  However,  it  is 
important  to  notice  that  the  surface  stresses  at  the  initial  stage  of  tempering  may  be 
tensile  or  compressive  depending  on  the  initial  tempering  temperature  (Tj).  At  low 
Ti  values  of  616<>C  and  648°C  (cases  A  and  B),  the  surface  initially  sustains  tensile 
stress  which  subsequently  changes  to  a  state  of  compression,  while  at  a  higher 
initial  temperature  of  738°C  (case  C),  it  passes  directly  from  an  initially  stress-free 
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Fig.  1.1    Residual  stress  distribution  in  glass  resulting  from  physical  (thermal) 
and  chemical  tempering  [1]. 
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Fig.  1.2   Transient  stress  distribution  in  glass  plate  during  tempering  from 
initial  temperatures  of  616'C  (A),  648*C  (B)  and  13  ft  (C)  [201. 
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condition  into  compression.  The  effect  of  heat  transfer  coefficient  and  specimen 
thickness  on  midplane  tension  in  a  flat  glass  has  also  been  reported  by  Gardon  [20] 
using  the  birefringence  technique.  At  a  constant  initial  temperature  of  650°C,  when 
the  heat  transfer  coefficient  in  air  was  increased  from  12.5  W/m2-K  (free 
convection)  to  552  W/m2-K  (forced  convection),  the  midplane  tension  was 
increased  by  a  factor  of  6.9.  For  the  same  initial  temperature,  an  increase  in 
specimen  thickness  from  2.1  mm  to  26.5  mm  resulted  in  an  increase  in  the 
midplane  tension  by  a  factor  of  8.5.  Although  the  surface  compression  which 
enhances  the  strength  of  tempered  specimens  had  been  taken  previously  as  2.0  to 
2.2  times  greater  than  the  midplane  tension,  the  actual  ratio  of  surface  compression 
to  midplane  tension  may  be  considerably  different  for  different  materials  and 
cooling  conditions. 

Gardon  [20]  has  also  shown  that  for  a  constant  glass  thickness  and  heat  transfer 
coefficient,  a  limiting  value  for  the  initial  cooling  temperature  exists  such  that  an 
increase  in  that  temperature  before  quenching  will  not  result  in  any  further  increase 
in  the  magnitude  of  induced  surface  compression  (Fig.  1.3).  However,  as 
previously  stated,  the  development  of  transient  tensile  stresses  at  the  initial  stage  of 
tempering  depends  on  the  initial  temperature.  For  moderate  quenching  conditions 
with  a  heat  transfer  coefficient  of  222  W/m2-K,  which  is  representative  of 
industrial  practice  for  tempering  of  6-mm-thick  glass,  an  initial  temperature  of 
648°C  produced  tensile  stresses  at  the  surface  during  the  initial  stage  of  tempering 
(Fig.  1.2),  while  tempering  from  an  initial  temperature  of  738°C  did  not  produce 
any  surface  tensile  stresses.  Fig.  1.3  also  indicates  that  for  a  fixed  value  of  initial 
temperature  and  specimen  thickness,  an  increase  in  the  heat  transfer  coefficient  will 
result  in  a  greater  degree  of  tempering.  These  results  show  that  the  beneficial  effect 
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Fig.   1 .3    Degree  of  temper  (nm/cm)  in  a  glass  plate  as  a  function  of  initial 
temperature  and  heat  transfer  coefficient  [20]. 
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of  quenching  is  strongly  related  to   the  heat  transfer  coefficient  and  the  initial 
temperature  employed  in  the  tempering  treatment. 

Strengthening  of  Ceramics 

By  comparison  of  the  color  and  texture  of  many  of  the  prehistoric  tools  and 
weapons  made  from  Ocala  chert,  anthropologists  have  empirically  determined 
that  prehistoric  Indians  performed  heat  treatment  on  this  ceramic  to  obtain 
smoother  fracture  surfaces  and  sharper  edges  [36,37].  Ceramic  materials  in 
general  have  desirable  properties  for  structural  design  such  as  thermodynamic 
stability,  resistance  to  corrosion,  low  density,  high  melting  point  and  reliable 
performance  at  high  temperatures.  However,  they  exhibit  a  high  degree  of 
brittleness,  low  impact  resistance,  and  high  sensitivity  to  surface  damage  which 
may  result  in  their  failure  without  prior  warning. 

Based  on  the  success  with  the  thermal  tempering  of  glass  which  produces 
beneficial  compressive  stresses  at  the  surface,  substantial  progress  has  also  been 
made  in  strengthening  ceramics.  Mechanical  strength,  resistance  to  thermal 
shock  failure  and  impact  of  some  polycrystalline  ceramics  such  as  alumina  have 
been  improved  by  thermal  tempering  [38].  Thermal  tempering  can  assist  in  the 
more  widespread  use  of  ceramics  in  engineering  applications  such  as  in  energy 
conversion  and  within  the  aerospace  industry.  High  temperature  strength  and  the 
smaller  density  of  ceramics  compared  with  metal  can  result  in  a  significant 
improvement  in  performance,  thermodynamics  and  overall  efficiencies  in  heat 

engines. 

Several  authors  have  discussed  factors  which  influence  the  design 
procedures  required  for  the  successful  applications  of  ceramics  in  thermal  and 
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structural  equipment  such  as  gas  turbines,  heat  exchangers,  automobile  engines 
and  refractory  metals  [39-43].  Although  a  wide  range  of  strengthening 
mechanisms  or  processes  such  as  annealing  and  minimizing  grain  size  have  been 
identified  for  ceramics,  in  most  cases  the  improvements  have  been  unimpressive 
except  for  the  results  obtained  by  using  compressive  surface  stresses  such  as 
glazing  and  thermal  tempering.  Most  pottery  and  chinaware  are  glazed  with  a 
layer  that  has  a  lower  thermal  expansion  coefficient.  During  cooling,  the  interior 
tends  to  contract  more  than  the  glaze,  which  put  the  glaze  layer  in  compression 
and  the  interior  into  tension. 

Major  Finding  from  Thermal  Tempering  of  Ceramics 

Studies  have  shown  that  tempering  improves  the  strength  of  various  ceramic 
materials  [44-46].  Kirchner  et  al.  [47],  conducted  experiments  by  quenching 
both  glazed  and  unglazed  alumina  rods  (96%  alumina),  using  gaseous  cooling 
media  such  as  air,  helium,  and  CO2  (forced  convective  cooling),  and  liquid 
cooling  media  such  as  silicone  oil  with  various  kinematic  viscosity,  kerosene  and 
water  (free  convective  cooling).  The  mean  value  of  the  flexural  strength  using  a 
four-point  bending  test  for  the  untreated  (control)  specimens  was  322.7  MPa. 
The  corresponding  values  for  specimens  quenched  in  low  viscosity  silicone  oil 
from  an  initial  temperature  of  1600°C  and  by  forced  helium  from  an  initial 
temperature  of  1500°C  were  728.3  MPa  and  522.4  MPa,  respectively. 
However,  quenching  the  specimens  in  water  caused  thermal  shock  failure. 
Shown  in  Table  1.1  is  a  summary  of  the  results  obtained  when  quenching 
alumina  in  various  media. 
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Table  1 . 1     Flexural  strength  of  96%  alumina  quenched  in  various  media  [471. 


QUENCHING  MEDIUM 


QUENCHING  NO.  OF  STRENGTH 

TEMPERATURE  (°C)       SPECIMENS  (MPa) 


As-received  (control) 

19 

322.7 

Silicone  oil  (y=  5.0  cs) 

1600 

5 

728.3 

Silicone  oil  (y=  12,500  cs) 

1550 

5 

593.8 

Water  (Thermal  Shock) 

1600 

3 

35.9 

10%  oil  +  90%  H2O 

1550 

3 

640.7 

Forced  ail- 

1500 

5 

533.8 

Forced  helium 

1500 

5 

522.4 

Significant  improvements  were  also  obtained  in  tensile  strength,  thermal 
shock  resistance  and  impact  resistance  of  tempered  compared  with  slowly  cooled 
specimens.  Thermal  shock  tests  were  performed  by  quenching  the  specimens  in 
water  at  room  temperature  and  measuring  the  remaining  flexural  strength  after 
one  thermal  shock  cycle.  Impact  resistance  was  measured  by  dropping  a  weight 
from  progressively  greater  heights  onto  a  specimen  until  a  height  was  reached  at 
which  the  specimen  failed.  The  residual  and  flexural  stress  profiles  for  96% 
alumina  quenched  in  silicone  oil  from  an  initial  temperature  of  1600°  C  are  shown 
in  Fig.  1.4a.  Compressive  surface  stresses  produced  by  tempering  have  reduced 
the  flexural  tensile  stress  from  724  MPa  to  330  MPa  (Fig.  1.4b). 

In  a  more  recent  study  by  Travitzky  et  al.  [48],  tiles  of  structural  alumina  at 
various  initial  temperatures  were  subjected  to  tempering  in  forced  air.   After 
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Fig.   1 .4    Stress  distribution  at  rupture  for  alumina  bars  quenched  in  silicone 
oil  from  an  initial  temperature  of  1600°C  [47]. 
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tempering,  a  diamond  saw  was  used  to  cut  the  tiles  to  produce  rectangular  bars 
for  mechanical  testing.  The  mean  value  of  flexural  strength  (three-point  bending) 
of  tempered  specimens  at  an  initial  temperature  of  140CTC  increased  to  430  MPa 
compared  with  the  225  MPa  reported  value  for  the  un tempered  samples.  This 
study  also  showed  that  the  flexural  strength  increased  significantly  with  initial 
tempering  temperature.  However,  at  initial  quenching  temperatures  above 
1400°C,  the  strength  no  longer  increased  (Fig.  1.5). 

Strengthening  of  Dental  Ceramics 

Background 

Dental  porcelain  frits  consist  of  fine  powders  that  are  formed  into  a  paste  by 
mixing  them  with  water.  The  desired  shape  is  then  fired  at  maturing 
temperatures  as  high  as  1300°C  to  produce  a  ceramic  restoration  with  color  and 
shade  similar  to  human  teeth.  Due  to  loss  of  water  and  a  sintering  process 
(fusion  of  the  porcelain  powder  to  form  a  continuous  mass),  the  porcelain  may 
shrink  during  firing  as  much  as  40  percent  by  volume.  The  firing  technique  has  a 
pronounced  effect  on  the  strength  and  quality  of  the  resultant  dental  restorations. 
The  major  mineral  constituent  of  a  high  temperature  maturing  porcelain  is 
feldspar  (Na20-Al203-6Si02  +K20-Al203-6Si02),  which  may  contain  as  much 
as  36%  alumina  (A1203)  by  volume  [49].  Similar  to  other  types  of  ceramic 
materials,  dental  porcelain  exhibits  brittle  fracture  and  because  of  surface 
irregularities,  it  has  low  tensile  strength. 

Shown  in  Table  1 .2  are  comparative  physical  properties  of  glass  and  two 
types  of  dental  ceramics  [5,  18,  50,  51]. 
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Table  1.2-  Physical  properties  of  glass  and  dental  ceramics. 


Specific 

Gravity 

(g/cm3) 


Thermal 
Diffusivity 
(mm2/sec) 


Modulus  of 

Elasticity 

(GPa)  ' 


Compressive 
Strength 
(MPa) 


Tensile 

Strength 

(MPa) 


Glass                     2.4 

0.57 

72.0 

383 

30-100 

Castable  Ceramic    2.7 

0.80 

70.3 

828 

34-69 

Porcelain                2.4 

0.64 

70.3 

172 

34-69 

To  overcome  the  limitations  associated  with  the  mechanical  properties  of 
porcelain  and  to  improve  the  structural  integrity  of  dental  crowns  and  bridges, 
Weinstein  et  al.  [52,53]  in  1962  developed  an  alloy  that  chemically  bonds  to 
porcelain.  Since  then,  construction  of  the  majority  of  dental  restorations  has  been 
based  on  porcelain- f used-to-metal  systems.  However,  structural  failure  of  these 
prostheses  as  well  as  all-ceramic  restorations  still  occur  despite  the  introduction 
of  more  compatible  systems  and  improvements  in  laboratory  techniques.  For  a 
glass-ceramic  molar  crown,  Moffa  [54]  has  reported  a  failure  rate  of  35.3% 
during  the  first  three  years  which  indicates  that  there  is  a  need  for  the 
development  of  new  dental  materials  with  improved  mechanical  and  thermal 
properties.  However,  at  the  present  time  it  seems  more  practical  to  focus  on  the 
development  of  better  laboratory  techniques  to  improve  the  design  of  dental 
restorations  using  the  existing  materials  rather  than  altering  their  chemical 
composition  to  produce  new  materials.  This  focus  on  existing  materials  is 
justified  because  all  conditions  required  for  biological  compatibility,  mechanical 
and  handling  properties,  esthetics  and  economics  should  be  satisfied  before  any 
new  material  can  be  used  clinically. 
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Influence  of  Crown  Design  on  Stress  Distribution 

Several  studies  have  analyzed  the  stress  distribution  in  dental  restorations 
due  to  mechanical  loading  in  order  to  improve  the  design  and  enhance  the 
structural  integrity  of  these  prostheses  [51,55-63].  In  our  previous  studies 
[51,58,59],  we  have  shown  that  the  stress  distributions  in  metal-ceramic  crown 
or  in  all-ceramic  crowns  are  not  significantly  influenced  by  small  changes  in 
metal  or  ceramic  thicknesses  at  some  locations.  Since  the  introduction  of  castable 
ceramics  in  1982,  several  studies  have  focused  on  optimizing  the  design  of  all- 
ceramic  crowns  based  on  finite  element  stress  analysis.  These  ceramics  are 
sometimes  preferred  to  metal-ceramic  crowns  for  restorations  of  anterior  teeth 
because  of  their  improved  aesthetic  potential. 

We  have  analyzed  [51]  the  stress  distribution  induced  by  applied  loads  in 
three  different  maxillary  ceramic  crown  designs  with  variable  thickness  of 
ceramic  in  the  incisal  region  (Fig.  1.6).  The  incisal  edges  of  ceramic  for  these 
designs  were  located  at  distances  of  1.0  mm  (case  I),  1.9  mm  (case  II)  and  4.0 
mm  (case  III)  from  the  incisal  edges  of  the  prepared  teeth.  Figs.  1.7  and  1.8 
show  the  influence  of  ceramic  thickness  and  loading  orientation  on  peak  tensile 
stress  in  the  ceramic  and  the  cement  layer.  The  locations  in  which  these  stresses 
were  obtained  are  also  shown  in  these  figures.  Comparison  of  the  maximum 
principal  stresses  induced  in  ceramic  and  dentin  for  cases  I  and  III  indicates  that 
the  increase  of  ceramic  thickness  from  1.0  to  4.0  mm  in  the  incisal  region  has  a 
negligible  effect  on  stress  distributions.  However,  the  influence  of  loading 
orientation  on  the  induced  stresses  is  significant.  A  200-newton  load  was 
distributed  at  the  lingual  region  along  vertical  and  horizontal  directions,  and  also 
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Fig.  1 .6    Two-dimensional  finite  element  model  of  a  prepared  central 

incisor  restored  with  each  of  three  ceramic  crown  designs  [51]. 
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Fig.   1 .7    Maximum  tensile  stress  (MPa)  for  crowns  with  ceramic  thickness 
of  1.0  mm  and  4.0  mm  produced  by  an  external  load  [51]. 
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Fig.   1.8    Maximum  compressive  stress  (MPa)  for  crowns  with  ceramic 
thickness  of  1.0  mm  and  4.0  mm  [51]. 
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at  angles  of  30°,  45°,  and  60°  to  the  horizontal  axis.  Maximum  tensile  stresses  in 
ceramic  and  cement  at  various  locations  in  the  crown  were  plotted  (Fig.  1 .9)  as  a 
function  of  the  angle  of  the  load  to  the  horizontal  axis.  When  the  load  was 
applied  along  a  horizontal  direction,  the  maximum  tensile  stresses  in  ceramic  and 
cement  were  288  MPa  and  47.4  MPa,  respectively.  These  values  were  much 
greater  than  the  published  tensile  strength  values  [59]  for  porcelain  and  cement. 
The  maximum  principal  stresses  in  ceramic  and  cement  increased  as  the  load 
orientation  changed  from  a  vertical  to  a  horizontal  orientation.  For  Dicor,  a 
tetrasilic  fluormica  glass-ceramic  (KMg2.5Si4O10F2),  Moffa  [54]  reported  a  failure 
rate  of  anterior  crowns  of  3.5%  during  the  first  three  years.  This  is  a  relatively 
small  rate  of  failure  which  indicates  that  such  a  horizontal  loading  condition 
analyzed  by  finite  element  methods  occurs  only  under  extraordinary  conditions. 

The  major  assumption  that  was  required  in  this  plane  stress  model  was  that 
all  vertical  planes  which  cross  through  the  crown  (sections  parallel  to  the  section 
shown  in  Fig.  1.6)  have  identical  stress  distributions  due  to  a  given  force.  In 
general,  a  plane-stress  distribution  corresponds  to  a  thin  plate  which  is  loaded  by 
forces  applied  at  the  plate  boundary,  parallel  to  the  plane  of  the  plate,  and 
distributed  uniformly  over  the  plate  thickness.  Such  an  ideal  case  does  not  exist 
in  a  ceramic  crown  under  the  influence  of  intraoral  forces.  A  three-dimensional 
model  would  yield  more  accurate  stress  values  compared  with  a  two-dimensional 
model.  However,  the  two-dimensional  model  employed  in  this  study  [51]  can  be 
used  for  a  relative  comparison  of  various  design  geometries. 

Hojjatie  and  Anusavice  [64]  have  developed  a  three-dimensional  finite 
element  model  to  determine  the  strength  of  posterior  ceramic  crowns  with  various 
designs  and  to  analyze  the  influence  of  loading  orientation  and  occlusal  thickness 
on  the  level  of  tensile  stresses  under  simulated  applied  loads.  The  results  of  the 
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Fig.    1 .9    Peak  values  of  principal  stresses  (Gmax)  in  ceramic  and  cement 
as  a  function  of  load  orientation  [51]. 
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finite  element  stress  analyses  of  ceramic  crowns  indicated  that  under  the  influence 
of  horizontally  applied  loads,  relatively  large  tensile  stresses  were  developed  in 
ceramic  despite  the  changes  in  the  crown  designs.  Therefore,  to  increase  the 
longevity  of  dental  restorations  under  intraoral  conditions,  a  combination  of 
strengthening  techniques  such  as  tempering  and  design  modifications  should 
also  be  pursued. 

Tempering  of  Dental  Ceramics 

The  potential  benefit  of  thermal  tempering  of  dental  ceramic  structures  was 
identified  several  years  ago  [49].  The  first  extensive  experimental  and  theoretical 
investigations  on  this  research  topic  have  only  recently  been  reported  [65-69]. 
These  studies  have  focused  on  analysis  of  the  effectiveness  of  this  technique 
when  it  is  used  for  strengthening  of  dental  porcelain.  In  one  of  these  studies 
[65],  it  was  reported  that  air  blasting  of  feldspathic  porcelain  can  significantly 
improve  the  resistance  to  crack  development  within  the  tempered  surface. 
However,  the  effect  of  transient  tensile  stresses  which  are  developed  at  the 
surface  of  porcelain  during  the  initial  stage  of  tempering  was  not  investigated.  If 
transient  tensile  stress  results  in  crack  initiation  and/or  propagation,  it  could 
negate  the  beneficial  effect  of  the  tempering  treatment. 

This  same  study  [65]  has  shown  that  the  mean  crack  size  induced  by  a 
microhardness  indenter  in  tempered  specimens  was  significantly  smaller  than 
those  for  the  conventionally  cooled  specimens.  This  suggests  that  the  surface 
compressive  stress  produced  by  tempering  decreases  the  growth  rate  of  surface 
cracks;  therefore,  this  procedure  can  improve  the  overall  strength  of  layered 
ceramics.  However,  several  variables  such  as  the  cooling  medium,  specimen 
geometry,  initial  tempering  temperature,  heat  transfer  coefficient,  and  other 
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thermal  properties  of  ceramic  materials  can  influence  the  effectiveness  of  the 
tempering  treatment. 

The  development  of  transient  tensile  stresses  in  ceramic  surfaces  during  the 
initial  stage  of  tempering  for  some  cooling  conditions  may  result  in  crack 
initiation  and  growth.  Also,  the  effect  of  residual  tensile  stresses  which  are 
produced  within  the  interior  of  tempered  ceramics  on  the  overall  strength  of 
restorations  must  be  determined.  Studies  on  tempering  of  glass  plates  indicate 
that  the  ratio  of  surface  compressive  stress  to  mid-plane  tensile  stress  increases 
with  the  initial  tempering  temperature  and  the  heat  transfer  coefficient  [20].  In 
addition,  the  initial  temperature  should  be  kept  low  enough  to  prevent  anatomical 
distortion  of  the  dental  prostheses  due  to  viscous  flow,  and  the  surface  heat 
transfer  coefficient  should  be  kept  small  enough  to  avoid  fracture  due  to  thermal 
shock.  The  transient  and  residual  stresses  can  be  calculated  based  on  viscoelastic 
analysis  or  elastic  analysis  using  analytical  approaches  and  finite  element 
methods. 

Stress  Analysis  bv  Viscoelasticity  and  Structural  Relaxation 

DeHoff  and  Anusavice  [66]  have  developed  an  analytical  model  to  calculate 
transient  and  residual  stress  distribution  in  dental  porcelain  plates  subjected  to 
cooling  rates  used  in  commercial  laboratories.  By  incorporating  linear 
viscoelasticity,  structural  relaxation  characteristic  of  porcelain,  and  an  assumed 
glass  transition  temperature  (Tg)  of  400°C  in  the  analytical  model,  they  have 
shown  that  for  symmetrical  cooling  of  a  2-mm-thick  porcelain  plate,  an  increase 
in  initial  cooling  temperature  from  700°C  to  800°C  does  not  change  the  stress 
distribution  significantly.  In  general,  relatively  low  compressive  stresses  were 
calculated  at  the  surfaces  of  porcelain.  At  an  initial  temperature  of  800°C  and  heat 
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transfer  coefficients  of  56.8,  113.6  and  170.3  W/m2-K,  the  stresses  were  26.4, 
17.6  and  8.5  MPa,  respectively.  The  viscoelastic  method  requires  identification 
of  a  glassy  and  a  liquid  state  for  porcelain  from  the  thermal  expansion/contraction 
data.  For  the  assumed  Tg  value  of  400°C,  the  following  linear  relations  were 
determined  from  the  experimental  dilatometry  data: 

Glassy  State    :  ag(t)  =  ao  +  aiT  at  Ta  <  T<  Tg  (1) 

Liquid  State     :  oq  (t)  =  a2  +  a3  T  at  Tg  <  T  <  Ti  (2) 

where  Ta  and  Ti  are  the  ambient  and  initial  cooling  temperatures  of  porcelain, 
respectively.  DeHoff  and  Anusavice  applied  the  procedures  reported  by  Lee  et 
al.  [33],  which  was  modified  by  Narayanaswamy  [28],  to  determine  the  transient 
stresses  in  a  biaxial  viscoelastic  slab  by  numerical  techniques.  A  summary  of  the 
procedures  used  in  the  study  by  DeHoff  and  Anusavice  [66]  is  given  in  Appendix 
A.  For  the  complete  description  of  the  mathematical  formulations  and  the 
numerical  technique,  one  should  refer  to  the  published  paper. 

Stress  Analysis  bv  Finite  Element  Method 

Transient  and  residual  tensile  stresses  produced  within  a  thermally  tempered 
ceramic  can  also  be  calculated  by  the  use  of  finite  element  analysis.  Some  general 
purpose  finite  element  programs  can  incorporate  the  theory  of  viscoelasticity  to 
determine  residual  stresses  in  ceramic  structures  with  complex  geometric  shapes. 
However,  the  finite  element  analyses  performed  in  this  study  are  based  on 
elasticity  theory  and  are  used  only  for  calculation  of  transient  stresses.  Although 
the  elastic  analysis  by  finite  element  method  assumes  that  the  viscoelastic  stresses 
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are  relaxed  during  cooling  through  the  glass  transition  region,  the  advantage  of 
this  method  over  the  viscoelastic  analysis  by  analytical  method  is  that  it  can  be 
used  for  a  single  layered  structure  as  well  as  multilayered  structures  with  any 
geometric  shape.  To  determine  the  transient  thermal  stresses  by  the  finite  element 
method,  first  the  unsteady-state  temperature  distributions  are  calculated  based  on 
an  energy  balance  within  the  model.  Then  the  calculated  temperatures  are 
incorporated  in  the  stress-strain  relations  corresponding  to  static  analysis.  A 
summary  of  the  finite  element  technique  for  temperature  and  stress  analyses  is 
given  in  Appendix  B. 

Objectives  of  the  Present  Study 

The  present  study  was  designed  to  optimize  the  tempering  process  for 
layered  porcelain  structures  as  a  function  of  the  thermal  characteristics  of  the 
ceramics  and  cooling  media  without  causing  transient  failure  or  inducing  high 
internal  residual  tensile  stresses.  The  specific  objectives  of  the  study  were  as 
follows: 

1 .  To  test  the  hypothesis  that  the  strength  of  feldspathic  porcelain  structures 
can  be  enhanced  by  thermal  tempering. 

2.  To  characterize  the  transient  temperature  and  stress  distribution  within 
bilayered  porcelain  disks  subjected  to  various  cooling  treatments. 

3.  To  test  the  hypothesis  that  thermal  tempering  of  dental  porcelain  in 

gaseous  cooling  media  is  more  effective  than  tempering  in  liquids. 

4.  To  analyze  the  influence  of  variables  such  as  cooling  method,  initial 

tempering  temperature,  thermal  contraction  mismatch,  and  specimen 
dimensions  on  the  transient  and  residual  stresses  produced  within 
ceramics  during  various  cooling  treatments. 


CHAPTER  n 

EXPERIMENTAL  METHOD  AND  ANALYSES 

Determination  of  Thermal  Data 

Measurement  of  Thermal  Expansion  Coefficients 

Five  specimens  in  each  group  of  two  opaque  feldspathic  porcelains  and  two 
body  feldspathic  porcelains  were  prepared  as  rectangular  bars  (5x5x50.8  mm) 
according  to  the  manufacturer's  recommendations  (J.F.  Jelenko  and  Co., 
Armonk,  NY)  and  were  used  to  measure  the  thermal  contraction  coefficient  (a)  in 
the  temperature  range  of  25-500°C.  The  measurements  were  conducted 
according  to  the  ASTM  standard  designation  E228-71  using  an  Orton  single 
push-rod  dilatometer  (Orton  Ceramic  Foundation,  Columbus,  OH).  The  percent 
contraction  versus  temperature  curves  for  each  specimen  were  obtained  by 
interfacing  the  dilatometer  to  a  recorder.  The  displacement  of  the  specimen 
because  of  temperature  change  was  transferred  through  a  fused  silica  pushrod  to 
a  linear  variable  differential  transducer  (LVDT).  The  LVDT  produced  a  voltage 
that  was  directly  proportional  to  the  displacement  of  the  specimen  under 
isothermal  conditions.  A  thermocouple,  which  was  in  contact  with  the  surface  of 
the  sample,  also  produced  a  voltage  proportional  to  the  specimen  temperature. 
The  voltages  produced  by  the  thermocouple  and  the  LVDT  were  amplified  and 
were  transferred  to  the  x  and  y  axes  of  the  chart  recorder,  respectively  (Fig.  2.1). 
The  a  values  in  the  temperature  range  between  25-500°C  were  determined  from 
the  contraction  curves  of  the  specimens  at  a  cooling  rate  of  3°C/min.  Shown  in 
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Fig.  2.1    Schematic  diagram  of  a  single  push-rod  dilatometer. 
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Fig.  2.2  is  the  percent  contraction  versus  temperature  for  two  opaque  porcelains 
(01  and  02)  and  two  body  porcelain  (Bl  and  B)  using  a  reference  temperature 
of  25 °C.    The  mean  values  of  the  linear  thermal  contraction  coefficients 
corresponding  to  five  specimens  in  each  group  for  the  opaque  porcelains  were 
OCoi  =  12.0  ppm/°C  and  a0i  =  14.2  ppm/°C.  The  corresponding  values  for  the 
body  porcelain  were  CCBi  =  11.0  ppm/°C  and  0CB1  =  13.5  ppm/°C.    These 
porcelains  were  employed  to  prepare  bilayered  circular  disk  specimens  with 
thermal  contraction  mismatches,  Aa  (oco-ocb),  of  +3.2  and  -1.5  ppm/°C.  Disk 
specimens  consisting  of  body  porcelain  only  (0Cb2=  13.5  ppm/"C)  were  used  as  a 
thermally  compatible  system  (Aa  =  0.0)  in  the  control  group  to  determine  the 
influence  of  thermal  contraction  incompatibility  on  transient  and  residual  stresses 
produced  by  various  cooling  methods.  These  porcelain  combinations  correspond 
to  the  largest  values  of  positive  and  negative  mismatch  types  studied  by 
Anusavice  et  al.  [65]. 

Preparation  of  Bilayered  Disk  Specimens 

One  set  of  bilayered  porcelain  disks  were  16  mm  in  diameter  and  2  mm  in 
thickness.  Each  disk  consisted  of  a  0.5  mm-thick  layer  of  opaque  and  a  1.5  mm- 
thick  layer  of  body  porcelain.  The  second  set  of  disks  were  32  mm  in  diameter 
and  6  mm  in  thickness.  Each  disk  consisted  of  1.5  mm  thick-layer  of  opaque 
porcelain  and  4.5  mm-thick  layer  of  body  porcelain.  In  both  disk  sizes  the 
thickness  ratios  of  body  to  opaque  porcelain  were  3  to  1  which  is  a  typical  value 
for  clinical  restorations.  The  opaque  and  body  porcelains  were  condensed  into 
slightly  oversized  cylindrical  brass  molds  to  compensate  for  sintering  shrinkage. 
After  condensing  the  specimen  with  a  Vibra  II  handpiece  (J.F.  Jelenko  and  Co., 
Armonk,  NY)  and  removing  the  excess  moisture  with  a  blotting  paper,  the 
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Fig.  2.2   Thermal  contraction  versus  temperature  profiles  of  opaque  and  body 
porcelains  using  a  reference  temperature  of  25 °C. 
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specimens  were  placed  near  the  open  door  of  a  preheated  Ney  Mark  IV  digital 
furnace  (J.  M.  Ney  Company,  Bloomfield,  CT)  for  10  minutes,  then  moved 
further  inside  the  furnace  at  an  initial  temperature  of  649°C  (1200T)  with  the 
furnace  door  closed.  The  furnace  temperature  was  then  raised  to  982°C 
(1800°F),  which  is  the  maturing  temperature  of  the  porcelain,  at  a  heating  rate  of 
20  to  55°C/min.  The  lower  heating  rate  was  applied  to  minimize  the  temperature 
gradient  within  the  6-mm-thick  disks  to  uniformly  fuse  the  porcelain  specimens 
and  to  prevent  thermal  shock  failure.  Firing  was  accomplished  in  a  vacuum  of  27 
in  Hg  to  minimize  surface  roughness  and  porosity  within  the  specimen.  After  the 
disks  were  heated  at  the  maximum  temperature  of  982" C  for  15  sec,  they  were 
taken  out  of  the  furnace  and  cooled  under  natural  (free)  convection  to  room 
temperature  (30°C).  Using  a  Buehler  polishing  wheel  (Buehler  Ltd.,  Lake  Bluff, 
EL),  each  porcelain  layer  was  ground  to  the  proper  dimensions  and  finished  with 
0.05  |im  alumina  abrasive. 

Temperature  Distribution  within  the  Dental  Furnace 

To  determine  the  steady-state  temperature  distribution  inside  the  furnace, 
four  chromel-alumel  thermocouples  (K  type)  were  calibrated  with  a  L  &  N 
reference  standard  thermocouple  (Leeds  &  Northrup  Company,  North  Wales, 
PA)  which  were  installed  at  various  locations  inside  the  furnace.  The 
thermocouples  were  interfaced  to  an  IBM  PC-AT  microcomputer  (IBM 
Corporation,  Boca  Raton,  FL)  through  a  pMEGA-4000  analog-to-digital 
converter  (Omega  Engineering,  Inc.,  Stamford,  CT)  which  was  plugged  into  the 
RS-232  serial  port  of  the  microcomputer.  The  furnace  temperature  was  set  at 
constant  temperatures  between  100  to  982°C  and  the  steady-state  temperature 
distribution  within  the  dental  furnace  for  each  set-point  temperature  was 
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measured  to  identify  the  regions  with  uniform  temperature.  Shown  in  Fig.  2.3 
are  the  locations  of  thermocouples  and  the  steady-state  temperature  distribution 
inside  the  heating  chamber  of  the  test  furnace  at  a  control  temperature  of  982°C. 
The  temperatures  at  the  right  side  of  the  heating  chamber  correspond  to  the  mean 
values  of  the  four  thermocouple  readings  at  each  position.  The  shaded  area 
which  is  identified  as  "hot  zone",  corresponds  to  the  region  inside  the  furnace 
with  minimum  thermal  gradients.  The  disk  specimens  were  placed  in  this  region 
for  firing. 

Measurement  of  Cooling  Profiles 

A  computer  program  in  BASIC  language  was  developed  to  record  the 
temperature  values  of  up  to  four  thermocouples  versus  any  arbitrary  time  interval 
which  was  controlled  by  using  the  IBM  PC- AT  internal  clock.  This  program 
used  the  commands  recognized  by  the  data  acquisition  unit  to  receive  temperature 
data  from  the  thermocouples  and  the  corresponding  time  values  from  the 
computer  clock.  These  data  were  stored  on  the  microcomputer  hard  disk  for 
plotting  and  further  analysis  of  the  cooling  profiles. 

To  determine  the  temperature  versus  time  profiles  for  each  cooling 
technique,  a  calibrated  thermocouple  was  installed  at  the  center  of  each  disk  on 
the  body  porcelain  surface  with  a  high  temperature  resistance  cement  (Omega 
Engineering,  Inc.,  Stamford,  CT)  and  the  thermocouple  was  interfaced  to  the 
microcomputer.  The  distance  between  the  center  of  the  thermocouple  tip  and  the 
disk  surface  was  0.5  mm.  Each  disk  was  supported  by  a  layer  of  refractory 
insulating  cloth  (Thermoz,  American  Dental  Supply,  Easton,  PA)  on  a  ceramic 
saggar  tray  and  was  placed  inside  the  furnace  at  the  "hot  zone"  region. 
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Fig.  2.3    Steady-state  temperature  distribution  in  the  test  furnace  at  a  control 
temperature  of  982°C. 
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Tempering  by  Gaseous  Cooling  Media 

To  determine  the  effectiveness  of  gaseous  cooling  media  and  cooling 
methods,  air,  nitrogen,  and  helium  were  employed  and  the  specimens  were 
cooled  by  free  convection  (FC),  tempering  by  vertical  flow  (TV)  and  tempering 
by  horizontal  flow  (TH)  from  an  initial  temperature  of  982°C.  Free  convection  in 
air  was  accomplished  by  bench  cooling  of  specimens  inside  the  test  laboratory. 
The  ambient  temperature  near  the  furnace  was  30°C.  Free  convective  cooling  in 
other  gases  was  performed  inside  a  relatively  large  environmental  chamber  which 
was  completely  sealed  from  the  outside  air.  The  temperature  inside  the 
environmental  chamber  during  each  test  was  monitored  by  two  digital 
thermocouples.  Shown  in  Fig.  2.4  is  a  schematic  diagram  of  the  tempering 
procedure  with  relative  positions  of  the  multiple-orifice  nozzle  and  the 
thermocouple  for  air  blasting  in  a  vertical  direction.  In  each  tempering  treatment 
the  gas  pressure  was  set  to  a  constant  value  of  0.34  MPa  (50  psi).  The  mean 
value  of  air  velocity  near  the  specimen  during  tempering  was  77  m/sec  as 
measured  by  a  hot  wire-anemometer  (Solomat  Corp.,  Stamford,  CT). 

Shown  in  Fig.  2.5  is  a  schematic  illustration  of  the  modes  of  heat  loss  from 
the  surface  during  the  tempering  procedure  using  horizontal  air  flow.  The  main 
sources  of  heat  transfer  during  cooling  treatment  with  gas  are  radiation  (Qr)  and 
convection  (Qc).  However,  heat  is  also  conducted  away  from  the  thermocouple 
junction  through  the  lead  wires.  This  heat  loss  can  cause  deviation  of  the 
thermocouple  output  from  the  actual  temperature  at  the  point  of  interest  [70].  To 
minimize  this  error,  the  lead  wires  near  the  thermocouple  junction  were  formed 
into  a  coil  (Fig.  2.5).  This  configuration  minimized  the  temperature  gradient 
between  the  thermocouple  sensor  and  the  lead  wire  and,  thus,  the  conduction  loss 
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Fig.  2.4    Schematic  illustration  of  vertical  air  flow  conditions  for  the  small 
disk  specimen. 
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Fig.  2.5    Schematic  mustration  of  horizontal  air  flow  conditions  and  modes 
of  heat  loss. 
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from  the  wire  to  the  surroundings.  For  cooling  treatment,  by  gaseous  media,  the 
bottom  surface  of  each  disk  was  placed  on  a  layer  of  insulating  material. 
Although  there  was  a  small  heat  loss  by  conduction  from  the  bottom  surface  to 
the  insulation,  for  the  determination  of  temperature  and  stress  distribution  by 
analytical  methods,  the  bottom  surface  was  assumed  to  be  perfectly  insulated. 

Tempering  by  Liquids 

To  analyze  the  effectiveness  of  tempering  in  liquid  media,  the  disks  were 
quenched  in  silicone  oils  (General  Electric  Company,  Waterford,  NY)  having 
kinematic  viscosities  of  50,  1000,  and  5000  centistokes.  These  oils  had  the  same 
specific  heat  (1,506  J/kg-K),  conductivity  (0.158  W/m-K),  and  density  (957 
kg/m3),  at  the  temperature  range  between  27°C  and  93°C.  The  silicone  oils  with 
variable  viscosities  were  chosen  based  on  the  findings  of  Kirchner  [38,47]  who 
quenched  alumina  in  silicone  oils  at  room  temperature  with  fluid  viscosities 
varying  from  5  to  12,500  centistokes.  He  reported  that  the  oils  with  lower 
viscosities  were  more  effective  than  those  of  higher  viscosities  because  of  their 
higher  heat  transfer  coefficient  which  produced  greater  compressive  stresses  at 
the  surface.  However,  quenching  media  with  very  large  heat  transfer 
coefficients,  such  as  water  will  often  promote  thermal  shock  failure.  Since 
quenching  the  bilayered  porcelain  disks  from  the  maturing  temperature  of  982°C 
in  room  temperature  silicone  oils  caused  thermal  shock  failure,  the  initial 
temperature  of  the  oil  was  increased  to  100°C  and  the  quenching  experiments 
were  conducted  from  specimens  at  initial  temperatures  of  850°C  or  lower. 
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Influence  of  Initial  Tempering  Temperature 

Travitzkey  et  al.  [48]  showed  that  initial  tempering  temperature  had  a 
significant  influence  on  the  flexural  strength  of  alumina  (Fig.  1 .4)  when  tempered 
by  an  air  jet.  In  the  present  study,  the  porcelain  disks  were  slow-cooled  (SC), 
fast-cooled  (FC)  and  tempered  (T)  in  air  or  in  silicone  oil  from  initial 
temperatures  of  650.  750,  850  and  982°C.  The  initial  temperature  of  air  was 
30°C  and  the  corresponding  value  for  silicone  oil  was  100°C.  The  effect  of  initial 
temperature  on  specimen  strength  using  silicone  oil  quenching  was  analyzed  via 
measurements  of  surface  cracks  induced  by  a  microhardness  indenter  and 
measurements  of  the  strength  by  means  of  a  biaxial  flexure  test.  Comparison  of 
the  data  was  performed  statistically  by  means  of  an  analysis  of  variance 
(ANOVA)  and  the  Tukey's  multiple  range  test  (HSD). 

Measurement  of  Heat  Transfer  Coefficients 

The  heat  loss  by  convection  from  the  disk  to  the  surrounding  for  each 
cooling  treatment  is  calculated  from: 

Qc=Ahc(Ts-Ta)  (3) 

where 

Qc  =  Convective  heat  loss  (W) 

A    =  Convective  surface  area  (m2) 

hc   =  Convective  heat  transfer  coefficient  (W/m2-K) 

Ts  =  Surface  temperature  (°C) 

Ta   =  Ambient  temperature  (°C) 

Equation  (3)  has  been  simplified  by  incorporating  all  of  the  complex 
convection  characteristics  in  the  coefficient  hc,  which  is  generally  a  function  of 
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fluid  flow,  thermal  properties  of  the  fluid  medium,  and  geometry  of  the  system 
[71].  A  common  method  to  determine  hc  is  to  combine  dimensional  analysis 
with  experimental  data.  It  is  important  to  emphasize  that  the  heat  transfer 
coefficient  is  influenced  by  the  shape  and  the  dimensions  of  the  body  under 
convection  and  not  by  its  material  properties.  According  to  the  principle  of 
similarity  (model  law),  the  behavior  of  two  systems  will  be  similar  if  the  ratios  of 
their  linear  dimensions  are  the  same.  Therefore,  for  the  experimental 
determination  of  the  heat  transfer  coefficient,  this  principle  was  applied  and  a 
copper  disk  with  the  same  diameter  as  the  porcelain  disk  was  used.  This 
simplified  both  the  experiments  and  the  calculations,  because  the  same  copper 
disk  was  used  to  measure  experimentally  the  cooling  profiles  for  various  cooling 
conditions.  Also,  because  of  the  large  thermal  conductivity  of  copper  and  the 
small  thickness  of  the  disk,  the  temperature  gradient  within  the  copper  was 
assumed  to  be  negligible  and  the  heat  transfer  coefficient  was  determined  by  a 
lumped  analysis  according  to  the  following  method: 

d  Ts 

Rate  of  energy  loss  by  copper  disk  =  Q  =  p  V  Cp  — gp  (4) 

where  p  ,  V  and  Cp  are  the  density,  volume  and  specific  heat  of  the  copper  disk, 
respectively.  At  any  instant  of  time,  the  rate  of  energy  loss  by  the  copper  disk  is 
the  same  as  the  convection  loss  from  the  disk  to  the  surrounding  fluid: 


A  T"1 

p  V  Cp-^1  =  Ahc(Ts-Ta)  (5) 

Using  the  initial  condition  T(0)  =  Tj ,  one  obtains  : 

.         -p  Cp  Lc  .     ,  Ts-Ta  ,,v 

hC=^1 In  (   Ti_Ta)  (6) 
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where  Lc  =  V/A  =  characteristic  length  of  the  disk,  t  =  time  and  A  =convective 
surface  area.    To  eliminate  Tj  ,  two  locations   (Ti  ,  ti)  and  (T2 ,  t/2)  on  the 
experimentally  determined  cooling  curve  were  used  [72].   Then  from  (6)  one 
obtains  : 

-  p  CP  Lc  T3  -Ta  -p  CPLC  T2-Ta 

hc= ij—  ln(T._Ta)  =  t2  ln(T-f^)  (7) 


By  eliminating  Ti,  it  can  be  reduced  to 


"c=    Vi   cMf^>  (8) 


The  condition  for  a  valid  lumped  analysis  [73]  is  that    Bj  =  — tt —  <  0.1 

where  K  is  the  thermal  conductivity  of  the  copper  disk  and  Bj  is  the  Biot 
number  .  For  an  unsteady  problem,  first  Bj  is  evaluated. 
Equations  (6)  and  (8)  are  commonly  used  transient  methods  for  determining  the 
heat  transfer  coefficient  by  experimentally  measuring  the  slope  of  the  cooling 
profile.  However,  this  method  is  associated  with  an  error  due  to  differentiation 
of  the  cooling  profile  [72].  In  this  study,  the  heat  transfer  coefficient  was  also 
calculated  from  other  analytical  and  numerical  approaches  and  the  values  were 
checked  by  comparing  the  calculated  and  experimentally  measured  cooling 
profiles. 

Analysis  of  Induced  Cracks  by  Microhardness  Indentation  Test 

The  microhardness  indentation  method  which  is  a  simple  technique  for  rapid 
evaluation  of  fracture  toughness  was  used  first  by  Evans  and  Charles  in  1976 
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[74].  They  concluded  that  the  Vickers  indentation  technique  should  enable  the 
calculation  of  fracture  toughness  data  within  an  accuracy  of  approximately  10 
percent  if  the  Young's  modulus  of  the  material  is  known.  The  fracture  toughness 
value  provides  a  measure  of  resistance  to  crack  initiation  and  extension  in  brittle 
materials.  The  fracture  toughness,  Kc ,  was  calculated  from  the  following 
relations  [75-77]  : 

Kc  =  XP/(CSc)3/2  W 

X  =  1.6xlCT8(E/H)1/2  (10) 

Where: 

Kc  =  fracture  toughness  (MPa-m1/2) 

%     =  dimensionless  indenter  constant 

P      =  indentation  load  =  4.9  N 

Csc  =  crack  size  for  slowly  cooled  disks  (m) 

E     =  elastic  modulus  of  porcelain  =  70.3  GPa 

H    =  hardness  of  porcelain  (GPa)  =  0.47  P/ap2 

ap    =  Projected  length  of  the  indenter  half-diagonal  indentation 

(  ap  =  21  x  10"6  m  for  a  load  of  4.9  N) 
For  determination  of  Kc,  the  surface  of  the  porcelain  disk  should  contain 
no  significant  residual  stresses  prior  to  indentation.  Therefore,  in  equation  (9), 
the  crack  size  measured  for  the  slowly  cooled  specimens  was  used.  In  this 
equation,  the  ratio  P/(CSc)3/2  was  approximately  independent  of  the  applied  load. 
Although  porosity  is  known  to  reduce  the  elastic  modulus  of  porcelain,  no 
correction  was  made  because  data  on  the  variation  of  E  with  porosity  are  not 

available  [73]. 

Traditionally,  residual  stresses  produced  at  the  surface  of  a  thermally  treated 
glass  have  been  determined  by  the  use  of  the  photoelastic  birefringence 
technique.  However,  this  optical  method  cannot  be  used  in  ceramics  because  of 
their  opacity.  One  of  the  most  convenient  methods  of  stress  determination  is  the 
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measurement  of  the  surface  cracks  induced  by  a  microhardness  indenter  [79-82]. 
Comparative  evaluation  of  crack  size  on  surfaces  of  tempered  and  untempered 
disks  will  provide  information  on  the  residual  compressive  stresses  which  are 
produced  because  of  the  tempering  treatment. 

The  specimen  surfaces  were  cleaned,  polished,  and  then  positioned  under 
the  diamond  pyramid  indenter  on  an  adjustable  mechanical  stage  corresponding  to 
a  Vickers  microhardness  indenter  (Tukon  Wilson,  Binghampton,  NY).  An 
indentation  load  of  4.9  N  (500  grams)  was  applied  through  the  mechanical 
control  of  a  dash  pot,  which  eliminates  shock  and  impact  and  vibration.  The 
dash  pot  causes  a  delay  on  the  movement  of  the  indenter  and  load  assembly  to 
allow  absorption  of  machine  disturbances  before  the  indenter  meets  the  specimen 
surface.  After  a  short  contact  between  the  indenter  and  the  specimen,  the  dash 
pot  automatically  raises  the  indenter  above  the  specimen.  By  using  the 
mechanical  stage,  the  specimens  were  placed  under  the  microscope  to  measure 
the  size  of  the  induced  indentation  crack  by  a  micrometer.  The  crack  sizes  were 
measured  within  60  seconds  after  indentation  to  minimize  error  due  to  crack 
growth  in  the  presence  of  environmental  moisture.  Twenty  indentations  were 
applied  on  the  surface  of  body  porcelain  in  each  disk  along  two  perpendicular 
diametral  lines.  Each  line  contained  10  indentations  at  the  surface  of  body 
porcelain  (Fig.  2.6).  The  crack  sizes  which  were  measured  for  each  of  the  three 
silicone  fluids,  three  mismatch  levels,  and  three  initial  cooling  temperatures  were 
compared  with  the  corresponding  crack  measurements  for  the  slow-cooled  (SC), 
fast-cooled  (FC)  and  tempered  (T)  specimens  in  air  from  an  initial  temperature  of 
982°  C  [Anusavice  et  al.,  65].  Analysis  of  variance  (ANOVA)  was  performed  to 
determine  if  the  crack  sizes  were  significantly  influenced  by  cooling  media, 
contraction  mismatch  and  initial  temperature.  If  the  influence  of  each  variable  or 


Fig.  2.6    Locations  and  directions  of  specimen  and  indentations  applied  at 
the  surface  of  the  body  porcelain. 
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their  combined  effects  was  significant,  a  Tukey's  multiple  contrast  analysis 
(HSD)  was  employed  at  a  significant  level  of  5%  to  determine  if  the  mean  values 
were  significantly  different.  The  mean  values  of  crack  sizes  for  tempered  (T), 
fast-cooled  (FC)  and  oil  quenched  (OQ)  specimen  were  used  to  calculate  the 
residual  surface  stresses  from  the  following  equation  which  was  developed  by 
Marshall  and  Lawn  [81]: 

a-  f  (g)"2  P-(^)3fl]  (11) 

,where  C$c  is  the  crack  size  for  slowly  cooled,  C  is  the  crack  size  for  the  fast- 
cooled  or  tempered  specimens,  and  Kc  is  the  fracture  toughness  calculated  from 
equation  (9). 

Biaxial  Flexure  Test 

Background 

The  biaxial  flexure  test  was  developed  by  Wachtman  et  al.  in  1972  [83]  for 
determining  the  strength  of  thin,  circular  disks  of  brittle  materials.  In  this  method 
the  specimen  is  supported  on  three  equally  spaced  balls  and  a  load  is  applied  to 
the  center  of  the  specimen  by  a  flat  piston.  Unlike  the  three-point  or  four-point 
bending  (uniaxial  tests),  the  strength  value  calculated  by  this  method  is  unaffected 
by  the  edge  condition  of  the  specimen.  In  the  biaxial  test,  the  three  points  of  the 
supports  are  located  near  the  periphery  and  equidistant  from  the  center  of  the 
disk.  Therefore,  the  area  of  the  maximum  tensile  stress  falls  at  the  center  of  the 
lower  face  of  the  disk  and  the  strength  should  be  independent  of  the  edge 
conditions.  The  breaking  load,  the  dimensions  and  elastic  constants  of  the 
specimen,  and  the  radii  of  the  support  and  the  load  are  used  to  compute  the 
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maximum  tensile  stress.  Wachtman  et  al.  [83]  employed  this  method  to 
determine  the  strength  of  eight  types  of  alumina  substrate  and  they  concluded  that 
if  the  mean  biaxial  flexure  value  of  five  specimens  is  used  as  a  measure  of 
strength,  a  coefficient  of  variation  of  about  7  percent  can  be  expected  for  similar 
materials.  This  is  acceptable  for  most  design  considerations. 

Flexural  Testing  of  Porcelain  Disks 

Based  on  the  results  from  the  indentation  test  for  the  five  mismatch  group 
reported  by  Anusavice  et  al.  [65],  disk  specimens  with  thermal  contraction 
differences  of  +3.2,  0.0  and  -1.5  ppm/°C  were  selected.  Six  disks,  16  mm  in 
diameter  and  2  mm  in  thickness  were  prepared  for  each  mismatch  group  and  each 
cooling  condition.  Also,  to  determine  the  influence  of  initial  temperature  on 
flexural  strength,  the  compatible  disks  (Ace  =  0.0)  were  treated  by  slow  cooling 
(SC),  fast  cooling  (FC)  and  tempering  (T)  methods,  each  from  initial 
temperatures  of  650,  750,  and  850°C.  Six  specimens  were  prepared  for  each 
group  and  each  disk  was  subjected  to  biaxial  flexure  testing  under  water  at  37°C. 
Shown  in  Fig.  2.7  is  the  schematic  drawing  of  the  biaxial  test  apparatus.  Each 
disk  was  placed  upon  three  equidistant  stainless  steel  balls  (1.6  mm  in  diameter) 
which  lie  upon  the  perimeter  of  a  circle  10  mm  in  diameter.  The  disks  were 
loaded  in  water  at  37°C  at  the  untempered  surface  by  means  of  a  steel  piston 
(with  a  3. 2- mm  diameter  flat  ground  along  the  contacting  surface)  at  a  loading 
rate  of  0.5  mm/min  until  fracture  occurred.  Since  it  is  expected  that  this  loading 
condition  causes  tensile  failure  at  the  surface  opposite  to  the  applied  load,  the 
tempered  surface  of  each  disk  was  placed  at  the  bottom  and  a  compressive  force 
was  applied  to  the  center  of  the  top  surface.  This  mode  of  flexure  caused  tensile 
failure  to  initiate  at  the  center  of  the  lower  surface. 


. 


Fig.  2.7    Biaxial  flexure  test  apparatus. 
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The  failure  stress,  Gf,  at  the  center  of  the  tempered  surface  was  calculated 

using  the  following  equations  developed  by  Marshall  [84]  for  the  piston-on- 

three-ball  flexure  test  described  by  Wachtman  et  al.  [83]: 

GF  =  XPf/a2  (12) 

X  =  (3/4:1)  [2(l+v)  Ln  (Rs  /  R0*)  +  (1-v)  (2RS2-R0*2)  /  2R2+(l+v)  (13) 

Ro*  =  (1.6Ro2  +  a2)i/2_0.675a  (14) 

where 

Gf  =  flexural  strength 

Pf  =  apphed  load  at  failure 

v    =  Poisson's  ratio  (0.28  for  dental  ceramics) 

Rs  =  radius  of  the  support  circle 

Ro  =  radius  of  the  piston  at  the  contacting  surface 

R    =  radius  of  disk 

a    =  thickness  of  the  disk 

The  failure  stress  data  were  analyzed  by  a  two-factor  ANOVA  and  Tukey's 
multiple  range  test  (HSD)  to  determine  the  significance  of  cooling  condition  and 
contraction  difference  on  biaxial  flexure  strength. 


CHAPTER  HI 

ANALYTICAL  AND  NUMERICAL  METHODS 

Calculation  of  Heat  Transfer  Coefficients 


Free  Convective  Cooling 
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The  convective  heat  transfer  coefficient  (hc)  as  a  function  of  specimen 
surface  temperature  for  each  cooling  method  was  calculated  from  the  appropriate 
empirical  correlations  available  for  each  cooling  medium  and  each  flow  direction. 
The  results  were  used  in  a  finite  element  thermal  analysis  to  determine  the 
temperature  distribution  within  each  specimen.  Preliminary  finite  element 
analyses  were  performed  using  hc  as  a  function  of  temperature.  However,  for 
more  simplification,  subsequent  analyses  were  performed  by  assuming  an 
average  heat  transfer  coefficient  for  each  condition.  The  hc  value  for  each 
cooling  method  was  calculated  from  the  following  relation  : 

h       NuKm  ™ 

hc  =  ^— -  (15) 

,where  Nu  =  Nusselt  number  =  ratio  of  conductive  thermal  resistance  to  the 
convective  thermal  resistance  of  the  cooling  medium 

Km  =  thermal  conductivity  of  the  cooling  medium  (W/m2-K) 
Lc   =  characteristic  length  of  the  specimen  (m) . 

For  the  free  convection  (FC)  case  the  value  of  Nu  was  obtained  from  the 

following  equations  [85]: 

Nu  =  0.54  Ra1/4  (16) 
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2g  (Ts  -  Ta)  Lc3 

Ra=    yK(Ts  +  Ta)  (17) 


where  Ra  =  Rayleigh  number 

g  =  acceleration  due  to  gravity  =  9.8  m/s2 

Ts  =  specimen  surface  temperature  at  the  center  of  body  porcelain  (°C) 

Ta  =  ambient  temperature  =  30°C 

Lc  =  characteristic  length  (m) 

y  =  kinematic  viscosity  of  the  cooling  medium  (m  /s) 

k  =  thermal  diffusivity  of  the  cooling  medium  (m2/s) 

Although  equation  (16)  has  been  developed  based  on  the  experimental  data 
for  a  relatively  large  plate  subjected  to  convective  cooling,  several  studies  have 
obtained  satisfactory  results  using  this  equation  or  a  similar  equation  for  disks  or 
plates  with  relatively  small  dimensions  [86-88].  In  the  present  study,  the 
calculated  heat  transfer  coefficient  from  equation  (16)  was  used  as  an  initial  value 
in  the  thermal  models  to  determine  the  temperature  distribution  numerically. 
Then  this  calculated  value  was  modified  to  produce  calculated  cooling  profiles 
comparable  with  the  experimental  data.  Based  on  the  gradient  of  the 
experimentally  determined  temperature  versus  time  profile,  various  surface 
temperatures  (Ts)  in  the  range  of  30-982°C  were  used  in  equation  (17).  For  each 
constant  value  of  Ts  ,  the  thermal  properties  of  the  cooling  media  were  evaluated 
at  the  following  temperature  [85]: 

Tftn=1^  +273.2  (18) 

where  Tfm  =  film  temperature  (°K).  For  each  constant  value  of  the  specimen 
surface  temperature  (Ts),  a  value  of  Ra  and  Nu  was  calculated  from  equations 
(17)  and  (16),  respectively.  Using  the  appropriate  values  of  Nu  and  Lc  in 
equation  (15),  the  convective  heat  transfer  coefficient  as  a  function  of  specimen 
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surface  temperature  was  determined  for  each  specimen  dimension  and  cooling 
medium. 

Characteristic  Length 

Many  discrepancies  exist  regarding  the  definition  of  characteristic  length  for 
bodies  subjected  to  natural  convection.  King  [89]  has  reported  the  following 
relation  for  the  characteristic  length  of  a  body  whose  horizontal  dimension  is  Lh 
and  vertical  dimension  is  Lv  : 

1   _  J_  +  1  (19) 

Lc       Lh      Lv 

He  indicated  that  the  results  of  this  equation  were  very  satisfactory  for  bodies 

such  as  small  short  cylinders,  planes  and  spheres.  This  relation  for  a  porcelain 

disk,  16  mm  in  diameter  (Lh)  and  2  mm  in  thickness  (Lv)  gives  a  characteristic 

length  of  1.78  mm.  Equation  (19)  was  quoted  for  many  years  in  heat  transfer 

textbooks.  However,  in  1983,  Sparrow  and  Ansari  [90]  conducted  experimental 

studies  on  free  convective  heat  transfer  for  multidimensional  bodies.    They 

concluded  that  King's  relation  (eq.  19)  was  inaccurate  and  its  citation  in  the 

modern  literature  [91-94]  has  been  propagated  from  older  textbooks  without 

challenge.  Lienhard  [95]  recommended  that  the  travel  length  of  a  fluid  in  the 

boundary  layer  be  used  as  a  characteristic  dimension.   For  the  porcelain  disk 

cooled  by  free  convection  from  the  top  and  sides,  and  insulated  from  the  bottom 

surface,  using  Lh  =16  mm  and  Lv  =  4  mm,  this  equation  gives  a  Lc  value  of  20 

mm.  Goldstein  et  al.  [96],  reported  the  following  commonly  used  relation  : 

LC  =  A/PS  (20) 
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,where  A  is  the  surface  area  and  Ps  is  the  perimeter.  Equation  (20),  which  was 
confirmed  by  Lloyd  and  Moran  [97],  is  expected  to  enable  horizontal  plates  as 
well  as  nonsymmetrical  surfaces  to  exhibit  a  common  correlation.  Using  this 
relation  the  value  of  Lc  for  the  disk  is  4  mm.  However,  Kreith  [71]  has  applied  a 
value  of  0.9  D  (D  =  disk  diameter),  for  the  characteristic  length  of  a  horizontal 
disk  which  gives  an  Lc  value  of  14.4  mm  for  the  porcelain  disk.  Several  authors 
have  used  the  Lc  value  to  be  the  same  as  the  disk  diameter  [86,87,90,96]  which 
was  16  mm  for  the  bilayered  disk.  In  view  of  these  discrepancies,  it  is  unlikely 
that  the  selection  of  a  single  equation  would  be  sufficient  to  determine  the 
characteristic  length  of  a  multidimensional  body.  Most  of  the  correlation 
formulas  available  for  free  convection  are  based  on  the  experimental  studies  of 
relatively  large  bodies.  However,  the  characteristic  length  of  D  gave  sufficiendy 
accurate  results  for  small  disks  with  variable  dimensions  in  the  range  of  1  to  20 
mm  [87].  For  the  analysis  of  the  convective  cooling  conditions  in  the  present 
study,  the  heat  transfer  coefficients  were  calculated  using  the  characteristic  length 
D  and  the  results  were  in  agreement  with  the  published  values  [73]. 

Tempering  by  Vertical  Row 

For  tempering  by  vertical  flow  (TV),  the  Nusselt  number  was  calculated 

from  [98]: 

Nu  =  0.93  Re0-5  Pr  °-4  (21) 

where   Re  =  Reynolds  number  based  on  radius  =  V  R/v 
R  =  Radius  of  the  porcelain  disk  (m) 
Pr  =  Prandtl  number  =  \i  Cp/K 
ju.    =  Dynamic  viscosity  of  the  cooling  media  (Pa-s) 
Cp  =  Constant  pressure  specific  heat  of  the  cooling  media  (kJ/kg-K) 
K  =  Thermal  conductivity  of  the  porcelain  disk  (W/m-K). 
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As  in  the  FC  case,  for  constant  values  of  surface  temperature,  the  thermal 
properties  of  cooling  media  were  evaluated  at  the  film  temperature  (18)  and  the 
corresponding  values  of  Nu  calculated  from  equation  (21)  were  employed  in  (15) 
using  the  specimen's  radius  as  the  characteristic  length. 

Tempering  by  Horizontal  Flow 

The  convective  heat  transfer  coefficient  for  tempering  by  a  horizontal  flow 
(TH)  was  approximated  by  the  equation  corresponding  to  laminar  flow  over  a  flat 

plate  [73]: 

Nu=  0.665  Re1/2Pr1/3  (22) 

where  Re  =  Reynolds  number  based  on  the  disk  diameter  =  V  D/v 

The  values  of  Nu  calculated  for  the  constant  surface  temperature  were 
applied  in  equation  (15),  using  D  as  the  characteristic  length,  and  the 
corresponding  hc  values  were  incorporated  in  equation  (3)  to  determine  the 
transient  heat  loss  due  to  convection. 

Radiative  Heat  Transfer  Coefficient 

If  the  cooling  medium  is  an  infrared  opaque  liquid,  all  of  the  heat  dissipated 
from  the  porcelain  disk  occurs  essentially  by  convection,  but  if  it  is  a  gas,  the 
exchange  of  radiant  heat  between  the  disk  and  the  surrounding  surfaces  must  be 
taken  into  account  in  computing  the  total  heat  loss.  It  is  convenient  to  include 
both  radiation  and  convection  in  expressing  the  total  heat  loss  from  the  disk 
surfaces.  This  is  accomplished  by  introducing  a  radiative  heat  transfer 
coefficient.  For  a  constant  surface  temperature  corresponding  to  an  instantaneous 
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time  of  cooling  in  free  convection  (FC)  and  Tempering  (T)  cases,  a  radiative  heat 
transfer  coefficient  is  defined  according  to  the  following  equation  [99]: 

.        e  QSB  (Ts4-Ta4)  .... 

hr  =  ds-Ta)  {Zi) 

where  e  =  emissivity  of  porcelain  =  0.99 

oSB  =  Stefan-Boltzmann  constant  =  5.669  x  10~8  W/m2-K4  . 

The  total  value  of  the  heat  transfer  coefficient  (ht)  due  to  combined  effects  of 
convection  and  radiation  losses  for  each  cooling  condition  was  obtained  from  : 


ht  =  hc  +  hr  (24) 

Closed-Form  Solutions 
Two-Dimensional  Model  (2-D) 

For  the  calculation  of  transient  and  residual  stresses  within  a  disk  during 
each  cooling  treatment,  one  needs  to  determine  initially  the  corresponding 
transient  temperature  distributions.  Since  the  disk  is  subjected  to  approximately 
uniform  cooling,  an  axisymmetric  thermal  model  can  be  used  in  which  the  heat  is 
transferred  from  the  disk  to  the  cooling  media  by  convection  and  radiation.  For  a 
further  simplification  of  the  analysis,  the  radiation  was  taken  into  account  via  a 
modified  convective  heat  transfer  coefficient.  For  gaseous  cooling  media,  the 
bottom  surface  of  the  disk  was  assumed  to  be  perfectly  insulated  (Fig.  2.4).  But 
for  quenching  in  liquids,  the  heat  is  also  transferred  from  the  bottom  surface. 
Under  this  condition  the  temperature  distribution  is  symmetric  with  respect  to  the 
midplane  of  the  disk  which  can  be  treated  as  an  insulating  boundary.  Therefore, 
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The  calculated  temperature  versus  time  values  obtained  from  the  one- 
dimensional  model  were  incorporated  in  the  viscoelastic  model  developed  by 
DeHoff  and  Anusavice  [66]  to  determine  the  tempering  (residual)  stresses  within 
each  disk.  However,  for  calculation  of  the  transient  stresses,  an  axisymmetric 
finite  element  model  was  developed.  Although  the  viscoelastic  model  considers 
the  stress  relaxation  behavior  of  porcelain,  it  can  be  applied  only  for  the  single 
layered  disks.  Also  this  model  is  based  on  the  assumption  that  temperature 
gradients  along  directions  parallel  to  the  surface  are  negligible,  while  in  the 
axisymmetric  finite  element  model  these  temperature  gradients  were  taken  into 
account. 

The  equation  governing  an  axisymmetric  temperature  distribution  T  (r,z,t) 
in  the  disk  shown  in  Fig.  2.5  is  : 


Tt^Trr  +  'pTr  +  Tzz)  (25) 


which  is  subjected  to  the  following  initial  and  boundary  conditions: 

T  (r,z,0)  =  Ti  (26) 

Tr  (0,z,t)  =  0 
-KT2  (r,a,t)  =  h  [T(r,a,t)  -Ta] 
-KTr(R,z,t)  =  h[T(R,z,t)-Ta] 
and  Tz(r,0,t)  =  0 
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the  same  solution  was  applied  for  gaseous  and  liquid  media,  except  that  in  the 
liquids,  only  half  of  the  disk  thickness  was  introduced  into  the  equations.  Since 
this  is  an  axisymmetric  problem  with  well  defined  initial  and  boundary 
conditions,  a  closed-form  solution  exists. 

In  this  study,  first  the  exact  solution  of  the  axisymmetric  model  was 
obtained.  Then  the  problem  was  simplified  by  the  use  of  a  one-dimensional 
model  (Fig.  3.1).  The  results  from  these  two  models  were  compared  with  the 
experimental  data  for  some  cooling  conditions. 


2-D 


Fig.  3.1    Simplification  of  the  heat  transfer  analysis  from  a  two-dimensional 
model  to  a  one-dimensional  analytical  model. 
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where  k  is  thermal  diffusivity  arid  K  is  thermal  conductivity  of  porcelain.  An 
average  heat  transfer  coefficient  (h)  is  assumed  for  the  top  and  side  surfaces.  To 
obtain  the  non-dimensional  form  of  this  equation  let: 


H=S       ,        C=a       -        T  =  ^R         and       9=T~TZ  (27) 

One  then  finds: 


ef-)ex  =  enTi+  jj  &n  +  <lh 2  eCC  <28) 


with 


e  (ti,C0)  =  i  (29) 

©n  C0.Ct)  =  ° 
en(i,^)+Bne(i,CsT)=o 

0^(r|,l,T)  +  Bi2e(ri)l,x)  =  O 
and  6^  01.0/0  =  0 


where  Bii  =    — g—       and  Bi2  =  — jr-  (30) 


are  the  Biot  numbers  based  on  disk  radius  and  disk  thickness,  respectively. 
Using  the  the  method  of  separation  of  variables  to  solve  the  PDE  (28),  one 
obtains: 


0(ti£,t)  =  A(ti).¥(C)-Q(t)  (31) 
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Substituting  (31)  into  (28)  and  applying  the  symmetric  (r  =  0)  and  the  insulation 
(z  -  0)  boundary  conditions  yields: 


0  Cn,C,x)  =  I     I  {  Am,n  exp  [-  |  (Ym2  +  ^n2)  T]  [J0(^n  n)l  cos  (     ~^-  Q} 
m=l  n=l 

(32) 


where    Ym  and    ^n    are  characteristics  of  the  transcendental  equations 
corresponding  to  the  top  and  right  side  boundary  conditions,  respectively.  J0  and 
Jl  are  the  Bessel  functions  of  the  first  kind  of  order  zero  and  order  one, 
respectively  (Fig.  3.2).  The  initial  condition  and  the  following  equations  [100] 
were  used  to  find  a  relation  for   Am,n  : 


1  i 

J[TlJo($nTl)]dn     =£   Ji($n)  (33) 


1    r,  2 


fft  Jo2(^nri)]dTl     =  |  [Jo   (^n)+  Ji    (^n)l 
Finally  one  obtains : 


(34) 


0  (ri,C,x)  =  4   I    Z  {    exp  [-  I  (Ym2  +  ^n2)  X]  }  (Cn)  (Cm)  (35) 

m=l  n=l 


.  .  ~  Jq  (sn  X|) ,_,. 

with  Cn   =       r         2  7£   v  ,  T  2  7*   x,  ^o) 

^n  [Jo    (Sn)  +  Jl     (qn)J 
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Fig.  3.2    Bessel  function  of  the  first  kind  of  order  zero  and  order  one. 
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sin  (-g-ym)  (cos  ("j^Ym  0) 

O  Q 

[  ^Ym  +  sin  (^Ym)  (cos  (^"Ym  ))] 


and  Cm  =        - ; (37) 


where  Yn   and  £,n  are  characteristics  of  the  following  transcendental  equations: 


.  /  a        \  a  Yn  /OON 

COt(RYn)=^Th  (38) 


Bil    Jo  (  £0  =    ^n  Jl  &0  •  (39) 

The  characteristics  of  (38)  can  be  easily  determined  by  a  numerical  or  graphical 
technique.  The  solution  of  (39)  is  obtained  by  the  Newton-Raphson  method  as 
follows: 
Let: 

F(  5n)=Bii  J0(§n)-   ^n  Jl  fei)  (40) 

then  by  using  the  relation  : 


5  [Jl(x))]=J0(x)-  £Jl(x)  (41) 


one  obtains: 


4[F(^n)]  =  F(§n)  =  -[^n  Jo(^n)  +  Bii  Ji(£n)] 


where  via  Newton-Raphson  one  obtains: 


(42) 


&,  (1+1)=  ^(D-  yd^)  c43) 
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For  a  known  value  of  Bii  and  an  initial  guess  for  £,n  ,  initial  values  for  F(£,n) 
and  F'  (c,n)  are  found  by  using  the  Integrated  Mathematical  and  Statistical 
Libraries  (IMSL),  subroutine  programs.  The  solution  by  the  Newton  -  Raphson 
method  (43)  is  strongly  influenced  by  the  initial  value  £,n  (I)-  Therefore,  to 
determine  the  value  of  E,n,  some  approximate  characteristic  values  for  (39)  are 
first  obtained  by  a  graphical  methods.  Then  each  approximate  value  is  used  as  an 
initial  value  in  (43).  Once  the  solutions  of  (38)  and  (39)  are  obtained,  they  are 
used  in  (36),  (37)  and  (35)  to  determine  the  transient  temperature  distributions  at 
the  locations  of  interest. 

A  more  simplified  one-dimensional  thermal  model  can  be  developed  by 
assuming  that  the  heat  transfer  from  the  side  of  the  disk  is  negligible.  It  is 
predicted  that  for  large  values  of  heat  transfer  coefficient  (tempering),  the 
difference  between  the  results  from  the  1-D  and  the  2-D  models  is  negligibly 
small  because,  under  this  condition,  the  heat  transfer  from  the  top  surface 
becomes  dominant. 

Simplified  One-dimensional  Model  (1-D) 

In  the  following  solution  it  is  assumed  that  the  transient  temperature 
distribution  within  the  porcelain  disk  is  similar  to  that  of  a  semi-infinite  plate 
(Fig.  3.1)  subjected  to  convective  cooling.  Here  the  governing  equation  is  : 


Tt  =  KT: 


zz 


(44) 


and  is  subjected  to  the  initial  condition 
and  boundary  conditions     at  left    (z  =  0) 

at  right  (z  =  a) 


T(z,0)  =  Ti 


(45) 


Tz  (0,t)  =  0 

Tz  (a,t)  + 1  [T(a,t)  -  TJ  =  0 
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where  again  h  is  an  average  heat  transfer  coefficient,  K  is  thermal  conductivity 
and   k  is  thermal  diffusivity  of  porcelain.    To  convert  equation    (44)  in 
nondimensional  form,  let: 

to  yield  :  0X  =  6^  (47) 


with  initial  condition  (I.C.)  :  6(T|,0)  =  1  (48) 

and  boundary  conditions  (left)    :  &rj  (0,x)  =  0 

(right):  erj(l,T)  +  Bie(l,T)  =  0  , 

where  again,  Bi  =  Biot  number  =  -rr- 

Many  heat  transfer  books  have  used  the  method  of  separation  of  variables  to 
solve  (47)  as  follows: 

Let  0  (r|,T)  =  A(r|)  •  Q(x)  (49) 

After  substituting  in  (47)  and  applying  the  initial  and  boundary  conditions  (48) 
one  finds  : 


o  /•«  ^     i    S  r        t  a  2   v  sin(|3nT)cos(pnr|) 

e  (h,t)  =  2  j  {exP  (-pn  x)  pn+[sin(pnT)C0S(PnTl)]  >  (50) 


where  (3n  are  characteristics  of  the  transcendental  equation  : 


tan(pn)  =  j|  (51) 
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To  find  the  nondimensional  temperature,  0,  first  one  must  solve  (51),  using  an 
iteration  technique  such  as  the  Newton-Raphson  method  to  get  the  characteristic 
values  of  (3n.  Then  at  each  point  in  space  (i\)  and  time  (x) ,  the  temperature  (9)  is 
determined  from  (50).  Because  of  the  negative  power  of  the  exponential 
function,  this  equation  usually  has  a  rapid  rate  of  convergence.  Therefore,  the 
closed  form  solution  by  the  separation  of  variable  approach  requires  the  solution 
of  equation  (50)  and  a  numerical  solution  of  equation  (51). 

An  alternative  approach  used  in  this  study  is  the  Laplace  transformation 
method. 

The  Laplace  transformation  is  defined  by: 

L  [f(x)]=  J[f(x)e     ]dx  =  F(s)  (52) 

o 

where  f  (x)  =  L  ~l  [F  (s)]  (53) 

Taking  L  of  both  sides  of  (47)  and  applying  the  initial  condition  will  yield  an 
ordinary  differential  equation  (ODE)  : 

©T|Tj - s  0  =-1  ,  where  0  =L (9)  (54) 


Solving  this  ODE  one  has: 


G  (T|,  s)  =  A  sinh  (Vs~ T|  )  +  B  cosh  (^[s~v\ )  +  1/s  (55) 

subjected  to  the  boundary  conditions  (left)    :      0  (0,  s)  =  0  (56) 

(right) :       ©n  +  Bi@=0     at  T]  =  1 
Applying  the  boundary  condition  and  taking  the  inverse  Laplace  transform  (ZT1) 
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of  both  sides  will  yield  : 


s[Vs  sinh(-vs  )+Bicos(Vs  )] 


For  given  values  of  Bi  and  t)  ,  the  L-1  term  is  easily  found  by  using  available 
computer  programs  such  as  IMSL  subroutines.  In  equation  (57),  s  is  a  parameter 
corresponding  to  the  nondimensional  time,  x  .  The  advantage  of  this  method 
over  the  previous  separation  approach  is  that  the  temperature  (0)  can  be  obtained 
by  solving  only  one  equation.  The  higher  efficiency  of  this  method  becomes 
more  apparent  when  the  temperatures  are  calculated  for  various  values  of  the  heat 
transfer  coefficient.  Also  the  computer  programming  of  this  method  is  less 
difficult  than  for  the  separation  of  variable  approach. 

Finite  Element  Models 

A.  Axisvmmetric  Model  of  Ceramic  Disks  for  Thermal  Stress  Analysis 

As  stated  earlier,  an  axisymmetric  finite  element  model  was  developed  to 
determine  the  transient  temperature  and  stress  distributions  within  the  bilayered 
porcelain  disks  by  considering  the  heat  loss  from  the  side  of  the  disk.  In  this 
model,  it  was  assumed  that  the  bottom  surface  (side  DC  in  Fig.  3.3)  was 
completely  insulated  and  the  disk  specimen  was  initially  at  a  uniform  temperature. 
The  shaded  cross  section  of  the  disk  was  divided  into  153  nodes  and  128 
quadrilateral  elements  and  the  boundary  conditions  shown  in  the  figure  were 
applied  in  the  thermal  analysis.  The  node  which  is  labeled  103  and  is  located 
near  point  A  on  the  y  axis  corresponds  to  the  thermocouple  location  for  each 
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Fig.  3.3    Axisymmetric  thermal  model  and  boundary  conditions  for  the 
finite  element  analyses. 
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disk.  The  total  heat  transfer  coefficient ,  ht ,  calculated  from  equation  (24)  was 
used  to  impose  the  boundary  conditions  along  the  sides  AB  and  BC. 

Shown  in  Fig.  3.4  are  the  node  and  material  type  numbers  and  the  boundary 
conditions  used  in  the  finite  element  stress  analyses.  The  nodes  located  on  the 
symmetric  axis  correspond  to  zero  displacement  along  the  horizontal  axis.  A  zero 
displacement  along  the  vertical  direction  was  also  imposed  on  the  node  1.  For 
the  bilayered  disk  the  thickness  ratio  of  body  (1)  to  opaque  porcelain  (2)  was 
assumed  to  be  3  (Fig.  3.4)  which  represents  the  clinical  conditions.  Shown  in 
Fig.  3.5  are  the  element  numbers  used  in  the  finite  element  model.  The  transient 
temperature  and  stress  distribution  within  disks  subjected  to  various  cooling 
conditions  were  calculated  by  the  use  of  the  finite  element  program,  ANSYS 
(Swanson  Analysis  System,  Inc.,  Houston,  PA).  To  improve  the  accuracy  of 
the  nodal  temperatures  and  to  minimize  the  computer  time  required  for  the 
unsteady-state  thermal  analysis,  the  automatic  time  step  optimization  option  of 
ANSYS  was  employed.  Therefore,  at  the  initial  stage  of  cooling  the  time  step 
(At)  had  the  smallest  value  of  0.5  sec  which  increased  gradually  as  the  slope  of 
the  cooling  curve  reduced  at  the  later  stages.  The  calculated  nodal  temperatures 
as  a  function  of  time  were  stored  in  a  computer  file  to  be  recalled  as  input  data  for 
the  axisymmetric  thermal  stress  analysis.  A  glass  transition  (Tg)  value  of  600°C 
was  assumed  for  each  porcelain  type  and  the  calculated  temperatures  near  this  Tg 
value  were  used  to  determine  the  transient  stress  distribution.  Shown  in  Table 
3.1  are  the  thermal  and  mechanical  properties  of  porcelains  used  in  the  finite 
element  models. 
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Fig.  3.5    Element  numbers  and  boundary  conditions  used  in  the  finite  element 
model  of  a  bilayered  ceramic  disk. 
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Table  3.1  Thermal  and  Mechanical  Properties  of  Porcelain  used  in  Finite 
Element  Analysis. 


PROPERTY 

UNIT 

VALUE 

Thermal  Conductivity 

W/m-K 

1.68 

Density 

kg/m3 

2.29  x 

103 

Specific  Heat 

J/kg-K 

1.68  x 

103 

Glass  Transition 

°C 

600 

Modulus  of  Elasticity 

GPa 

70.3 

Poisson's  Ratio 

— 

0.28 

Thermal  Contraction  Mismatch 

ppm/°C 

-1.5,  0.0, 

+3.2 

B.  Three-Dimensional  Model  of  Dental  Ceramic  Crowns  for  Static  Analysis 

Stress  analyses  were  performed  for  ceramic  thicknesses  of  0.5  mm  (case  I),  1 .5 
mm  (case  II),  and  3.0  mm  (case  III)  at  the  occlusal  region  (Fig.  3.6).  The  elastic 
modulus  values  used  for  ceramic,  dentin  and  pulp  were  7.03  x  10  ,  1.86  x  10 
and  2.0  MPa,  respectively,  and  the  corresponding  Poisson's  ratios  were  0.28, 
0.31  and  0.45,  respectively.  The  analyses  were  performed  on  an  IBM  PC- AT 
microcomputer  using  the  SAP86  (Number  Cruncher  Microsystems,  Inc.,  San 
Francisco,  CA)  finite  element  program.  All  the  components  of  normal,  shear  and 
principal  stresses  and  displacements  along  the  three  directions  of  the  global 
coordinate  system  were  examined  by  use  of  stress  contour  maps. 

Shown  in  Fig.  3.7  is  the  three-dimensional  finite  element  model  of  the 
molar  crown  which  consisted  of  586  nodes  and  456  solid  elements.  All  the 
nodes  in  the  x-y  plane  which  correspond  to  the  root  portion  of  dentin,  were 
assumed  to  be  fixed.   A  distributed  load  of  600  Newtons  which  is  an  average 
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Fig.  3.6    Mesial  view  of  a  mandibular  second  molar  showing 
three  design  cases. 
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Fig.  3.7  Three-dimensional  finite  element  model  of  the  molar  crown. 
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biting  force  on  a  molar  tooth,  was  applied  at  the  occulusal  region  along  vertical 
(Fig.  3.8)  and  horizontal  directions. 

C  Axisvmmetric  Model  of  Dental  Ceramic  Crowns  for  Thermal  Stress  Analysis 

Shown  in  Fig.  3.9  is  the  finite  element  model  of  the  molar  crown  used  in 
the  stress  analysis.  Because  the  crown  is  axisymmetric  with  respect  to  a  vertical 
axis  passing  through  AB,  only  half  of  the  model  (Fig.  3.10)  consisting  of  164 
nodes  and  136  elements  was  used  in  the  analysis.  The  bilayered  ceramic  crown 
consisted  of  an  approximately  1 .5  mm  thick  layer  of  opaque  (a  =  14.2  pprnfC) 
and  a  4.5-mm-thick  layer  of  body  porcelain  (a  =  1 1 .0  ppm/°C).  In  the  thermal 
analysis,  the  inside  surface  was  assumed  to  be  completely  insulated  and  the 
outside  surface  was  subjected  to    uniform  convective  cooling.   The  transient 
temperature  distributions  were  determined  for  free  convective  cooling 
(h  =70  W/m2-K)  and  tempering  by  forced  air  (h  =560  W/m2-K).    The 
temperatures  were  calculated  for  an  initial  time  step  (At)  of  0.5  sec  and  the  time 
step  optimization  option  of  the  ANSYS  program  was  used  to  increase  the  At 
value.    The  calculated  temperatures  were  incorporated  in  the  stress  analysis 
program   to  determine   transient  stresses   for  assumed  glass   transition 
temperatures,  (Tg  ),   of  600°C  and  also  400°C.   The  transient  stresses  were 
calculated  relative  to  these  Tg  values  for  the  crowns  cooled  from  the  initial 
temperatures  of  850°C  and  750°C.  The  node  A  located  on  the  symmetric  axis 
(Fig.  3.10)  was  assumed  to  be  fixed  along  the  x  and  y  directions.  Points  C  and 
D  (Figs  3.9  and  3.10)  identify  the  cuspal  region  of  the  crown  which  has  the 
maximum  ceramic  thickness.  The  temperature  versus  time  corresponding  to  each 
cooling  treatment  was  plotted  at  these  points  using  the  postprocessing  options  of 
the  ANSYS  finite  element  program. 


74 


600  N 


Fig.  3.8    Locations  of  the  applied  loads  in  the  finite  element  model. 
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Fig.  3.9    Cross  section  of  a  ceramic  molar  crown  used  for  finite  element 
thermal  stress  analysis. 
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Fig.  3.10  Axisymmetric  finite  element  model  of  the  bilayered  ceramic  molar 
crown  consisted  of  opaque  and  body  porcelain  layers. 


CHAPTER  IV 

RESULTS 

Thermal  Analyses 

Thermal  Contraction  Curves 

The  percent  thermal  contraction  data  as  a  function  of  temperature  for  the  two 
opaque  and  two  body  porcelain  types  (Fig.  2.2)  were  obtained  from  the  cooling 
curves.  Each  curve  corresponds  to  the  mean  values  obtained  from  five 
specimens  in  the  temperature  range  between  25°C  to  500°C,  using  a  reference 
temperature  of  25  °C.  The  average  coefficient  of  linear  thermal  contraction  (a)  for 
each  porcelain  group  (five  specimens  in  each  group)  which  was  calculated  from 
the  slope  of  these  cooling  curves,  agreed  well  with  the  values  supplied  by  the 
manufacturer  (J.F.  Jelenko  and  Co.,  Armonk,  NY).  Although  the  transition  of 
dental  ceramics  from  solid  to  liquid  occurs  in  a  temperature  range  (glass  transition 
region)  which  depends  on  the  type  of  porcelain  and  thermal  parameters  such  as 
heating  and  cooling  rates  used,  in  this  study  a  constant  glass  transition 
temperature  (Tg)  of  600°C  was  assumed  for  the  finite  element  analysis. 

Temperature  Distribution  within  the  Dental  Furnace 

A  schematic  illustration  of  the  furnace  heating  chamber  and  the  steady-state 
temperature  distribution  inside  this  chamber  at  982° C  was  shown  in  Fig.  2.3.  The 
mean  temperature  values  of  the  sides  (S\  and  S2),  top  (T)  and  bottom  (B) 

thermocouples  at  each  distance  from  the  furnace  door  indicate  that  there  is  a 
relatively  larger  temperature  gradient  from  the  center  toward  the  furnace  door 
compared  to  the  back  wall  which  is  located  at  a  distance  of  1 15  mm  from  the  door. 
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Although  the  furnace  thermocouple  is  controlled  at  a  temperature  of  982  C,  this 
thermal  gradient  exists  because  of  the  heat  losses  from  the  front  and  back  walls  of 
the  chamber  and  the  nonuniform  arrangements  of  the  electric  heating  elements.  At 
each  distance  there  was  a  larger  difference  between  the  top  (T)  and  bottom  (B) 
thermocouples  compared  to  the  difference  between  Si  and  S2  (see  Fig.  2.3).  The 

average  temperature  values  at  the  side  (S)  and  the  overall  mean  (M)  corresponding 
to  the  four  thermocouples  as  a  function  of  distance  are  plotted  in  Fig.  4.1.  The 
shaded  region  which  is  located  between  the  distance  of  65  mm  and  75  mm  is 
related  to  the  maximum  temperature  value  of  each  curve  and  is  identified  as  the  hot 
zone  region  of  the  heating  chamber.  To  minimize  the  initial  temperature  gradient 
within  the  specimen,  each  disk  was  placed  inside  the  furnace  near  this  region  and 
15  sec  after  the  desired  set  point  temperature  was  reached,  it  was  rapidly  removed 
from  the  furnace  for  various  cooling  treatments. 

Experimental  Cooling  Profiles 

Shown  in  Fig.  4.2  are  the  measured  temperature  versus  time  profiles  of  a 
2-mm-thick  disk  subjected  to  free  convection  (FC),  tempering  by  horizontal  flow 
(TH)  and  tempering  by  vertical  flow  (TV)  from  an  initial  temperature  of  982°C, 
using  air  as  a  cooling  medium.  The  largest  slope  of  the  cooling  curve  corresponds 
to  the  TV  case  in  which  the  disk  temperature  at  the  thermocouple  location  (0.5  mm 
from  the  surface)  reached  the  ambient  temperature  (30°C)  after  50  sec.  However, 
because  the  bottom  surface  was  still  at  a  higher  temperature,  the  air  blasting  was 
continued  for  approximately  90  sec  until  the  entire  specimen  reached  ambient 
temperature.  The  time  required  for  the  specimen  to  reach  the  ambient  temperature 
was  90  sec  for  TH  and  8.5  minutes  for  FC  conditions.  The  glass  transition 
temperature  (To)  of  600°C  for  TV,  TH  and  FC  cases  were  reached  after  4.5,  16 
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Fig.  4.1  Temperature  versus  distance  from  the  door  of  the  test  furnace  at 
a  control  temperature  of  982°C. 
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Fig.  4.2    Temperature  versus  time  profiles  of  a  2-mm-thick  disk  for  free 
convective  cooling,  tempering  by  horizontal  and  vertical  flows. 
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and  56  sec,  respectively.  As  the  experimentally  determined  cooling  curves  for  the 
6  mm-thick  disk  (Fig.  4.3)  also  indicate,  tempering  by  vertical  air  flow  produced  a 
more  rapid  cooling  rate  than  the  horizontal  flow.  For  this  specimen  which  has  a 
diameter  of  32  mm,  the  Tg  value  was  reached  after  4.5  sec  (TV),  34  sec  (TH)  and 
72  sec  (FC).  A  comparison  among  the  cooling  profiles  of  a  2-mm  thick  disk 
subjected  to  slow  cooling  inside  the  furnace  (SC),  fast  cooling  (FC),  and 
tempering  by  vertical  flow  of  air  (TV)  are  shown  in  Fig.  4.4.  For  the  disks  that 
were  slow-cooled  from  an  initial  temperature  of  982°C,  the  ambient  temperature  of 
30°C  was  reached  after  approximately  4  hours. 

Shown  in  Fig.  4.5  is  the  comparison  between  the  capabilities  of  air  and  helium 
gases  for  cooling  of  a  2-mm-thick,  1 6-mm-diameter  ,  and  6-mm  thick,  32-mm- 
diameter  specimens  by  free  convection.  In  both  specimens,  the  slope  of  the 
cooling  curves  for  helium  gas  were  larger  than  the  values  for  air.  As  was 
expected,  because  of  the  difference  between  their  heat  capacities,  the  time  required 
for  the  6  mm-thick  disk  to  reach  ambient  temperature  was  much  greater  than  that 
for  the  2  mm-thick  disk.  The  experimentally  measured  temperature  versus  time 
profiles  of  a  2  mm-thick  disk  subjected  to  free  convection  (FC)  and  tempering  (TV) 
in  helium  and  air  are  shown  in  Fig.  4.6.  For  the  TV  case,  there  was  a  negligible 
difference  between  the  cooling  curves  obtained  from  helium  and  air.  Both  of  the 
cooling  media  led  to  an  ambient  disk  temperature  of  30  C  after  about  30  sec.  In  the 
free  convection  case  (FC),  the  helium  gas  was  more  effective  than  air  but  the 
difference  between  the  slopes  of  the  curves  near  the  glass  transition  temperature 
(Tg)  was  relatively  small.  Because  of  the  similarity  between  the  thermal 
characteristics  of  nitrogen  and  oxygen  with  air,  additional  cooling  experiments 
were  carried  out  only  for  nitrogen  gas. 
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Fig.  4.3  Temperature  versus  time  profiles  of  a  6-mm-thick  disk  for  three 


air  cooling  conditions. 
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Fig.  4.4     Temperature  versus  time  profiles  in  a  2-mm-thick  disk  using 
slow  cooling,  free  convective  cooling  and  tempering  in  air. 
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Fig.  4.5    Temperature  versus  time  profile  of  a  2-mm  (t  =  2)  and  a  6-mm 
(t  =  6)  thick  disk  cooled  by  free  convection  in  air  and  helium. 
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Fig.  4.6    Cooling  profiles  of  a  2-mm-thick  disk  subjected  to  free  convection 
and  tempering  in  helium  and  air. 
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Shown  in  Fig.  4.7  are  the  cooling  profiles  for  2  mm-thick  and  6  mm-thick 
specimens  tempered  by  vertical  (TV)  and  horizontal  (TH)  flow  of  nitrogen.  This 
figure  also  shows  that  vertical  tempering  is  more  effective  than  horizontal 
tempering  and  the  difference  is  more  pronounced  for  the  6  mm-thick  disk  which 
has  a  larger  diameter  (32  mm).  For  the  horizontal  flow  condition,  the  thickness  of 
the  thermal  boundary  layer  increases  with  the  distance  from  the  nozzle, 
consequently,  the  heat  transfer  or  the  effectiveness  of  tempering  is  reduced  along 
this  direction. 

Heat  Transfer  in  Gaseous  and  Liquid  Cooling  Media 

The  heat  transfer  coefficient  (h)  is  a  function  of  the  cooling  medium  and  the 
method  of  cooling.  For  free  convective  cooling  of  a  body  in  air,  h  is  in  the  range 
between  5  to  50  W/m2-K,  while  for  forced  convection  the  value  can  increase  up  to 
500  W/m2-K  (Fig.  4.8).  When  water  is  used  as  cooling  medium,  the  value  of  h 
can  be  as  high  as  4,000  for  free  convection  and  15,000  W/m2-K  for  forced 
convection.    Shown  in  Fig.  4.9  is  the  comparison  among  the  equivalent  heat 
transfer  coefficients  (h)  for  a  2-mm  thick  porcelain  porcelain  disk  cooled  by  free 
convection  (FC),  radiation  and  tempering  by  vertical  air  flow  (TV).  The  h  values 
for  FC  and  TV  conditions,  which  were  determined  from  equations  (16)  and  (21), 
are  approximately  independent  of  specimen  surface  temperature.   However,  the 
equivalent  heat  transfer  coefficient  for  radiation  (hr)  which  is  calculated  from 
equation  (23)  and  exists    in  both  FC  and  TV  cooling  cases  is  significantly 
influenced  by  the  surface  temperature.  Therefore,  the  total  heat  transfer  coefficient 
(ht  in  eq.  24)  is  a  temperature  dependent  property,  during  the  initial  stage  of 
cooling.     As  shown  in  Fig.  4.10,  the  major  mode  of  heat  loss  at  higher 
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temperatures,  which  for  the  cooling  condition  in  this  study  is  called  "free 
convection"  or  "fast  cooling"  (FC),  is  actually  radiative  heat  transfer.  At  the  initial 
stage  of  cooling,  the  heat  loss  by  radiation  is  approximately  7  times  greater  than  the 
convection  heat  loss.  As  the  specimen  temperature  drops  below  300°C,  the 
amount  of  heat  losses  by  convection  and  radiation  become  comparable.  The  ht 
values  in  the  temperature  range  between  30  and  982°C  which  were  obtained  by 
superposition  of  convective  and  radiative  heat  transfer  coefficients  were  determined 
for  each  cooling  condition  and  used  in  the  thermal  models  to  calculate  the  transient 
temperature  distributions  within  each  disk. 

As  stated  in  Chapter  III,  the  heat  transfer  coefficient  is  independent  of  the 
material  properties  of  the  body.  However,  the  thermal  history  of  the  body  is 
controlled  by  the  heat  transfer  coefficient.  Shown  in  Fig.  4.1 1  is  the  calculated 
temperature  versus  time  profile  determined  from  equation  6  for  a  thin  copper  disk 
for  some  arbitrarily  chosen  values  of  the  heat  transfer  coefficient. 

To  compare  the  cooling  capabilities  of  various  media  and  to  determine  an 
average  h  value  associated  with  each  cooling  medium,  a  thin  copper  disk,  1  mm  in 
thickness  and  16  mm  in  diameter,  was  employed.  Experimentally  measured 
cooling  profiles  for  free  convection  in  air,  silicone  oils  with  kinematic  viscosities 
of  50  and  1000  centistokes,  and  water  are  shown  in  Fig.  4.12.  The  largest 
temperature  gradient  was  obtained  in  the  water  quenched  specimens  in  which  the 
ambient  temperature  was  reached  after  7  sec.  The  silicone  oil  with  a  kinematic 
viscosity  of  50  centistokes  resulted  in  a  steeper  cooling  profile  because  it  is 
associated  with  a  larger  heat  transfer  coefficient  compared  with  that  of  oil  at  a  larger 
viscosity. 

To  determine  the  influence  of  porcelain  and  silicone  oil  soaking  temperatures  in 
producing  thermal  shock  failure,  several  quenching  experiments  were  conducted 
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by  reducing  the  initial  temperature  of  the  disk  and/or  increasing  the  initial 
temperature  of  the  silicone  oil.  When  porcelain  disk  specimens  were  quenched 
from  an  initial  temperature  of  982° C  in  the  silicone  oils  at  room  temperature, 
thermal  shock  failure  occurred  in  almost  all  bilayered  specimens.  But  when  the 
initial  temperature  was  850°C  and  the  oil  was  at  100°C,  thermal  shock  occurred  in 
approximately  5%  of  the  bilayered  (Act  =+3.2  or  Aa  =  -1.5)  porcelain  disks.  In 
general,  the  quenching  experiments  with  the  single  layer  disks  (Aa  =  0)  were  more 
successful  (  fewer  thermal  shock  failures)  than  that  of  the  bilayered  (Aa  =  -1.5  or 
Aa  =  +3.2)  disks,  probably  because  the  incompatibility  tensile  stresses  were  not 
produced  in  the  first  group. 

Determination  of  the  Average  Heat  Transfer  Coefficients 

DeHoff  and  Anusavice  [66]  assumed  a  constant  heat  transfer  coefficient  (h)  of 
57  W/m2-K  in  the  solution  of  a  semi-infinite  plate  (eq.  A-l),  to  determine  the 
transient  temperature  distribution  in  a  single-layer  porcelain  plate  subjected  to  free 
convective  cooling  from  an  initial  temperature  of  800"C.  Although  the  value  of  h 
was  chosen  based  on  thermocouple  measurements  reported  in  a  study  by  Twiggs  et 
al.  [101],  the  calculated  cooling  profile  was  not  compared  with  any  experimental 
data.  In  the  present  study,  the  h  value  of  57  W/m2-K  was  used  in  equations  (35) 
and  (57)  to  determine  the  temperature  versus  time  profile  in  a  porcelain  disk 
subjected  to  free  convective  cooling  (FC)  from  an  initial  temperature  of  982°C. 
The  results  from  1-D  (57)  and  2-D  (35)  thermal  models  and  the  comparison  with 
the  experimental  data  are  shown  in  Fig.  4.13.  The  heat  transfer  coefficient  of  57 
W/m2-K  applied  in  the  1-D  model,  overestimates  the  temperatures.   However, 
when  it  was  applied  in  the  2-D  solution,  better  agreement  between  the  model  and 
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the  experiment  was  obtained  near  the  glass  transition  region.  This  difference  is 
because  of  the  heat  loss  from  the  side  of  the  disk  specimen  which  is  not  taken  into 
account  in  the  1  -D  analysis. 

The  experimental  cooling  profile  for  a  copper  disk  was  used  to  determine  the 
mean  value  of  the  heat  transfer  coefficient  for  free  convective  cooling  of  a  porcelain 
disk.  The  h  values  of  60  and  80  W/m2-K  were  used  in  the  lumped  analysis  (eq.  6) 
to  calculate  the  transient  temperatures  and  the  results  were  compared  with  the 
experimental  data  (Fig.  4.14).  These  results  indicate  that  the  cooling  profile  of  a 
copper  disk  subjected  to  free  convection  can  be  predicted  by  using  a  heat  transfer 
coefficient  between  60  and  80  W/m2-K.  When  an  h  value  of  70  W/m2-K  was 
applied  in  equation  (57),  the  calculated  temperature  versus  time  profile  agreed  well 
with  that  of  the  experimental  values.  Shown  in  Fig.  4.15  is  the  comparison 
between  the  1-D  model  and  the  experimental  data  for  a  porcelain  disk  using  the 
non-dimensional  forms  of  time  and  temperature  (eq.  46).  To  determine  the  cooling 
profile  from  the  2-D  analysis  with  the  heat  transfer  coefficient  of  70  W/m2-K,  the 
characteristic  values  of  (38)  and  (39)  were  calculated  by  the  Newton-Raphson 
method.  The  approximate  values  of  characteristics  corresponding  to  the 
transcendental  equation  (39)  were  found  by  graphical  methods  (see  Fig.  4.16)  and 
were  applied  in  equation  (43).  For  the  conditions  used  in  this  analysis,  the  first  10 
positive  characteristics  of  equations  (38)  and  (39)  yielded  acceptable  accuracy  in 
the  temperatures  obtained  from  equation  (35). 

Shown  in  Fig.  4.17  is  a  comparison  between  the  results  from  the  1-D  and  2-D 
analyses  with  the  experimental  data  corresponding  to  a  thermocouple  which  was 
embedded  along  the  symmetric  axis  of  the  porcelain  disk  at  a  distance  of  0.5  mm 
from  the  surface.  This  figure  indicates  that  a  heat  transfer  coefficient  of  70 
W/m2-K  is  suitable  for  predicting  the  thermal  history  based  on  the  1-D  analysis. 
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Fig.  4.14  Experimental  and  the  calculated  cooling  profiles  for  a  thin  copper 
disk  subjected  to  free  convection  using  a  heat  transfer  coefficients 
of  60  W/m2°K  and  80  W/m2.  K. 
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a  heat  transfer  coefficient  of  70  W/m2-  K. 
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However,  the  2-D  model  here  yields  lower  temperatures  than  the  measured  values. 
The  results  in  Figs.  4.13  and  4.17  suggest  that  the  cooling  profiles  of  the  porcelain 
disk  can  be  predicted  with  variable  accuracy  depending  on  the  type  of  model  (1-D 
or  2-D)  and  the  value  of  the  heat  transfer  coefficient  used.  For  free  convective 
cooling  the  heat  loss  from  the  side  cannot  be  neglected  compared  with  the  loss 
from  the  surface.  Thus,  the  h  value  of  57  W/m2-K  used  by  DeHoff  and  Anusavice 
[66]  in  a  1-D  model  of  a  porcelain  plate  will  result  in  a  more  accurate  cooling  curve 
when  it  is  applied  in  a  2-D  model  of  a  porcelain  disk.  However,  to  make  the 
computations  less  difficult,  further  analyses  were  performed  based  on  the  results 
from  the  1-D  model  for  a  higher  heat  transfer  coefficient  value  of  70  W/m2-K  (Fig. 
4.17). 

Shown  in  Fig.  4.18  is  a  comparison  between  the  cooling  profiles  for  a 
tempered  disk  by  a  vertical  air  flow  from  the  1-D  and  the  2-D  models.  Similar  to 
the  previous  method,  an  average  value  of  h  was  obtained  using  the  slope  of  the 
cooling  profile  from  a  copper  disk.  The  mean  values  of  the  total  heat  transfer 
coefficient  used  in  the  1-D  thermal  model  for  specimens  quenched  in  air,  silicone 
oil  and  water  were  560,  720  and  1590  W/m2-K,  respectively.  As  Fig.  4.18 
indicates,  for  a  large  value  of  h  (  560  W/m2-K  or  higher),  the  results  from  the  1-D 
and  2-D  models  are  similar  above  the  glass  transition  region.  This  supports  our 
hypothesis  that,  for  such  rapid  cooling  condition,  the  heat  transfer  from  the  top 
surface  of  the  disk  is  much  greater  than  that  from  the  side  so  that  side  heat  loss  can 
be  considered  negligible.  Therefore,  the  one-dimensional  model  is  sufficient  for 
calculation  of  the  temperature  and  from  that  the  stress  distributions  within  the 
tempered  disks  in  air  and  in  silicone  oil.  The  stresses  for  these  conditions  were 
calculated  based  on  the  results  from  equation  (57). 
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Fig.  4.18    Cooling  profiles  calculated  from  the  1-D  and  2-D  models  for 
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Stress  Analyses  using  Viscoelasticitv  and  Structural  Relaxation 
A.  Gaseous  Cooling  Media 


Us- 


Computer  programs  in  Basic  language  were  developed  by  P.  H.  DeHoff 
(Mechanical  Engineering  and  Engineering  Science  Department,  University  of 
North  Carolina  at  Charlotte)  on  an  IBM  microcomputer  to  determine  the  transient 
and  residual  stresses  induced  within  a  porcelain  plate  during  a  cooling  treatment. 
These  programs  which  are  referred  to  as  "fictive  temperature"  and  "stress" 
programs,  incorporate  linear  viscoelasticity  and  structural  relaxation  behavior  of 
porcelain  and  were  used  by  DeHoff  and  Anusavice  [66]  to  calculate  the  thermal 
stresses  in  a  single  layer  porcelain  plate  symmetrically  cooled  on  both  sides  from 
an  initial  temperature  of  800°C.  A  summary  of  the  procedures  employed  in  this 
analysis  is  given  in  chapter  I  and  in  Appendix  A.  However,  for  a  more  complete 
description  of  the  analytical  and  numerical  techniques,  one  should  refer  to  the 
original  study  [66].  In  the  following  calculations  of  the  thermal  stresses,  the  edge 
effect  of  the  porcelain  disk  is  neglected  and  it  is  assumed  that  the  stress 
distributions  produced  within  a  porcelain  disk  are  similar  to  that  within  a  porcelain 
plate. 

The  temperature  versus  time  data  which  were  calculated  from  the  1-D  thermal 
model  (eq.  57)  for  each  cooling  condition  were  incorporated  in  the  "fictive 
temperature"  program  and  the  transient  fictive  temperature  distribution  for  that 
cooling  condition  was  calculated.  The  fictive  temperature  is  the  quenching 
temperature  of  glass  corresponding  to  a  particular  structural  state  (see  eq.  A-3). 
Using  procedures  utilized  by  DeHoff  and  Anusavice  [66],  a  glassy  and  a  liquid 
state  for  the  body  porcelain  were  defined  by  a  linear  regression  fit  to  the  measured 
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thermal  contraction  data.  To  compare  the  results  of  this  study  with  the  results  of 
analyses  performed  in  the  original  paper  [66],  the  contraction  data  at  temperatures 
greater  than  400°C  were  used  for  the  liquid  and  at  temperatures  below  400°C  for 
the  glassy  states  (elastic  solid  state)  of  porcelain.  The  transient  fictive  temperatures 
calculated  for  each  disk  and  the  thermal  contraction  values  for  the  liquid  and  glassy 
states  were  applied  in  the  "stress"  program  to  determine  the  induced  transient  and 
residual  stresses.  It  should  be  emphasized  that  in  this  study,  the  thermal 
contraction  coefficients  were  measured  at  a  constant  cooling  rate  of  3°/min. 
However,  this  property  can  be  significantly  influenced  by  the  rate  of  cooling 
employed  in  the  dilatometer. 

Shown  in  Fig.  4.19  is  the  actual  and  fictive  transient  temperatures  at  the 
surface  and  the  insulated  face  of  a  2-mm-thick  disk  subjected  to  free  convective 
(FC)  cooling  in  air  from  an  initial  temperature  of  982°C.  The  fictive  temperatures 
reached  an  asymptotic  value  of  approximately  600°C  after  70  sec.  Above  this 
temperature,  there  is  only  small  difference  between  the  actual  and  the  fictive 
temperatures.  Also,  at  70  sec  after  cooling,  there  is  a  negligible  temperature 
gradient  within  the  disk.  In  the  "stress"  program,  the  actual  and  fictive 
temperature  distribution  within  the  disk  at  any  given  time  are  fit  by  quadratic 
equations.  The  corresponding  transient  stresses  at  the  top  and  bottom  (insulated) 
surfaces  as  a  function  of  time  are  shown  in  Fig.  4.20.  At  approximately  300  sec 
after  cooling  when  the  disk  reached  a  constant  temperature  (Fig.  4.19),  the 
stresses  also  approached  a  fixed  value  (Fig.  4.20).  The  magnitude  of  the 
compressive  stresses  within  the  top  surface  and  tensile  stresses  at  the  bottom 
surface  increased  rapidly  at  about  70  sec  after  initial  cooling.  The  residual  stresses 
within  the  disk  versus  the  distance  from  the  surface  are  shown  in  Fig.  4.21.  The 
maximum  compressive  and  tensile  stress  values   of  24.2  and  11.1  MPa  which 
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Fig.  4.19    Actual  and  fictive  temperature  versus  time  profiles  at  the  surface 
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4.20    Transient  stresses  versus  time  at  the  surface  and  the  insulated  face 
of  a  2-mm-thick  body  porcelain  disk  cooled  by  free  convection 
from  an  initial  temperature  of  982°C. 
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Fig.  4.21    Residual  stress  distribution  in  a  2-mm- thick  body  porcelain  disk 
subjected  to  free  convective  cooling  from  an  initial  temperature 
of 982°C. 
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correspond  to  a  ratio  of  2.1,  were  produced    at  the  top  and  bottom  surfaces, 
respectively. 

Figs.  4.22  to  4.24  correspond  to  the  temperature  and  stress  data  for  free 
convective  cooling  of  a  6-mm-thick  disk.  Compared  with  the  disk  with  a  thickness 
of  2  mm,  larger  actual  and  fictive  temperature  gradients  are  produced  within  this 
disk  at  any  given  time  (Fig.  4.22).  A  fictive  temperature  of  570°C  was  reached  at 
the  surface  after  180  sec.  During  the  initial  stage  of  cooling,  relatively  large  tensile 
stresses  were  induced  at  the  surface  (Fig.  4.23).   At  about  30  sec  after   initial 
cooling,  these  stresses  rapidly  decreased  in  magnitude  and  became  compressive  in 
nature.    However,  as  stated  in  Chapter  I,  these  tensile  stresses  are  highly 
undesirable.  Experimental  cooling  of  several  6-mm  thick  disks  resulted  in  thermal 
shock  failure,  specifically  for  the  bilayered  (Aoc  =+3.2  or  -1.5)  specimens.  The 
failure  occurred  during  various  stages  of  free  convective  cooling  in  which  pieces  of 
porcelain  were  violently  separated  from  the  disk.  The  cause  of  the  thermal  shock 
failure  was  the  development  of  relatively  large  tensile  stresses  at  the  top  surface 
during  the  initial  stage  and  at  the  interior  during  a  later  stage  of  cooling.  These 
stresses  were  developed  because  of  the  combined  effect  of  the  thermal  contraction 
incompatibility   and  temperature  gradient  within  the  6-mm-thick  disks.    The 
residual  stress  distribution  as  a  function  of  distance  from  the  top  surface  is  shown 
in  Fig.  4.24.   This  figure  shows  that  if  a  disk  does  not  break  because  of  initial 
tensile  stresses  at  the  surface,  it  will  have  a  compressive  stress  layer  approximately 
2.42  mm  in  depth.  Since  the  preparation  and  cooling  treatments  of  a  6-mm-thick 
disks  were  associated  with  several  failures,  additional  experiments  were  performed 
only  for  the  2-mm  thick  specimens. 

To  determine  the  influence  of  the  initial  temperature  on  the  effectiveness  of 
tempering,  the  single  layer  body  porcelain  disks,  (Aoc  =  0.0),  having  a  2  mm 
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Fig.  4.22    Actual  and  fictive  temperature  versus  time  profiles  at  the  surface. 
and  the  insulated  face  of  a  6-mm-thick  body  porcelain  disk 
cooled  by  free  convection  from  an  initial  temperature  of  982"C. 
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Fig.  4.23    Transient  stresses  versus  time  at  the  surface  and  the  insulated 
face  of  a  6-mm-thick  body  porcelain  disk  cooled  by  free 
convection  from  an  initial  temperature  of  982°C. 
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thickness,  were  tempered  by  a  vertical  air  jet  (TV)  from  initial  temperatures  of  982, 
850,  750  and  650°C.  Shown  in  Figs  4.25  to  4.37  are  the  results  of  thermal  stress 
analyses  for  these  tempered  specimens.  Six  disks  in  each  group  were  prepared  and 
the  effectiveness  of  the  tempering  treatment  was  analyzed  experimentally  by  the 
biaxial  flexure  test.  These  results  will  be  presented  in  a  separate  section.  As  stated 
earlier,  the  average  heat  transfer  coefficient  of  560  W/m2-K,  was  determined  based 
on  the  slope  of  the  cooling  curve  from  a  copper  disk.  Shown  in  Fig.  4.25  are  the 
actual  and  fictive  temperature  versus  time  profiles  at  the  surface  and  the  insulated 
face  of  a  2-mm-thick  body  porcelain  disk  tempered  from  an  initial  temperature  of 
982°C.  Compared  with  the  profiles  from  the  free  convective  cooling  (Fig.  4.19), 
larger  temperature  gradients  for  actual  and  fictive  curves  were  produced  within  the 
tempered  specimens  independent  of  the  initial  tempering  temperature  (Figs.  4.25, 
4.28,  4.31  and  4.34).  The  asymptotic  values  obtained  for  fictive  temperatures 
were  not  significantly  influenced  by  the  initial  tempering  temperature.  In  all 
cooling  profiles,  the  tempered  disks  reached  an  ambient  temperature  of  30°C  in 
about  60  seconds.  This  is  consistent  with  the  experimentally  determined  cooling 
profiles  presented  earlier. 

Shown  in  Fig.  4.26  are  the  transient  stresses  versus  time  at  the  surface  and 
insulated  face  of  a  2-mm-thick  body  porcelain  disk  tempered  from  an  initial 
temperature  of  982°C.  Approximately  10  sec  after  initiation  of  tempering,  the 
transient  compressive  stresses  produced  at  the  surface  increased  in  magnitude  and 
in  about  45  sec  reached  a  constant  value.  This  was  also  true  for  the  disks  tempered 
from  lower  initial  temperatures  (Figs.  4.29,  4.32  and  4.35).  However,  the 
magnitude  of  the  transient  stresses  at  any  given  time  was  influenced  by  the  initial 
soak  temperature.  The  transient  stress  profiles  indicate  that  an  increase  in  the 
surface  compression  corresponds  to  a  greater  amount  of  undesirable  interior  tensile 


UJ 


en 


Q_ 


20  30  40 

TIME  (SEC) 


Fig.  4.25    Actual  and  fictive  temperature  versus  time  profiles  at  the  surface 
and  the  insulated  face  of  a  2-mm-thick  body  porcelain  disk 
tempered  from  an  initial  temperature  of  982°C. 
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Fig.  4.26  Transient  stresses  versus  time  at  the  surface  and  the  insulated 
face  of  a  2-mm-thick  body  porcelain  disk  tempered  from  an 
initial  temperature  of  982°C. 
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Fig.  4.27     Residual  stress  distribution  in  a  2-mm-thick  body  porcelain  disk 
tempered  from  an  initial  temperature  of  982°C. 
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Fig.  4.28    Actual  and  fictive  temperature  versus  time  profiles  at  the  surface 
and  the  insulated  face  of  a  2-mm-thick  body  porcelain  disk- 
tempered  from  an  initial  temperature  of  850°C. 
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Fig.  4.29  Transient  stresses  versus  time  at  the  surface  and  the  insulated 
face  of  a  2-mm-thick  body  porcelain  disk  tempered  from  an 
initial  temperature  of  850°C. 
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Fig.  4.30     Residual  stress  distribution  in  a  2-mm-thick  body  porcelain  disk 
tempered  from  an  initial  temperature  of  850°C. 
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Fig.    4.31    Actual  and  fictive  temperature  versus  time  profiles  at  the 

surface  and  the  insulated  face  of  a  2-mrn-thick  body  porcelain 
disk  tempered  from  an  initial  temperature  of  750°C. 
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Fig.  4.32     Transient  stresses  versus  time  at  the  surface  and  the  insulated 
face  of  a  2-mm-thick  body  porcelain  disk  tempered  from  an 
initial  temperature  of  750°C. 
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Fig.    4.33    Residual  stress  distribution  in  a  2-mm-thick  body  porcelain  disk 
tempered  from  an  initial  temperature  of  750"C. 
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Fig.  4.34    Actual  and  fictive  temperature  versus  time  profiles  at  the  surface 
and  the  insulated  face  of  a  2-mm-thick  body  porcelain  disk 
tempered  from  an  initial  temperature  of  650°C. 
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Fig.  4.35    Transient  stresses  versus  time  at  the  surface  and  the  insulated  face 
of  a  2-mm-thick  body  porcelain  disk  tempered  from  an  initial 
temperature  of  650°C. 
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Fig.  4.36    Residual  stress  distribution  in  a  2-mm-thick  body  porcelain  disk 
tempered  from  an  initial  temperature  of  650°C. 
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Fig.  4.37    Residual  compressive  and  tensile  stresses  at  the  surface  of  a 

2-mm-thick  body  porcelain  disk  as  a  function  of  initial  tempering 
temperature. 
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stresses.  However,  in  general  the  surface  compressive  stresses  produced  in  a  disk 
at  any  given  time  are  greater  in  magnitude  than  the  interior  tension.  The  residual 
stress  distribution  in  a  2-mm-thick  body  porcelain  disk  tempered  from  an  initial 
temperature  of  982°C  is  shown  in  Fig.  4.27.  The  maximum  compressive  and 
tensile  stresses  produced  at  the  surface  and  the  insulated  face  were  171 .4  MPa  and 
62.6  MPa,  respectively.  Under  this  condition,  the  magnitude  of  the  surface 
compression  was  2.7  time  greater  than  that  of  the  maximum  tension.  The 
beneficial  compressive  stresses  induced  at  the  surface  of  the  tempered  disks  from 
initial  temperature  of  850,  750  and  650°C,  were  165.8  MPa  (Fig.  4.30),  155.9 
MPa  (Fig.  4.33)  and  102.0  MPa,  respectively.  These  stresses  were  greater  than 
their  corresponding  maximum  tensile  stresses  by  factors  of  2.3,  2.6  and  1.7,  for 
the  initial  tempering  temperatures  of  850,  750,  and  650°C,  respectively.  As  stated 
earlier,  the  corresponding  factor  for  the  free  convective  cooling  was  2.1.  This 
indicates  that,  the  ratio  of  the  surface  compression  to  the  maximum  tension  is  not 
strongly  influenced  by  the  initial  tempering  temperature.  However,  a  small  change 
in  this  ratio  can  make  a  significant  difference  in  the  strength  and  fracture 
characteristics  of  porcelain.  Gardon  [24],  introduced  the  idea  of  "modulated 
quenching"  to  control  this  ratio  and  to  obtain  the  desired  strength  and  fracture 
characteristics  in  a  glass  plate. 

Shown  in  Fig.  4.37  is  the  residual  compressive  stresses  produced  at  the 
surface  of  a  2-mm-thick  porcelain  disk  as  a  function  of  initial  tempering 
temperature.  An  increase  of  the  soak  temperature  from  650  to  750°C  will  result  in 
a  relatively  large  increase  in  the  compressive  stress  value.  However,  any  further 
increase  in  the  initial  temperature  does  not  significantly  increase  the  surface 
compression.  Also,  a  greater  increase  in  the  ratio  of  the  surface  compression  to  the 
maximum  tension  corresponds  to  an  increase  in  the  initial  tempering  temperature 
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from  650  to  750°C.  Since  the  structural  distortion  of  the  specimens  due  to  viscous 
flow  may  occur  at  high  temperatures,  for  the  same  tempering  effectiveness  a 
smaller  soak  temperature  is  more  desirable.  In  this  study  an  optimum  value  for  the 
soak  temperature  is  determined  from  this  result  and  the  data  obtained  from  the 
biaxial  flexure  test. 

B.  Liquid  Cooling  Media 

Shown  in  Fig.  4.38  to  4.40  are  the  results  for  the  2-mm-thick  porcelain  disk 
quenched  in  silicone  oil  from  an  initial  temperature  of  850°C.  The  mean  value  of 
heat  transfer  coefficient  for  this  condition  was  720  W/m2-K  which  was  applied 
both  on  the  top  and  bottom  surfaces.  Because  of  the  symmetric  cooling  condition, 
the  midplane  of  the  disk  was  treated  like  an  insulated  boundary  and  the  transient 
temperature  distributions  were  obtained  from  the  results  of  the  Laplace 
transformation  approach  (equations  46  and  57)  using  an  "a"  dimension  of  1  mm. 
At  the  initial  stage  of  quenching,  a  relatively  large  temperature  gradient  from  the 
surface  to  the  midplane  was  produced  (Fig.  4.38)  because  of  the  large  value  of  the 
heat  transfer  coefficient.  The  fictive  and  the  actual  temperatures  approached  their 
asymptotic  values  after  5  and  25  sec,  respectively.  Since  the  entire  disk  reached 
the  oil  temperature  of  100°C  in  about  25  sec  (Fig.  4.38),  the  transient  stresses  at 
the  surface  and  midplane  subsequently  stayed  constant  (Fig.  4.39).  Shown  in 
Fig.  4.40  is  the  symmetric  distribution  of  the  residual  stresses  within  the  disk  as  a 
function  of  the  distance  from  the  center.  A  maximum  compressive  stress  value  of 
107.8  MPa  was  produced  at  the  surfaces.  The  largest  residual  tensile  stress  value 
which  occurred  at  the  midplane  of  the  disk  was  44.3  MPa.  Under  this  condition, 
the  ratio  of  the  surface  compression  to  the  maximum  interior  tension  was  2.4  and 
the  thickness  of  the  compressive  layer  at  each  side  was  approximately  0.4  mm. 
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Fig.  4.38    Actual  and  fictive  temperature  versus  time  profiles  at  the  surface 
and  midplane  of  a  2-mm-thick  body  porcelain  disk  quenched  in 
silicone  oil  from  an  initial  temperature  of  850°C. 
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Fig.  4.39     Transient  stresses  versus  time  at  the  surface  and  midplanc  of 

a  2-mm-thick  body  porcelain  disk  quenched  in  silicone  oil  from 
an  initial  temperature  of  850°C. 
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Fig.  4.40    Residual  stress  distribution  in  a  2-mm-thick  body  porcelain  disk 
quenched  in  silicone  oil  from  an  initial  temperature  of  850°C. 
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For  the  same  initial  temperature  of  850°C,  the  maximum  stresses  produced  in  the 
oil  quenched  specimens  were  smaller  than  those  of  the  tempered  specimen  in  air. 
However,  the  results  of  a  2-mm-thick  disk  tempered  from  one  side  and  insulated 
from  another  side  are  the  same  as  a  4-mm-thick  disk  symmetrically  tempered  from 
both  sides.  Therefore,  because  of  a  larger  effective  thickness,  at  a  given  time,  the 
temperature  gradient  produced  within  the  disks  during  n  on -symmetrical  cooling  in 
air  was  greater  than  the  corresponding  values  during  symmetrical  cooling  in 
silicone  oil.  It  should  be  emphasized  that  the  silicone  oil  temperature  was  kept  at 
100°C  to  prevent  thermal  shock  failure.  To  analyze  the  effect  of  a  cooling  medium 
with  a  very  large  heat  transfer  coefficient  on  the  temperature  and  stress 
characteristics  of  a  single  layer  porcelain,  the  results  of  quenching  a  2-mm-thick 
disk  in  water  were  determined  and  are  shown  in  Figs.  4.41  and  4.42.  The  initial 
temperatures  of  the  disk  and  water  were  850°C  and  100°C,  respectively.  The  mean 
value  of  the  heat  transfer  coefficient  applied  in  this  analysis  was  1590  W/m2-K 
which  was  determined  from  the  slope  of  the  cooling  curve  for  the  copper  disk  (Fig. 
4.12).  Approximately,  15  sec  after  quenching,  the  entire  disk  reached  the  water 
temperature  of  100°C  (Fig.  4.41).  Midplane  tension  and  the  surface  compression 
reached  constant  values  of  91.0  MPa  and  243.0  MPa,  respectively  (Figs.  4.42  and 
4.43).  Comparison  between  Figs.  4.40  and  4.43  indicates  that  under  this 
condition,  the  depth  of  the  residual  compressive  stress  layer  is  not  influenced  by 
the  heat  transfer  coefficient  and  is  approximately  0.4  mm  for  both  cases. 
However,  under  this  condition,  thermal  shock  failure  occurred  because  of  the  large 
tensile  stresses  which  developed  within  the  surface  at  the  initial  stage  or  at  the 
interior  during  later  stages  of  quenching. 
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Fig.    4.41     Actual  and  fictive  temperature  versus  time  profiles  at  the  surface 
and  midplane  of  a  2-mm-thick  body  porcelain  disk  quenched  in 
water  from  an  initial  temperature  of  85CTC. 
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Fig    4  42    Transient  stresses  versus  time  at  the  surface  and  midplane  of 
a  2-mm-thick  body  porcelain  disk  quenched  in  water  from 
an  initial  temperature  of  850°C. 
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Fig.  4.43    Residual  stress  distribution  in  a  2-mm-thick  body  porcelain  disk 
ouenched  in  water  from  an  initial  temperature  of  850"C. 
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Stress;  Analyses  using  Finite  Element  Method 
A  Transient  Stress  Analysis  of  Bilavered  Ceramic  Disks 

The  stress  analysis  by  viscoelastic  and  structural  relaxation  theory  which  was 
presented  earlier  only  determined  the  stresses  for  a  single-layered  ceramic  disk  by 
using  a  specific  stress  condition  of  a  two-dimensional  model.  This  model 
assumed  that  the  stress  components  along  the  horizontal  and  vertical  directions 
were  identical.  In  this  section,  the  influence  of  a  thermal  contraction  mismatch  on 
the  transient  stresses  within  a  bilayered  disk  was  analyzed  by  using  an 
axisymmetric  model.  However,  this  model  assumed  that  the  viscoelastic  stresses 
were  relaxed  during  cooling  through  the  glass  transition  region. 

Shown  in  Fig.  4.44  is  a  comparison  between  the  experimental  and  the 
calculated  cooling  profiles  for  free  convective  cooling  of  a  2-mm-thick  disk  from 
an  initial  temperature  of  982°C  in  helium.  At  the  time  that  the  temperature  at  the 
thermocouple  location  (node  number  103  in  Figs.  3.3  and  3.4)  approached  the 
glass  transition  temperature  (Tg  =  600°C),  the  difference  between  the  experimental 
and  the  calculated  temperatures  was  negligibly  small.  Shown  in  Fig.  4.45  are  the 
temperature  contours  near  this  glass  transition  temperature  which  are  perpendicular 
to  the  insulated  surface  (bottom)  and  the  symmetrical  (left  side)  boundaries.  The 
largest  temperature  value  (599.8°C)  corresponds  to  node  1  and  the  regions  adjacent 
to  it  (see  Fig.  3.4  for  the  identification  of  node  numbers)  and  the  smallest  value 
(496. 9°C)  corresponds  to  the  upper  right  corner  node  (number  153).    This 
temperature  distribution  was  obtained  for  a  time  of  41  sec  after  the  initiation  of 
cooling.   The  vector  plot  of  the  thermal  fluxes  are  shown  in  Fig.  4.46.   These 
vectors,  which  are  perpendicular  to  the  thermal  contours,  determine  the  direction 
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Fig.  4.  44      Experimental  and  calculated  cooling  profiles  for  a  2-mm-thick 
disk  subjected  to  free  convective  cooling  in  helium. 
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Fig.  4.45     Temperature  distribution  near  glass  transition  temperature  within 
a  2-mm-thick  disk  subjected  to  free  convective  cooling  in  helium. 


A 


f 
f 
/ 


i 
? 


I 

/ 

/ 

/ 


7-  /!■  /•  7:  A./L/L 


/ 

/ 


/ 

/ 


A  /■  /■  S  z1 

/■    /*  /*  s*   /* 

/*•    y*  /*  y    /* 

S*  S*  >*      ^* 


/t 

sr 


Fig.  4.46    Vector  plot  of  thermal  fluxes  for  free  corrective  cooling 
of  a  2-mm-thick  disk  in  helium. 
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and  magnitude  of  the  heat  loss  at  each  location.  The  nodes  which  are  located  on 
the  symmetrical  axis  and  adjacent  to  it  are  not  significantly  affected  by  the  heat  loss 
from  the  side  of  the  disk.  The  temperature  gradients  at  these  regions  are  parallel  to 
the  vertical  axis. 

Shown  in  Figs.  4.47  and  4.48  are  principal  stress  contours,  Gl  and  G2,  within 
a  bilayered  porcelain  disk  with  a  mismatch  level  of  +3.2  ppm/°C.  The  peak  value 
of  the  transient  tensile  stress  (<J1)  was  78.6  MPa  (contour  E)  which  occurred  near 
the  edge  of  the  disk.  The  peak  transient  compressive  stress  value  was  68.2  MPa 
which  occurred  at  the  top  surface  and  near  the  insulated  boundary.  These  stresses 
developed  because  of  the  combined  effect  of  the  temperature  gradient  within  the 
disk  and  the  thermal  contraction  difference  between  the  porcelain  layers. 

Shown  in  Fig.  4.49  are  the  temperature  contours  near  Tg  for  free  convective 
cooling  of  a  2-mm-thick  disk  from  an  initial  temperature  of  850°C  using  air  as 
cooling  medium.  In  this  study,  similar  values  of  thermal  conductivity  (K),  specific 
heat  (Cp)  and  density  (p)  were  used  for  all  porcelain  types.  Therefore,  the 
temperature  distributions  for  single-layered  and  bilayered  disks  are  identical. 
However,  because  of  the  difference  in  thermal  contraction  mismatch,  the  transient 
stress  distribution  for  each  porcelain  combination  was  different,  shown  in  Figs. 
4.50  and  4.51  are  the  principal  stress  contours,  Gl  and  G2,  for  free  convective 
cooling  of  a  bilayered  disk  with  a  mismatch  level  of  +3.2  ppm/°C.  Comparison  of 
these  figures  with  Figs.  4.47  and  4.48  indicates  that  the  difference  among  the 
stress  distributions  for  disks  subjected  to  cooling  treatments  in  helium  and  in  air  is 
negligibly  small. 

Shown  in  Figs.  4.52  and  4.53  are  the  contours,  CI  and  G2,  for  a  bilayered 
porcelain  with  a  mismatch  level  of  -1.5  ppm/°C  subjected  to  free  convective 
cooling  in  air.  The  peak  values  of  principal  tensile  stress  (Gl)  and  compressive 
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Fig.  4.50  Transient  principal  stress  (Gl  in  MPa)  produced  within  a 
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Fig.  4.51  Transient  principal  stress  (o2  in  MPa)  produced  within  a 
2-mm-thick  disk  (Aoc  =  +3.2  ppm/'C)  subjected  to  free 
convective  cooling  in  air  (Tj  =  850°C). 
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(02)  stress  were  55.1  MPa  and  49.3  MPa,  respectively.  Although  these  stresses 
were  smaller  than  the  corresponding  values  in  a  +3.2  ppmfC  mismatch  case, 
comparison  of  Figs.  4.50  and  4.52  indicates  that  under  this  condition  the  tensile 
stresses  produced  at  the  top  surface  are  relatively  larger.  Shown  in  Fig.  4.54  are 
the  superimposed  plots  of  the  deformed  and  undeformed  conditions  for  a  disk  with 
a  mismatch  level  of  Aa  =  +3.2  ppm/°C.    The  maximum  displacement  along  the 
vertical  direction  was  0.02  mm  which  occurred  at  point  B  (node  153  in  Fig.  3.4). 
Since  the  incompatibility  stresses  are  not  present  in  a  single-layered  specimen  (Aa 
=  0.0),  the  disk  is  not  subjected  to  bending.  Under  this  condition,  the  maximum 
vertical  displacement  was  0.012  mm  (Fig.  4.55).   Shown  in  Figs.  4.56  and  4.57 
are  stress  contours,  Ol  and  02  for  the  free  convective  cooling  of  a  2-mm  thick 
body  porcelain  disk.  Comparison  of  these  figures  with  Fig.  4.49  indicates  that  in 
the  absence  of  incompatibility  stresses,  the  stress  contours  are  in  the  same  direction 
as  the  temperature  contours  and  are  perpendicular  to  the  insulated  and  symmetrical 
boundaries.     Comparison  of  Figs.  4.56  and  4.57  with  Figs.  4.50  and  4.51 
indicates  that  the  peak  Ol  and  02  values  for  a  mismatch  level  of  +3.2  ppm/°C  were 
at  least  two  times  larger  than  the  corresponding  values  for  the  body  porcelain 
disks.  Therefore,  under  the  free  convective  cooling  condition,  the  influence  of  the 
thermal  contraction  mismatch  on  the  transient  stresses  was  significant.  Shown  in 
Figs.  4.58  and  4.59  are  the  Ol  and  02  stress  distributions  at  the  surface  of  a  body 
porcelain  as  a  function  of  disk  radius.   The  minimum  and  maximum  transient 
tensile  stresses  (oi)  occurred  at  the  center  and  edge  of  the  disk,  respectively  (Fig. 
4.58).  However,  the  maximum  compressive  stress  (02)  at  the  surface  occurred  at  a 
point  with  radius  of  5.0  mm  (Fig.  4.59).  The  principal  stresses  at  the  edge  of  the 
disk  as  a  function  of  distance  from  the  bottom  surface  are  shown  in  Figs.  4.60  and 
4.61.  The  transient  tensile  stress  distribution  (Ol)  along  this  direction  (Fig.  4.60) 


Fig.  4.54    Superimposed  plots  of  the  deformed  and  undeformed  disk 

(Aa  =  +3.2  ppmfC)  under  the  influence  of  transient  thermal 
and  incompatibility  stresses. 
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Fig.  4.58    Transient  principal  stress  (Gl  in  MPa)  produced  at  the  surface  of 
a  2-mm-thick  disk  (Aa  =  0)  subjected  to  free  convective  cooling 
in  air  as  a  function  of  radius  (Ti  =  850°C). 
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Fig.  4.59    Transient  principal  stress  (C2  in  MPa)  produced  at  the  surface  of 
a  2-mm-thick  disk  (Act  =  0)  subjected  to  free  convective  cooling 
in  air  as  a  function  of  radius  (Ti  =  850°C). 
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is  parabolic  in  shape  with  a  minimum  value  of  31  MPa  at  the  bottom  (node  17  in 
Fig.  3.5)  and  maximum  value  of  44  MPa  at  the  top  surface  (node  153).  The  peak 
value  of  G2  stress  was  13.0  MPa  (Fig.  4.61)  which  occurred  at  the  bottom  edge. 

Shown  in  Fig.  4.62  are  the  temperature  contours  near  Tg  (600°C)  within  a  6- 
mm-thick  disk  subjected  to  free  convective  cooling  in  air  from  an  initial  temperature 
of  982°C.    These  temperature  occurred  146  sec  after  the  initiation  of  cooling. 
Compared  with  the  2-mm-thick  disk  (Fig.  4.49),  this  disk  was  subjected  to  a  larger 
temperature  gradient  along  the  thickness  and  along  the  radius.   Shown  in  Figs. 
4.63  and  4.64  are  stress  contours,  the  Gl  and  02  for  the  bilayered  disk  with  a 
mismatch  level  of  3.2  ppm/°C.  The  gradient  of  principal  tensile  stress  (Fig.  4.63) 
increased  with  the  radius  and  it  was  maximum  at  the  edge  of  the  disk.  However,  at 
this  location  the  gradient  of  the  principal  compressive  stress  (02)  was  minimum. 
Although  the  transient  stresses  at  this  stage  of  cooling  are  not  significantly  greater 
than  the  corresponding  stresses  in  a  2-mm-thick  disk,  the  interaction  of  these 
stresses  with  the  microcracks  and  flaws,  which  are  probably  more  likely  to  be 
present  in  a  6-mm-thick  disk  due  to  laboratory  techniques,  can  cause  thermal  shock 
failure.  As  mentioned  earlier,  several  experiments,  which  were  conducted  for  free 
convective  cooling  of  the  6-mm-thick  disks,  resulted  in  thermal  shock  failure  at 
various  stages  of  cooling.  Therefore,  tempering  analysis  was  performed  only  for 
the  2-mm-thick  disks. 

Shown  in  Fig.  4.65  is  the  temperature  versus  time  profiles  for  three  selected 
nodes  located  on  axis  of  symmetry  in  a  2-mm-thick  disk  subjected  to  tempering  by 
air  from  an  initial  temperature  of  850°C.  At  the  initial  stage  of  tempering,  a 
significant  temperature  gradient  occurred  along  the  disk  thickness  (between  nodes 
1  and  137).  This  gradient  reduced  in  magnitude  at  a  later  time.  The  temperatures 
at  node  103  which  correspond  to  the  thermocouple  location  agreed  well  with  the 
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Fig.  4.61    Transient  principal  stress  (a2  in  MPa)  produced  at  the  edge  of 
a  2-mm-thick  disk  (Aa  =  0)  subjected  to  free  convective  coolim 
in  air  as  a  function  of  thickness  (Ti  =  850°C). 
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Fig.  4.62    Temperature  distribution  near  Tg  (600°C)  for  free  convective 
cooling  of  a  6-mm-thick  disk  in  air  (Tj  =  982°C). 
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Fig.  4.63    Transient  principal  stress  (al  in  MPa)  produced  within  a 

6-mm-thick  disk  (Aoc  =  3.2  ppm/°C)  subjected  to  free  convective 
cooling  in  air  (T{  =  982°C). 


CO 


A 

— 80.457 

B 

—58.605 

C 

— 36.754 

D 

--14.902 

E 

-6.949 

F 

-28.801 

Fig.  4.64    Transient  principal  stress  (G2  in  MPa)  produced  within  a 

6-mm-thick  disk  (Act  =  3.2  ppm/°C)  subjected  to  free  convective 
cooling  in  air  (Ti  =  982°C). 


160 


o 

o 
H 


TIME  (SEC) 


Fig.  4.65    Temperature  versus  time  profiles  for  three  selected  nodes  on  the 
symmetrical  axis  in  a  2-mm-thick  disk  subjected  to  tempering 
in  air  (Ti  =  850°C). 
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experimental  data  and  the  analytical  models  (eqs.  57  and  35).  The  axisymmetric 
temperature  distribution  which  occurred  at  8  sec  after  the  initiation  of  tempering,  is 
shown  in  Fig.  4.66.  A  maximum  temperature  differential  of  1 13°C,  occurred  near 
the  center  of  the  specimen.  At  this  stage  of  tempering,  a  temperature  differential  of 
302°C  was  produced  between  the  hottest  and  the  coldest  point  of  the  disk.  Shown 
in  Figs.  4.67-4.69  are  the  Gl  and  G2  values  and  the  stress  intensity  contours  (SI), 
respectively,  for  a  tempered  disk  with  a  contraction  mismatch  level  of  3.2  ppm/°C. 
Stress  intensity  is  defined  as  twice  the  maximum  shear  stress  (Tmax)-  The  peak 
values  of  tensile  stresses   (Fig.  4.67)  and  stress  intensity  (Fig.  4.69)  were  152.9 
MPa  and  163.8  MPa,  respectively,  which  occurred  at  the  edge  of  the  disk  (contour 
F).    The  peak  value  of  compressive  stress  at  the  surface  (Fig.  4.68)  was  88.0 
MPa.    Comparison  of  these  stresses  with  the  corresponding  values  for  the 
mismatch  level  of -1.5  ppm/°C  (Fig.  4.70  through  4.72)  indicates  that  the  contours 
are  different  in  shape  and  direction,  however,  theix  peak  values  are  not  significantly 
different.  For  instance,  the  largest  Gl  value  for  the  +3.2  ppm/°C  mismatch  level 
was  153  MPa  (Fig.  4.67)  and  the  corresponding  value  for  the  -1.5  ppm/°C  level 
was  152  MPa  (Fig.  4.70).    This  indicates  that  in  the  tempered  specimen  the 
transient  incompatibility  stresses  are  dominated  by  the  stresses  developed  because 
of  the  temperature  gradient. 

To  determine  the  contribution  of  the  temperature  gradient  to  transient  stresses, 
the  principal  stress  contours  were  plotted  (Figs.  4.73^.76)  for  a  tempered  body 
porcelain  disk  with  no  contraction  mismatch  (Aa  =  0.0).  Similar  to  the  free 
convective  cooling  condition,  the  insulated  boundary  for  the  thermal  analysis  was 
treated  as  a  symmetrical  boundary  for  stress  analysis.  The  peak  value  of  the 
transient  tensile  stress  (Gl)  was  142.4  MPa  (Fig.  4.73)  which  occurred  at  the 
surface  near  the  edge  of  the  disk.  This  location  was  also  associated  with  the  largest 
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Fig.  4.66    Temperature  distribution  near  Tg  (600°C)  in  a  2-mm-thick  disk 
subjected  to  tempering  in  air  (Ti  =  850°C). 
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Fig.  4.67    Transient  principal  stress  (ol  in  MPa)  produced  within  a 

2-mm-thick  disk  (Aa  =  3.2  ppm/"C)  subjected  to  tempering 
in  air  (Ti  =  850°C). 
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Fig.  4.68  Transient  principal  stress  (o2  in  MPa)  produced  within  a 

2-mm-thick  disk  (Aot  =  3.2  ppm/°C)  subjected  to  tempering 
in  air  (Ti  =  850°C). 
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Fig.  4.69    Transient  stress  intensity  (SI  in  MPa)  produced  within  a 

2-mm-thick  disk  (Aa  =  3.2  ppm/°C)  subjected  to  tempering 
in  air  (Ti  =  850°C). 
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Fig.  4.70    Transient  principal  stress  (Gl  in  MPa)  produced  within  a 

2-mm-thick  disk  (Aa  =  -1.5  ppm/°C)  subjected  to  tempering 
in  air  (Tj  =  850°C). 
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Fig.  4.71   Transient  principal  stress  (G2  in  MPa)  produced  within  a 

2-mm-thick  disk  (Aa  =  -1 .5  ppm/°C)  subjected  to  tempering 
in  air  (Tj  =  850°C). 
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IYansient  principal  stress  (ol  in  MPa)  produced  within  a  2-mm- 
thick  disk  (Aa  =  0)  subjected  to  tempering  in  air  (Tj  =  850°C). 
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Fig.  4.74    Transient  principal  stress  (G2  in  MPa)  produced  within  a  2-tnm- 
thick  disk  (Aa  =  0)  subjected  to  tempering  in  air  (Tj  =  850°C). 
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Fig.  4.75    Transient  principal  stress  (o3  in  MPa)  produced  within  a  2-mni- 
thick  disk  (Aa  =  0)  subjected  to  tempering  in  air  (Tj  =  850*C 
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thick  disk  (Act  =  0)  subjected  to  tempering  in  air  (Tj  =  850*C). 
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stress  intensity  value  of  141.7  MPa  (Fig.  4.76).  The  stress  contours  were  along 
the  same  directions  as  the  temperature  contours  (Fig.  4.66)  The  locations  of 
maximum  02  and  03  gradients  (Figs.  4.74,  4.75)  were  associated  with  the  regions 
of  the  largest  temperature  gradient.    Shown  in  Fig.  4.77  and  4.78  are  the 
distribution  of  01  and  02  at  the  surface  as  a  function  of  the  disk  radius.  Similar  to 
the  free  convection  case,  the  01  stress  has  a  parabolic  shape  with  a  minimum  value 
of  72  MPa  at  the  center  and  maximum  value  of  155  MPa  at  the  edge  of  the  disk 
(Fig.  4.77).  The  largest  02  value  was  84  MPa  which  occurred  at  a  radius  of  4.4 
mm.  The  regions  with  a  larger  radius  are  influenced  by  the  heat  loss  from  the  side 
and  are  associated  with  a  smaller  temperature  gradient.  Therefore,  the  principal 
compressive  stresses  (02)  at  the  surface  for  free  convective  cooling  (Fig.  4.59) 
and  tempering  (Fig.  4.78),  which  is  more  influenced  by  temperature  gradient, 
decreases  toward  the  edge  of  the  disk.    Shown  in  Figs.  4.79  and  4.80  are  the 
principal  stresses  (01  and  02)  at  the  edge  of  the  disk  as  a  function  of  distance  from 
the  insulated  (bottom)  surface.    These  stresses  are  also  parabolic  in  shape. 
However,  at  the  surface  where  01  has  a  maximum  of  156  MPa  (Fig.  4.79),  02  has 
a  minimum  of  8  MPa  (Fig.  4.80).  The  minimum  and  maximum  values  of  01  and 
02  which  occurred  at  the  bottom  surface  were  1 16  MPa  and  71  MPa,  respectively. 
Shown  in  Fig.  4.81  are  the  temperature  contours  near  the  Tg  value  of  600°C 
for  a  2-mm-thick  body  porcelain  disk  (Aoc  =  0)  quenched  from  an  initial 
temperature  of  850°C  in  silicone  oil  at  a  temperature  of  100°C.  Since  the  disk  was 
subjected  to  uniform  cooling,  the  temperature  distribution  is  symmetrical  along  a 
horizontal  plane  which  passes  through  the  half-thickness  of  the  disk.  This  plane  is 
referred  to  as  midplane  of  the  disk.   Shown  in  Fig.  4.82  is  a  vector  plot  of  the 
thermal  fluxes  which  represent  the  magnitude  and  direction  of  the  heat  loss  at  each 
location.   Shown  in  Figs.  4.83^.85  are  the  principal  stresses  (01  and  02)  and 
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g.  4.78    Transient  principal  stress  (o2  in  MPa)  at  the  surface  of  a  2-mm- 
thick  disk  (Aa  =  0)  subjected  to  tempering  in  air  (Ti  =  850° C) 
as  a  function  of  radius. 
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Fig.  4.79    Transient  principal  stress  (<Xl  in  MPa)  at  the  edge  of  a  2-mm- 
thick  disk  (Aoc  =  0)  subjected  to  tempering  in  air  (Ti  =  850°C) 
as  a  function  of  thickness. 
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Fig.  4.80  Transient  principal  stress  (o2  in  MPa)  at  the  edge  of  a  2-... 

thick  disk  (Aoc  =  0)  subjected  to  tempering  in  air  (Ti  =  850°C) 
as  a  function  of  thickness. 
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Fig.  4.8 1     Temperature  distribution  near  Tg  (600°C)  for  a  2-mm-thick  body 
porcelain  disk  (Act  =  0)  quenched  from  an  initial  temperature 
of  850°C  in  silicone  oil  at  a  temperature  of  100°C. 
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Fig.  4.82    Vector  plot  of  (hernial  fluxes  for  a  2-mm-thick  body  porcelain 
disk  (Aa  =  0)  quenched  from  an  initial  temperature  of  850°C 
in  silicone  oil  at  a  temperature  of  I(X)°C. 
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Fig.  4.83    Transient  principal  stress  (Gl  in  MPa)  within  a  2-inm-thick  disk 

(Aot  =  0)  quenched  from  an  initial  temperature  of  850°C  in  silicone 
oil  at  a  temperature  of  100°C. 
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Fig.  4.84    Transient  principal  stress  (o2  in  MPa)  within  a  2-mm-thick  disk 

(Aa  =  0)  quenched  from  an  initial  temperature  of  850V  in  silicone 
oil  at  a  temperature  of  100  V. 


A 

— 42.569 

B 

--24 .092 

C 

--5. 614 

D 

=12.863 

E 

=  31  .34 

F 

-49.818 

<p 


R 

R 

R 

R 

R 

n 

R 

R 

R             1 

A 

A 

A 

A 

A 

A 

A 

A 

A             j 

A 

A 

A 

A 

A 

A 

A 

A 

A                I 

_fl_ 


.B_ 


A 

A     , 

A 

...  A_. 

A 

A     , 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A       - 

R 

R     . 

R 

R 

R... 

R__ 

...    R 

R     . 

.R    . 

n 

c 

D 
E 
F 


=16.864 
=39.036 
=  61  .209 
=83.381 
=105.554 
=127. 727 


Fig.  4.85    Transient  stress  intensity  (SI  in  MPa)  within  a  2-mm-thick  disk 

(Aa  =  0)  quenched  from  an  initial  temperature  of  850°C  in  silicone 
oil  at  a  temperature  of  100°C. 
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stress  intensity  (SI)  contours  for  this  cooling  condition.  These  stresses  are 
parabolic  in  shape  and  similar  to  the  the  temperature  contours,  are  symmetrical 
along  the  midplane  of  the  disk.  The  peak  values  of  ai  and  SI  stresses  were  130.8 
MPa  (Fig.  4.83)  and  127.7  MPa  (Fig.  4.85),  respectively,  which  occurred  at  the 
edge  of  the  disk.  The  maximum  02  value  was  49.8  MPa  (Fig.  4.84)  which  was 
tensile  in  nature  and  occurred  at  the  surface.  Shown  in  Figs.  4.86  through  4.88 
are  the  transient  stresses  produced  at  the  surface  as  a  function  of  radius.  The 
magnitude  of  stresses,  oi  (Fig.  4.86)  and  SI  (Fig.  4.88)  increased  with  an  increase 
in  radius  and  the  largest  values  occurred  at  the  edge  of  the  disk.  However,  the 
largest  value  of  C2  occurred  at  a  radius  of  6.2  mm.  Shown  in  Figs.  4.89  through 
4.91  are  the  transient  stresses  developed  at  the  edge  of  the  disk  as  a  function  of 
distance  from  the  bottom  surface.  These  stresses  which  are  produced  because  of 
the  temperature  gradients  in  the  absence  of  incompatibility  stresses,  are  parabolic  in 
shape  and  are  symmetrical  with  respect  to  the  midplane  of  the  disk. 

The  minimum  values  of  stresses,  Gl  and  SI  were  103.0  MPa  (Fig.  4.89)  and 
110.0  MPa  (Fig.  4.91),  respectively.  The  maximum  value  of  G2  stress  was  39.0 
MPa  which  occurred  at  the  midplane  of  the  disk. 

B.  Static  Analysis  of  Dental  Ceramic  Crowns 

The  maximum  tensile  stress  contours  for  a  ceramic  thickness  of  0.5  mm 
(case  I  design)  due  to  vertically  and  horizontally  applied  loads  are  shown  in  Figs. 
4.92  and  4.93,  respectively.  The  largest  principal  tensile  stresses  induced  in 
ceramic  due  to  a  distributed  load  of  600  N  applied  in  a  cuspal  region  were  1 1 .4 
MPa  and  171.9  MPa  for  vertical  and  horizontal  loading  orientations,  respectively. 
Stresses  which  developed  in  the  facial  and  lingual  margins  were  primarily 
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Fig.  4.86    Transient  principal  stress  «7l  in  MPa)  at  the  surface  a  2-mm-thick 
disk  (Aoc  =  0)  quenched  in  silicone  oil  as  a  function  of  radius. 
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Fig.  4.S7    Transient  principal  stress  (o2  in  MPa)  at  the  surface  a  2-mm-thick 
disk  (Ace  =  0)  quenched  in  silicone  oil  as  a  function  of  radius. 
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Fig.  4.88    Transient  stress  intensity  (SI  in  MPa)  at  the  surface  a  2-mm-thick 
disk  (Aa  =  0)  quenched  in  silicone  oil  as  a  function  of  radius. 
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Fig.  4.89    Transient  principal  stress  (<Jl  in  MPa)  at  the  edge  of  a  2-mm-thick 
disk  (Aa  =  0)  quenched  in  silicone  oil  as  a  function  of  thickness. 
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Fig.  4.90    Transient  principal  stress  (d2  in  N4Pa)  at  the  edge  of  a  2-mm-thick 
disk  (Act  =  0)  quenched  in  silicone  oil  as  a  function  of  thickness. 
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Fig.  4.91    Transient  stress  intensity  (SI  in  MPa)  at  the  edge  of  a  2-mm-thick 
disk  (Aoc  =  0)  quenched  in  silicone  oil  as  a  function  of  thickness. 
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Fig.    4.92    Principal  stress  contours  (Gl  in  MPa)  in  a  ceramic  dental  crown 
due  to  vertically  applied  loads. 
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Fig.  4.93    Principal  stress  contours  (oi  in  MPa)  in  a  ceramic  dental  crown 
due  to  horizontally  onented  loads. 
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compressive  under  vertical  loads.  However,  tensile  stresses  developed  when  the 
load  was  applied  horizontally.  The  stress  distributions  within  the  crowns  three 
occlusal  thicknesses  were  different  only  near  the  site  of  loading.  Shown  in  Figs. 
4.94  and  4.95  are  the  principal  compressive  and  maximum  shear  stress  contours 
due  to  a  horizontally  oriented  load  of  600  N.  The  result  of  this  study  showed  that 
the  variation  of  occlusal  thickness  of  ceramic  between  0.5  mm  and  3.0  mm  has 
only  a  minima]  effect  on  the  induced  tensile  stress  state.  The  three-dimensional 
model  also  showed  that  the  dominant  factor  in  development  of  tensile  stress  at  the 
surface  of  dental  ceramic  is  the  orientation  of  applied  forces.  The  influence  of 
these  stresses  in  crack  development  and  growth  can  be  reduced  by  thermal 
tempering. 

C.  Transient  Stress  Analysis  of  Dental  Ceramic  Crowns 

Shown  in  Fig.  4.96  is  the  vector  plot  of  thermal  fluxes  for  free  convective 
cooling  (h=70  W/m2-K)  of  the  ceramic  crown  from  an  initial  temperature  of 
850°C.  The  crown  was  assumed  to  be  completely  insulated  from  the  inner 
surface.  The  temperature  versus  time  profiles  for  the  inner  and  outer  surfaces 
near  the  cuspal  regions  are  shown  in  Fig.  4.97.  The  thermal  contours  near  an 
assumed  glass  transition  temperature  of  600°C  are  shown  in  Fig.  4.98.  This 
temperature  distribution  occurred  12.5  sec  after  initiation  of  the  free  convective 
cooling.  Although  the  crown  is  initially  subjected  to  uniform  cooling,  because  of 
the  variation  in  thickness  and  surface  area  there  is  a  temperature  gradient  from 
cuspal  to  marginal  regions.  When  the  temperature  at  the  cuspal  region  was 
approximately  600°C  (contour  I),  the  temperature  at  the  margin  was  543.3 °C 
(contour  A).    Shown  in  Figs.  4.99  and  4.100  are  the  corresponding  principal 
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Fig.  4.94    Principal  stress  contours  (o2  in  MPa)  in  a  ceramic  dental  crown 
due  to  horizontally  oriented  loads. 
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Fig.  4.95    Maximum  shear  stress  contours  (Tmax  in  MPa)  in  a  ceramic  denta 
crown  due  to  horizontally  oriented  loads. 
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Fig.  4.96  Vector  plot  of  thermal  fluxes  for  a  bilavered  dental  ceramic 
crown  subjected  to  free  convective  coolin". 
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Fig.  4.97  Temperature  versus  time  profiles  at  the  inner  and  outer  surfaces 
of  the  bi layered  ceramic  molar  crown  (Aoc  =  +  3.2  ppm/°C) 
subjected  to  free  cooling  from  an  initial  temperature  of  850°C. 
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Fig.  4.99    Transient  stress  contours  (ai  in  MPu)  in  a  bi layered  ceramic 

molar  crown  (Act  =  +  3.2  ppm/°C)  subjected  to  free  convective 
cooling  from  an  initial  temperature  of  850°C. 
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Fig.  4.100    Transient  stress  contours  (a2  in  MPa)  in  a  biiayered  ceramic 

molar  crown  (Aa  =  +  3.2  ppm/'C)  subjected  to  free  convective 
cooling  from  an  initial  temperature  of  850°C. 
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tensile  (<Jl)  and  compressive  (g2)  transient  stresses  induced  within  a  ceramic 
crown  with  a  thermal  contraction  mismatch  of  3.2  ppm/°C.  The  maximum  tensile 
and  compressive  were  121.3  MPa  and  163.8  MPa,  respectively,  at  the  margins 
and  near  the  occlusal  region  respectively.  Fig.  4.101  shows  the  contours 
corresponding  to  a  horizontal  displacement.  The  maximum  displacement  (0.023 
mm)  occurred  at  the  margin  which  had  the  largest  distance  from  the  fixed  node 
located  at  the  symmetric  axis  (node  A  in  Fig.  3.9). 

The  temperature  versus  time  profiles  at  the  inner  and  outer  surfaces  of  a 
ceramic  crown  near  the  cuspal  regions  are  shown  in  Fig.  4.102.  The  crown  was 
subjected  to  thermal  tempering  (h  =  560  W/m2-K)  from  an  initial  temperature  of 
850° C.   The  temperature  of  entire  crown  reached  the  ambient  temperature  of 
30°C,  approximately  after  20  sec.    Shown  in  Fig.  4.103  is  the  temperature 
contours  near  the  glass  transition  temperature  of  600°C.  At  5  sec  after  initiation 
of  cooling,  the  difference  between  the  temperature  at  the  cuspal  and  marginal 
regions  was  approximately  238°C.  The  transient  principal  stresses  induced  in  the 
crown  because  of  the  thermal  contraction  mismatch  of  +3.2  ppm/°C  and  due  to 
the  temperature  gradients  produced  by  tempering  treatment  are  shown  in  Figs. 
4.104,  4.105  and  4.106.  The  maximum  principal  tensile  stress  (Gl)  value  was 
178.1  MPa  (Fig.  4.104)  which  occurred  at  the  marginal  region.  The  maximum 
values  of  G2  and  G3  were  183.2  MPa  and  218.8  MPa,  respectively,  which  were 
compressive  in  nature  and  occurred  near  the  occlusal  regions.   Shown  in  Fig. 
4.107  are  the  contours  corresponding  to  the  displacements  along  the  horizontal 
direction. 

To  determine  the  influence  of  the  initial  tempering  temperature  on  the 
magnitude  of  principal  stresses,  the  transient  stresses  were  also  calculated  for  an 
initial  tempering  temperature  of  750°C.  Shown  in  Fig.  4.108  is  the  temperature 
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Fig.  4.101    Displacement  contours  (mm)  along  horizontal  direction  in  a 
bilayered  ceramic  molar  crown  (Aa  =  +  3.2  ppmf  C)  subjected 
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Fig.  4.102    Temperature  versus  time  profiles  for  the  inner  and  outer  surfaces 
of  the  bilayered  ceramic  molar  crown  (Aa  =  +  3.2  ppm/'C) 
subjected  to  tempering  from  an  initial  temperature  of  85CXC. 
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Fig.  4.104    Transient  stress  contours  (oi  in  MPa)  in  a  bilavered  ceramic 
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Fig.  4.108    Temperature  contours  (°C)  near  the  glass  transition  temperature 
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gradient  near  the  glass  transition  temperature  which  occurred  2  sec  after  initiation 
of  tempering.  Since  the  temperatures  were  nearly  similar  to  the  values 
corresponding  to  an  initial  temperature  of  850°C  (Fig.  4.104),  the  stress 
distributions  for  these  two  cases  were  also  similar.  Shown  in  Figs.  4.109 
through  4.1 12  are  the  stress  and  displacement  contours  for  an  initial  tempering 
temperature  of  750°C.  Comparison  of  these  figures  with  the  corresponding 
figures  for  an  initial  temperature  of  850°C  indicates  that  the  reduction  of  initial 
tempering  temperature  from  850°C  to  750°C  does  not  cause  a  significant 
difference  in  the  magnitude  of  stresses. 

These  results  show  that,  for  a  ceramic  crown  with  a  relatively  large  thermal 
contraction  differential  of  3.2  ppm/°C,  the  thermal  tempering  treatment  produces 
relatively  large  transient  tensile  stresses.  However,  transient  tensile  stresses 
above  the  strength  of  dental  ceramic  were  also  developed  under  free  convective 
cooling  (Figs.  4.99  and  4.100).  Therefore,  for  the  conditions  assumed  in  this 
analysis,  failure  due  to  incompatibility  stresses  will  probably  occur  even  for  the 
conventional  cooling  of  a  crown  with  this  level  of  thermal  contraction  mismatch. 
It  is  important  to  emphasize  that  the  transient  stresses  were  calculated  for  an 
assumed  glass  transition  temperature  of  600°C.  However,  the  transition  of  dental 
porcelain  from  a  liquid  to  a  solid  state  occurs  over  a  temperature  range  and  not  at 
a  fixed  temperature.  So,  to  determine  the  influence  of  glass  transition 
temperature  on  the  magnitude  of  stresses,  the  transient  stresses  were  also 
calculated  for  an  assumed  glass  transition  temperature  (Tg)  of  400°C.  This  is  the 
same  Tg  value  which  was  used  in  the  study  by  DeHoff  and  Anusavice  [66]. 
Shown  in  Fig.  4.113  are  the  temperature  distributions  near  the  glass  transition 
temperature  of  400°C  within  a  crown  tempered  from  an  initial  temperature  of 
850°C.  This  temperature  distribution  corresponds  to  a  time  of  4.5  sec  after  the 
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Fig.  4.110    Transient  stress  contours  (c:  in  MPa)  in  a  bilayered  ceramic 
molar  crown  (Aa  =  +  3.2  ppm/°C)  tempered  from  an 
initial  temperature  of  750°C. 
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Fig.  4.1 12    Displacement  contours  (mm)  along  horizontal  direction  in  a 
bi layered  ceramic  crown  molar  crown  (Aa  =  +  3.2  ppm/°C) 
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initial  tempering.  Under  this  assumption,  a  smaller  temperature  gradient  is 
produced  within  the  ceramic  which  results  in  smaller  transient  stresses.  Shown 
in  Figs.  4.114  through  4.116  are  the  principal  stresses  induced  within  ceramic 
with  the  thermal  contraction  mismatch  level  of  +3.2  ppm/°C.  The  maximum 
tensile  stress  at  the  margin  was  reduced  to  105.2  MPa  (Fig.4.114).  The 
maximum  values  of  02  and  G3  stresses  were  125.4  MPa  (Fig.  4.1 15)  and  149.5 
MPa  (Fig.  4.116),  respectively. 

To  determine  the  contribution  of  the  temperature  gradient  produced  because 
of  tempering  treatment  to  transient  stresses,  the  principal  stress  (Fig.  4.117 
through  4.119)  and  stress  intensity  (Fig.  4.120)  contours  were  plotted  for  a 
tempered  crown  with  no  contraction  mismatch.  For  an  assumed  glass  transition 
temperature  of  600°C,  the  peak  values  of  Gl  and  G2  and  03  were  86.6  MPa  (Fig. 
4.117),  43.5  MPa  (Fig.  4.118)  and  47.2  MPa  (Fig.  4.119),  respectively.  The 
peak  value  of  the  stress  intensity  (SI)  was  71.2  MPa  (Fig.  4.120).  When  the 
temperature  distributions  near  an  assumed  glass  transition  temperature  of  400"C 
was  used  (Fig.  4.113),  because  of  smaller  temperature  gradient  the  transient 
stresses  were  reduced  significantly.  Shown  in  Figs.  4.121  through  4.124  are  the 
stress  distributions  within  a  single-layered  ceramic  crown  for  this  condition.  The 
peak  value  of  cl,  a2,  o3  and  SI  stresses  were  48.5  MPa  (Fig.  4.121),  29.6  MPa 
(Fig.  4.122),  34.6  MPa  (Fig.4.123)  and  40.1  MPa  (Fig.  4.124),  respectively. 
These  results  show  the  glass  transition  temperature  has  a  significant  influence  on 
the  magnitude  of  transient  stresses  produced  within  a  ceramic  crown  during 
tempering  treatment. 
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Fig.  4.1  14    Transient  stress  contours  (Gl  in  MPa)  in  a  bilayered  ceramic  molar 
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Fig.  4.1 15    Transient  stress  contours  (02  in  MPa)  in  a  bilayered  ceramic  molar 
crown(Aa  "  +  3-2  Ppm/°C)  tempered  from  an  initial  temperature 
of  850  C  assuming  a  glass  transition  temperature  of  400°C. 
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Fig.  4.1 16    Transient  stress  contours  (03  in  MPa)  in  a  bilayered  ceramic  molar 
crown  (Aa  =  +  3.2  ppm/°C)  tempered  from  an  initial  temperature 
of  850°C  assuming  a  glass  transition  temperature  of  4()0°C. 
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Fig.  4.1 1 7    Transient  stress  contours  (01  in  MPa)  for  a  tempered  ceramic 
crown  (Ace  =  0)  from  an  initial  temperature  of  850°C  usin" 
a  glass  transition  temperature  of  600°C. 
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Fig.  4.1 18    Transient  stress  contours  ((72  in  MPa)  for  a  tempered  ceramic 
crown  (Aa  =  0)  from  an  initial  temperature  of  850*C  usin* 
a  glass  transition  temperature  of  600°C.  '     B 
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Fig.  4.120    Transient  stress  intensity  (SI  in  MPa)  lor  a  tempered  ceramic 
crown  (Aoc  =  0)  from  an  initial  temperature  of  850'C 
using  a  glass  transition  temperature  of  6008C. 
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Fig.  4.121    Transient  stress  contours  (Gi  in  MPa)  for  a  tempered  ceramic 
crown  (Aa  =  0)  from  an  initial  temperature  of  850'C 
using  a  glass  transition  temperature  of  400°C. 
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Transient  stress  contours  (<J2  in  MPa)  for  a  tempered  ceramic 
crown  (Aoc  =  0)  from  an  initial  temperature  of  850°C 
using  a  glass  transition  temperature  of  400°C. 
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Fig.  4.123    Transient  stress  contours  (CJ3  in  MPa)  for  a  tempered  ceramic 
crown  (Aoc  =  0)  from  an  initial  temperature  of  850°C 
using  a  glass  transition  temperature  of  400°C. 
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using  a  giass  transition  temperature  of  400°C. 
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Analysis  of  Crack  Size  Induced  by  a  Microhardness  Indenter 
A.  Statistical  Analysis 

Three-factor  ANOVA  revealed  that  the  influence  of  cooling  rate,  thermal 
contraction  mismatch,  and  initial  quenching  temperature  and  their  interaction  effects 
on  induced  crack  size  (c)  was  highly  significant  (p  =  0.001).  Throughout  this 
paper,  the  dimension  'c'  (Fig.  3.8),  refers  to  one  half  of  the  measured  crack 
length,  or  the  'crack  size'.  Shown  in  Table  4.1  are  the  results  of  the  Tukey's 
multiple  range  test  analysis  for  the  comparison  of  the  mean  values  of  measured 
crack  size  in  specimens  quenched  in  silicone  oil  at  various  kinematic  viscosities. 
As  stated  earlier,  these  oils  had  the  same  values  of  thermal  conductivities,  specific 
heat  and  density.  To  avoid  thermal  shock  failure,  the  oil  quenched  experiments 
were  performed  in  silicone  fluids  at  100°C.  When  the  single-layered  specimens 
(Aa  =  0.0  ppm/°C)  with  an  initial  temperature  of  850°C  were  quenched  in  silicone 
oil,  the  differences  among  the  crack  sizes  at  different  fluid  viscosities  were  not 
significant  (p  >  0.05).  There  was  no  significant  difference  in  crack  sizes  at 
viscosities  of  50  and  5,000  centistokes  for  mismatch  levels  of  0.0  and  -1 .5 
ppm/°C  mismatch  levels  (Table  4.1).  Shown  in  Table  4.2  are  the  mean  crack  sizes 
for  disks  quenched  from  an  initial  temperature  of  750°C.  These  values  for  disks 
with  mismatch  levels  of  +3.2  ppm/°C  and  —1.5  ppm/°C  quenched  in  silicone  oils  at 
viscosities  of  50,  1000,  and  5,000  centistokes  were  not  significantly  different. 

In  Table  4.3,  the  mean  crack  size  for  specimens  quenched  in  silicone  oil 
with  a  viscosity  of  50  centistokes  at  each  initial  cooling  temperature,  was  compared 
for  the  three  contraction  mismatch  cases.  At  the  initial  temperature  of  750°C,  there 
were  no  significant  difference  among  the  induced  crack  sizes  of  the  porcelain 
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Table  4. 1        Mean  crack  size  (pun)  for  bilayered  disks  tempered 
in  silicone  oil  from  an  initial  temperature  of  850°C. 


MISMATCH 


SILICONE  OIL  VISCOSITY  (CS) 


(ppm/°C)  5000  1000  50 


+3.2  33.8  33.4  27.2 


0.0  12J_ 28^ 264 

-!-5  _4212         NS        34^ 38^ 


NS  =  Not  Significant 


Table  4.2        Mean  crack  size  (pm)  for  bilayered  disks  tempered 
in  silicone  oil  from  an  initial  temperature  of  750°C. 


SILICONE  OIL  VISCOSITY  (CS) 
MISMATCH 

(ppm/°C)  5000  1000  50 

+3-2  312 37JS 3L8 

0-0  293 3L5  37.3 

-1-5  36.6  38.8  33.7 
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Table  4.3        Mean  crack  size  (pm)  for  bilayered  disks  tempered  in 
silicone  oil  at  kinematic  viscosity  of  50  centistokes. 


INITIAL 
TEMPERATURE 

CC) 


+3.2 


MISMATCH  (ppm/°C) 


0.0 


-1.5 


850 
750 
650 


27.2 


28.2 


38.2 


31.8 

37.3 

33.7 

37.4 

NS 

29.2 

32.5 

t 

♦ 

NS  =  Not  Significant 


Table  4.4         Mean  crack  size  (pm)  and  95%  confidence  intervals  for 
bilayered  disks  tempered  in  silicone  oil  and  in  air. 


MISMATCH 
(ppmTC) 


INITIAL  TEMPERING  TEMPERATURE  (°C) 


650 


750 


850 


SILICONE  OIL  (50  CS) 


+3.2 

37.4  ±  2.4 

31.8  ±2.1 

27.2  ±  2.3 

0.0 

29.2  ±  2.0 

37.3  ±2.9 

28.2  ±  1.4 

-1.5 

32.5  ±2.0 

33.7  ±  2.4 

38.2  ±  4.8 

982 


AIR 


44.9  ±  2.7 
38.3  ±  1.4 
38.7  ±  2.6 
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combinations.  For  the  contraction  mismatch  level  of  +3.2  ppm/°C,  the  crack  size 
decreased  from  37.4  |im  to  27.2  (im  when  the  initial  tempering  temperature  was 
increased  from  650°C  to  850°C.  Shown  in  Table  4.4  are  the  mean  crack  size  and 
95%  confidence  intervals  for  specimens  quenched  in  silicone  oil  from  initial 
temperatures  of  650,  750,  and  850°C.  These  values  are  compared  with  the 
corresponding  values  of  tempered  specimens  by  air  blasting  from  an  initial 
temperature  of  982°C.  For  the  contraction  mismatch  levels  of -1.5  ppm/°C,  the 
mean  crack  size  of  specimens  quenched  in  oil  from  850°C  were  not  significantly 
different  than  those  of  specimens  which  were  tempered  in  air.  It  should  be 
emphasized  that  the  2-mm-thick  specimens  quenched  in  silicone  oil  were  cooled 
symmetrically  from  both  sides,  while  under  the  air  cooling  condition,  one  surface 
was  insulated.  As  previously  stated,  under  the  same  cooling  rate,  the  temperature 
distribution  in  a  2-mm- thick  disk  cooled  on  one  side  and  perfectly  insulated  from 
another  side  is  the  same  as  that  in  a  4-mm-thick  disk  symmetrically  cooled  on  both 
sides.  However,  quenching  experiments  in  silicone  oil  were  conducted  using  2- 
mm-thick  disks  rather  than  4-mm-thick  disk  to  make  the  experiments  more 
clinically  relevant  and  to  reduce  the  risk  of  thermal  shock  failure.  Small  values  of 
95%  confidence  intervals  (Table  4.4)  indicate  that  the  surface  of  body  porcelain  in 
each  disk  is  relatively  homogeneous  and  that  20  indentations  per  specimen  are 
sufficient  for  statistical  analysis  of  tempering  treatment. 

Shown  in  Table  4.5  are  the  mean  values  of  crack  sizes  for  slow-cooled, 
fast-cooled  (cooled  by  free  convection)  and  tempered  specimens  in  air  from  an 
initial  temperature  of  982°C.  For  all  contraction  mismatch  levels,  the  crack  sizes 
corresponding  to  tempered  specimens  were  significantly  smaller  than  the  values  of 
the  slow-cooled  (SC)  and  fast-cooled  (FC)  specimens.  The  mean  values  of  crack 
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Table  4.5        Mean  crack  size  (|im)  for  bilayered  disks  cooled 
in  air  from  an  initial  temperature  of  982°C. 


MISMATCH 


COOLING  CONDITION 


(ppm/°C)  SC  FC 


+3.2  51.6  51.8  44.9 


51.6 

51.8 

57.5 

51.0 

46.0 

45.8 

0.0  57.5  51.0  38.3 


-1.5  46.0  45.8  38.7 


Table  4.6     Surface  stresses  in  body  porcelain  disks  cooled  by  free  convection 
(FC)  and  tempering  (T)  in  air,  and  quenched  in  silicone  oil  (OQ). 


MISMATCH 

SURFACE  STRESS  (MPa) 

(ppm/°C) 

FC 

T 

OQ 

+3.2 

-4.2 

-34.7 

-231.9 

0.0 

-11.9 

-84.1 

-218.7 

-1.5 

-29.6 

-81.0 

-84.6 
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size  for  the  SC  and  FC  conditions  were  not  significantly  different  for  all  three 
mismatch  conditions. 

B.  Calculation  of  Surface  Compressive  Stresses 

The  calculated  surface  stresses  in  body  porcelain  for  fast-cooled,  tempered 
and  oil-quenched  specimens  are  summarized  in  Table  4.6.  Relatively  small 
differences  in  crack  size  for  specimens  cooled  in  air  and  silicone  oil  (Table  4.4) 
resulted  in  large  differences  in  calculated  surface  stresses  for  the  corresponding 
mismatch  levels.  For  a  mismatch  level  of  +3.2  ppm/°C,  the  ratio  of  the  mean  crack 
size  for  the  specimens  tempered  in  air  (44.9  (im)  to  the  corresponding  value  of 
those  tempered  in  silicone  oil  (27.2  |im)  is  1.6,  while  the  compressive  surface 
stress  in  oil  (231.9  MPa)  is  6.7  times  greater  than  the  corresponding  value  in  air 
(34.7  MPa).  This  is  due  to  the  high  sensitivity  of  equations  (9)  and  (11)  to  the 
input  values  of  mean  crack  size.  Although  the  surface  compressive  stresses  in 
disks  which  were  quenched  in  oil  were  greater  than  that  specimens  tempered  in  air, 
quenching  in  silicone  oil  from  high  temperatures  were  associated  with  a  higher  risk 
of  thermal  shock  failure  compared  with  air  blasting. 

Evaluation  of  the  Biaxial  Flexure  Test  Data 

Tables  4.7  and  4.8  compare  the  biaxial  flexure  strength  for  specimens 
cooled  in  air  from  various  initial  temperatures.  The  increase  of  initial  cooling 
temperatures  (Tj)  from  650  to  850°C  did  not  have  a  major  effect  on  the  flexural 
strength  of  SC  and  FC  specimens.  However,  for  specimens  tempered  in  air,  the 
increase  of  Tj  from  650  to  850°C  resulted  in  a  25%  increase  in  flexural  strength 


233 


Table  4.7        Mean  flexural  strength  (MPa)  of  body  porcelain  cooled  in  air 
and  quenched  in  silicone  oil  (OQ)  from  various  initial 
temperatures. 


COOLING 


INITIAL  TEMPERATURE  (°C) 


TREATMENT  650  750  850 


SC  66.3  65.2  63.5 


FC  62.8  58.0  53.2 


T  80.2  100.0  107.6 


OQ*  42.1  59.3  83.1 


*  Oil  quenched  specimens  with  20  indentations  on  the  surface  of  body  porcelain. 


Table  4.8        Mean  flexural  strength  (MPa)  and  95%  confidence  intervals  of 
body  porcelain  disks  cooled  in  air  by  various  cooling  methods. 


INITIAL 


COOLING  CONDITION 


TEMPERATURE  (°C)  SC  FC 


982  46.4  +  5.4  38.9  ±  5.0  116.3  ±  14.8 

850  63.5  ±4.9  53.2+8.6  107.6±11.5 


46.4  ±  5.4 

38.9  ±5.0 

63.5  ±  4.9 

53.2  ±8.6 

65.2  ±  4.2 

58.0  ±  4.6 

750  65.2  ±4.2  58.0  ±  4.6  100.0  ±2.1 


650  66.3  ±  7.0  62.8  ±  7.2  80.2  ±  6.2 
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(Table  4.7).  There  was  no  significant  difference  (p  >  0.05)  between  the  flexural 
strength  of  specimens  tempered  in  air  at  750  and  850°C.  The  flexural  strength  of 
the  oil  quenched  specimens  are  also  given  in  Table  4.7.  It  should  be  emphasized 
that  these  specimens  were  used  for  the  microhardness  indentation  analyses  prior  to 
the  flexure  test.  The  indentations  at  the  surface  of  the  body  porcelain  may  have 
reduced  these  strength  values.  Relatively  small  values  of  95%  confidence  intervals 
were  obtained  for  the  flexural  strength  of  the  air  cooled  specimens  (Table  4.8). 
For  the  comparison  of  the  three  cooling  treatments  and  the  initial  temperatures, 
these  values  are  plotted  in  Fig.  4.125.  The  results  of  Tukey's  multiple  range  test 
(Table  4.8)  indicate  that  the  flexural  strength  of  disks  tempered  in  air  were 
significantly  greater  than  the  corresponding  values  for  slow-cooled  and  fast-cooled 
specimens.  Shown  in  Table  4.9  is  the  comparison  between  the  tempering  stresses 
at  the  surface  of  body  porcelain  calculated  from  the  viscoelasticity  theory  [66],  and 
the  biaxial  flexure  test  data  as  a  function  of  initial  tempering  temperature.  The 
compressive  surface  stresses  from  the  biaxial  flexure  test  were  obtained  by 
subtracting  the  mean  strength  values  of  the  slow-cooled  specimens  from  the 
corresponding  values  of  the  tempered  specimens  at  each  temperature.  These 
results  indicate  that  the  surface  stresses  predicted  from  the  viscoelasticity  and  the 
structural  relaxation  analysis  are  much  greater  than  those  calculated  from  the 
flexure  test  data.  For  the  body  porcelain  disks  tempered  from  an  initial  temperature 
of  982°C,  the  surface  compressive  stress  calculated  from  the  flexure  test  data  was 
69.4  MPa,  while  the  value  obtained  from  the  viscoelastic  model  was  171.4  MPa 
(Table  4.9).  The  corresponding  value  obtained  from  the  microhardness  indentation 
data  (using  eq.  11)  was  84.1  MPa  (Table  4.6).  Therefore,  under  this  condition, 
the  data  obtained  from  indentation  measurements  showed  better  agreement  with  the 
stress  values  obtained  from  the  flexure  test  than  the  values  determined  from  the 
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Fig.  4.125    Mean  flexure  strength  for  disk  specimens  cooled  in  air  from 
various  initial  temperatures. 


Table  4.9        Comparison  of  the  tempering  stresses  in  body  porcelain 
calculated  from  viscoelastic  model  and  biaxial  flexure  test. 


INITIAL  TEMPERING 
TEMPERATURE  fC) 


CALCULATED  STRESS  (MPa) 


Viscoelastic  Model 


Flexure  Test 


982 

850 

750 
650 


171.4 

165.8 

155.9 
102.0 


69.4 

44.1 

34.8 
14.0 
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viscoelastic  model.  Shown  in  Fig.  4.126  are  the  mean  values  of  biaxial  flexure 
strength  of  tempered  disks  as  a  function  of  initial  tempering  temperature.  The 
strength  value  increases  with  temperature  between  650  and  982°C.  However  the 
influence  of  initial  temperature  decreases  at  higher  values. 


600 


700  BOO  ,  .     900 

INITIAL  TEMPERATURE    (C) 


1000 


Fig.  4. 126    Mean  biaxial  flexure  strength  for  tempered  disks  as  a  function  of 
initial  tempering  temperature. 


CHAPTER  V 
DISCUSSIONS  AND  CONCLUSIONS 


The  strengthening  effect  of  tempering  results  from  the  temperature  gradient 
which  is  introduced  over  the  thickness  of  ceramic  specimens  during  cooling.  To 
optimize  the  effectiveness  of  tempering  in  each  specimen,  it  was  necessary  to 
obtain  the  temperature  versus  time  profiles  under  each  cooling  condition. 
Experimental  and  analytical  results  showed  that  the  most  important  variables  which 
influenced  the  cooling  profiles  and  the  temperature  gradient  within  the  specimens 
were  the  initial  cooling  temperature,  specimen  thickness  and  the  heat  transfer 
coefficient  associated  with  the  cooling  media.  In  general,  at  a  constant  initial 
tempering  temperature,  an  increase  in  thickness  or  a  higher  value  of  heat  transfer 
coefficient  resulted  in  a  greater  temperature  gradient  within  the  specimen,  thereby 
producing  a  larger  magnitude  of  beneficial  compressive  stresses  at  the  surface. 
However,  the  amount  of  undesirable  tensile  stress  at  the  interior  also  increased. 

Cooling  Profiles 

The  major  findings  of  the  experimentally  and  analytically  determined  cooling 
profiles  can  be  summarized  as  follows: 

1.  Tempering  by  a  vertical  flow  of  air  resulted  in  a  greater  temperature 

gradient  compared  with  tempering  by  a  horizontal  flow  of  air. 
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2.  The  advantage  of  tempering  by  a  vertical  air  jet  compared  with  a  horizontal 
air  jet  was  more  pronounced  in  larger  diameter  specimens.  Tempering  by 
a  horizontal  flow  may  result  in  a  nonuniform  local  heat  transfer 
coefficient  because  of  the  change  of  boundary  layer  thickness.  Although 
vertical  tempering  is  also  associated  with  a  local  heat  transfer  coefficient 
that  may  be  nonuniform,  this  situation  can  be  corrected  by  lateral 
movement  of  impinging  jets  relative  to  the  ceramic. 

3.  Although  helium  had  a  greater  cooling  potential  compared  with  nitrogen 
and  air,  the  slopes  of  the  temperature  versus  time  profiles  for  a  2-mm- 
thick  disk  subjected  to  tempering  by  a  vertical  flow  of  either  gas  were 
similar.  This  similarity  in  cooling  profiles  corresponds  to  similar  stress 
distributions  for  these  cooling  media.  Therefore,  the  availability  of 
compressed  air  in  any  dental  laboratory  and  its  low  cost,  support  the 
suitability  of  air  tempering  treatment  of  dental  porcelain  compared  with 
helium  and  nitrogen  gases. 

4.  Experimental  measurement  of  the  temperature  distribution  within  the  dental 

furnace  indicated  that  the  specimen  location  inside  the  furnace  could  have 
a  major  influence  on  the  effectiveness  of  cooling  treatments  employed  for 
dental  restorations. 
In  the  thermal  models  it  was  assumed  that  each  specimen  was  initially  at  a 
uniform  temperature  before  cooling.   However,  as  shown  in  Figs.  2.3  and  4.1, 
initially  a  temperature  gradient  may  be  associated  for  each  specimen,  depending  on 
the  size  of  the  specimen  and  its  location  inside  the  furnace.  Disks  with  a  diameter 
of  16  mm  and  a  thickness  of  2  mm  placed  at  the  hot  zone  region  corresponds  to  a 
more  uniform  initial  temperature  than  disk  with  a  diameter  of  32  mm  diameter  and 
a  thickness  of  6  mm.  However,  in  the  thermal  models  it  was  assumed  that  each 
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disk  was  at  a  uniform  initial  temperature  and  had  a  perfect  insulation  layer  at  the 
bottom  surface. 

Each  cooling  treatment  was  associated  with  a  convective  and  a  radiative  heat 
transfer  coefficient.  As  shown  in  Fig.  4.9,  the  calculated  heat  transfer  coefficients 
for  free  and  forced  convection  (tempering)  were  approximately  independent  of 
surface  temperature.  However,  in  this  study  a  radiative  mode  of  heat  transfer  was 
present  in  each  cooling  condition.  Based  on  the  principle  of  superposition,  the 
coefficient  corresponding  to  radiation  was  added  to  the  convective  heat  transfer 
coefficient  (eq.  24).  Thus,  the  total  film  coefficient  (h)  was  influenced  by  the 
specimen  surface  temperature  at  the  initial  stage  of  cooling  which  becomes 
relatively  constant  at  lower  temperatures.  For  the  calculation  of  transient  tensile 
stresses  by  finite  element  and  viscoelastic  analyses,  an  average  heat  transfer 
coefficient  over  the  entire  temperature  range  was  determined.  This  coefficient, 
which  was  independent  of  the  thermal  properties  of  the  body,  was  calculated  using 
a  thin  copper  disk  with  the  same  diameter  as  the  porcelain  disk. 

The  findings  based  on  1-D  and  2-D  thermal  models  supported  the  proposed 
hypothesis  that,  under  tempering  conditions,  the  influence  of  heat  loss  from  the 
edge  of  the  disk  model  on  the  temperature  distributions  was  negligible.  The 
temperature  versus  time  data  calculated  from  the  1-D  model  for  each  disk  were 
used  in  a  viscoelastic  model  which  was  originally  developed  for  calculation  of 
thermal  stresses  within  a  porcelain  plate  symmetrically  cooled  from  both  sides 
[66].  In  general,  only  small  differences  were  obtained  between  the  calculated  and 
experimental  data.  These  differences  were  probably  due  to  the  assumptions  of 
perfect  insulation  at  the  bottom  surface,  constant  initial  temperatures  for  the 
specimen  and  cooling  media,  and  constant  values  used  for  the  thermal  properties  of 
porcelain  and  the  heat  transfer  coefficients.  Because  of  symmetry  conditions,  the 
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results  of  the  viscoelastic  model  for  a  2-mm-thick  disk  cooled  from  the  top  surface 
and  perfectly  insulated  at  the  bottom  surface  are  the  same  as  the  results  of  a  4-mm- 
thick  disk  symmetrically  cooled  from  both  sides. 

Viscoelastic  Analysis 

Experimental  techniques  such  as  photoelasticity  and  analysis  of  strain 
changes  after  hole-drilling  procedures  are  routinely  used  for  residual  stress 
measurement  of  glass  and  metal  based  structures.  However,  no  direct  and  non- 
destructive experimental  procedure  is  available  to  measure  the  transient  and  residual 
stresses  in  layered  porcelains  because  of  the  inherent  opacity  of  dental  porcelain 
and  the  high  temperature  in  which  stresses  begin  to  develop  during  cooling.  The 
analytical  methods  available  are  based  on  either  the  viscoelastic  or  the  elastic 
behavior  of  porcelain.  In  the  viscoelastic  models  which  are  usually  developed  for 
simple  and  uniform  geometric  shapes,  it  is  assumed  that  porcelain  exhibits  time- 
dependent  viscoelastic  behavior  above  a  temperature  (Tg).  The  elastic  models  are 
based  on  the  assumption  that  a  glass  transition  temperature  (Tg)  exists  below  which 
porcelain  behaves  like  an  elastic  body.  These  models  do  not  consider  the  stresses 
which  are  developed  above  Tg  and  because  of  that  may  be  associated  with  some 
error.  In  the  present  study,  the  model  developed  by  DeHoff  and  Anusavice  [66] 
was  used  to  determine  the  transient  and  residual  stresses  within  porcelain  because 
this  model  includes  the  viscoelastic  effect  of  porcelain  and  accounts  for  stresses 
which  are  developed  above  Tg.  Although,  the  viscoelastic  model  was  originally 
developed  for  a  ceramic  plate  which  may  have  a  different  stress  distribution 
compared  with  a  porcelain  disk  because  of  edge  effects,  the  results  of  the  thermal 
models  indicated  that  under  tempering  conditions,  the  edge  effects  were  negligibly 
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small.  Also  under  convective  cooling,  the  edge  effect  was  reduced  by  adjusting  the 
heat  transfer  coefficient.  The  major  findings  of  the  viscoelastic  analysis  can  be 
summarized  as  follows: 

1 .  The  viscoelastic  model  predicted  relatively  large  transient  and  residual 
tensile  stresses  within  2-mm-thick  and  6-mm-thick  disks,  respectively, 
cooled  by  tempering  from  high  initial  temperatures  in  liquid  media  and  by 
free  convection  in  air.  These  results  were  consistent  with  the  experimental 
cooling  of  several  disks  which  resulted  in  thermal  shock  failure  under  these 
cooling  conditions. 

2.  The  initial  tempering  temperature  had  a  significant  influence  on  the 
magnitude  of  stresses  and  on  the  ratio  of  maximum  residual  compressive  to 
tensile  stresses  at  temperatures  below  750°C.  However,  tempering  from 
initial  temperatures  above  750°C  did  not  improve  the  residual  stress 
distributions  significantly. 

3.  In  general,  larger  beneficial  compressive  stresses  were  associated  with 
greater  values  of  undesirable  tensile  stresses,  independent  of  cooling  rate, 
specimen  thickness,  and  initial  tempering  temperature.  At  initial  cooling 
temperatures  above  850°C,  the  ratio  of  beneficial  compression  to 
undesirable  tension  was  not  strongly  influenced  by  thermal  parameters. 
However,  it  is  possible  that  a  small  variation  in  this  ratio  significantiy 
influences  the  strength  and  fracture  characteristics  of  porcelain,  and  thereby 
controls  the  effectiveness  of  tempering. 

4.  The  depth  of  compressive  stress  was  relatively  insensitive  to  the  rate  of 
cooling  and  initial  tempering  temperature. 

It  should  be  emphasized  that  the  model  proposed  considers  only  a  single  layer 
viscoelastic  material  which  behaves  in  a  linear  fashion.  However,  the  behavior  of 
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porcelain  over  the  entire  range  of  cooling  temperatures  considered  in  this  study, 
may  not  be  linear.  Also,  because  of  the  lack  of  reliable  thermal  contraction 
properties  of  dental  porcelain  at  high  temperatures,  a  glass  transition  temperature  of 
400°C  was  assumed  in  the  viscoelastic  model. 

Finite  Element  Analyses 

To  determine  the  influence  of  loading  orientation,  ceramic  thickness, 
cooling  rate,  thermal  contraction  mismatch,  initial  tempering  temperature  and  glass 
transition  temperature  on  the  stresses  developed  within  a  dental  ceramic  structure, 
these  stresses  were  also  obtained  from  finite  element  analysis  which  was  based  on 
elastic  theory.  The  results  of  the  finite  element  stress  analyses  for  the  bilayered 
ceramic  disks  and  crowns  are  summarized  as  follows: 

A.  Transient  Stress  Analysis  of  Bilayered  Ceramic  Disks 

1 .  Transient  tensile  stresses  produced  within  the  6-mm-thick  disk  were  higher 
than  that  of  the  2-mm-thick  disk  and  this  difference  increased  with  an 
increase  in  the  heat    transfer  coefficient. 

2.  For  a  contraction  mismatch  level  of  +3.2  ppm/°C,  oil  quenching  and 
tempering  from  high  initial  temperatures  produced  relatively  large  transient 
tensile  stresses  at  the  edge  of  the  bilayered  disks. 

3.  In  general,  the  calculated  transient  tensile  stresses  at  the  edge  of  the  disks 
were  greater  than  the  published  tensile  strength  value  for  porcelain.  This 
suggests  that  during  tempering  treatment,  the  2-mm-thick  disk  is  susceptible 
to  thermal  shock  failure.  However,  the  tempering  experiments  did  not 
support  this  finding.    The  transient  tensile  stresses  calculated  from  the  finite 
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element  method  for  a  single  layered  specimens  (Aa  =  0.0)  were  also  greater 
than  the  corresponding  values  from  the  viscoelastic  models.  The  large 
predicted  stresses  from  the  finite  element   analysis  are  probably  due  to  the 
nature  of  the  elastic  analysis  which  does  not  account  for  stress  relaxation. 

B.  Static  Analysis  of  Dental  Ceramic  Crowns 

1 .  Increase  in  ceramic  thickness  from  0.5  mm  to  3.0  mm  had  a  negligible 
influence  on  the  stress  distribution. 

2.  For  a  vertically  oriented  applied  load  of  600  N,  the  magnitudes  of  resulting 
tensile  stresses  were  relatively  small.  However,  when  the  load  was  applied 
along  a  horizontal  direction,  the  tensile  stresses  which  developed  within 
ceramic  rose  above  the  failure  level  of  ceramic.  Therefore,  the  orientation  of 
the  applied  load  is  an  important  factor  in  development  of  large  tensile 

stresses. 

C.  Transient  Stress  Analysis  of  Dental  Ceramic  Crowns 

1 .  Relatively  large  transient  stresses  were  produced  within  the  bilayered 
ceramic  crowns  with  a  thermal  contraction  mismatch  of  3.2  ppm/°C 
when  they  were  subjected  to  free  convective  cooling  and  tempering 
treatment.  Under  this  condition,  the  influence  of  thermal  contraction 
mismatch  on  the  magnitude  of  tensile  stresses  was  more  pronounced  than 
the  influence  of  cooling  rate. 

2.  Reduction  of  initial    tempering  temperature  from  850°C  to  750°C  did  not 
have  a  major  influence  on  the  transient  stresses  significantly. 
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3.       The  glass  transition  temperature  had  a  significant  influence  on  the  magnitude 
of  transient  stresses. 

Quenching  In  Silicone  Oil 

In  our  previous  study  [65],  it  was  shown  that  tempering  treatment 
significantly  reduce  the  potential  for  crack  growth  in  the  surface  of  feldspathic 
porcelain.  However,  we  indicated  that  development  of  surface  tensile  stresses 
during  the  initial  stage  of  cooling  may  be  of  potential  concern.  The  experimental 
results  of  the  present  study  indicate  that  the  transient  tensile  stresses  which  were 
produced  during  tempering  of  a  2-mm-thick  disk  specimen  in  air,  were  generally 
below  the  tensile  strength  of  porcelain. 

The  experiments  performed  by  Kirchner  et  al.  [47]  showed  that  the  mean 
flexural  strength  of  96%  alumina  rods  increased  substantially  when  they  were 
quenched  in  silicone  oils  with  kinematic  viscosities  ranging  from  5  to  12,500 
centistokes.  They  reported  that  the  highest  strength  was  achieved  by  quenching  the 
rods  into  low  viscosity  silicone  oil  at  30  °C  (Table  1.1).  In  this  medium,  the 
flexural  strength  of  two-inch-long  bars  determined  from  four-point  bending  test 
were  increased  by  a  factor  of  2.5  compared  with  the  untempered  specimens,  while 
tempering  by  forced-air  resulted  in  a  strengthening  factor  of  1.6.  In  the  present 
study,  the  effectiveness  of  silicone  oil  as  a  quenching  medium  for  dental  ceramics 
was  tested  by  the  use  of  the  indentation  technique.  Several  quenching  experiments 
were  performed  to  determine  the  maximum  and  minimum  initial  quenching 
temperatures  for  feldspathic  ceramic  specimens  (850°C)  and  for  silicone  oil 
(100°C),  respectively,  to  prevent  thermal  shock. 

In  general,  there  was  no  statistically  significant  difference  between  the  mean 
values  of  induced  crack  size  for  specimens  quenched  in  silicone  oil  compared  with 
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air- tempered  specimens.  The  silicone  oils  used  in  this  study  had  the  same  thermal 
conductivity,  specific  heat  and  density,  however,  they  were  different  in  viscosities. 
Data  reported  in  Tables  4.1  and  4.3  indicated  that  the  silicone  oil  viscosity  in  the 
range  of  50  to  5,000  centistokes  had  no  major  effect  on  strengthening  of  ceramics, 
while  in  the  Kirchner  study  [47]  in  which  alumina  specimens  were  quenched  in 
silicone  oil,  a  change  in  oil  viscosity  from  5  to  350  centistokes  resulted  in  a 
decrease  in  flexural  strength  by  a  factor  of  1 .2. 

The  kinematic  viscosity  of  oil  (y)  has  an  influence  on  the  heat  transfer 
coefficient  (h).  As  y  increases,  h  becomes  smaller,  however  the  temperature 
distribution  and  thus  the  stress  distribution  within  the  tempered  ceramics  depends 
also  on  the  initial  quenching  temperature  and  specimen  thickness.  Perhaps,  in  the 
Kirchner  study  the  initial  specimen  temperature  (1600°C)  and  specimen  thickness 
(3.2  mm)  were  large  enough  for  the  oil  viscosity  to  change  the  residual  stresses. 
However,  in  the  present  study  with  a  initial  specimen  temperature  of  850°C,  an  oil 
temperature  of  100°C,  and  a  specimen  thickness  of  2  mm,  the  change  in  heat 
transfer  coefficient,  because  of  the  oil  viscosity,  was  not  large  enough  to  have  a 
major  influence  on  residual  stresses. 

Microhardness  Indentation  Analysis 

Tables  4.1  and  4.3  show  only  a  small  difference  among  the  mean  crack  sizes 
for  specimens  with  the  three  mismatch  levels  which  were  tempered  in  oil. 
Therefore,  as  reported  in  our  previous  studies  [Anusavice  et  al.  65],  the 
compatibility  stresses  at  the  surface  of  body  porcelain  are  dominated  by  tempering 
stresses. 

Although  the  mean  values  of  crack  size  for  disks  quenched  in  oil  were  not 
significantly  smaller  than  the  corresponding  values  of  specimens  tempered  in  air 
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(Table  4.4),  the  surface  compressive  stresses  calculated  from  equation  (11)  were 
much  higher  for  specimen  quenched  in  silicone  oil.  However,  tempering  in 
silicone  oil  from  initial  temperatures  greater  than  850°C  was  associated  with  a 
higher  risk  of  thermal  shock  failure  for  the  thermally  incompatible  specimens. 
Quenching  of  a  2-mm-thick  bilayered  disk  (Aoc  =  -1.5  or  +3.2  ppmfC)  from 
982°C  in  various  silicone  fluids  at  room  temperature  resulted  in  approximately  10% 
thermal  shock  failure.  However,  when  the  the  initial  specimen  temperature  was 
reduced  to  850°C  and  the  oil  temperature  was  increased  to  100°C,  because  of 
smaller  temperature  gradients,  thereby  smaller  transient  stresses  during  cooling,  no 
failure  occurred. 

The  major  findings  of  the  experimental  analysis  by  using  the  microhardness 
indentation  technique  can  be  summarized  as  follows: 

1 .  The  mean  crack  sizes  at  the  surface  of  body  porcelain  were  significantly 
smaller  for  the  tempered  specimens  compared  with  the  conventionally  air- 
cooled  specimens. 

2.  The  interaction  effect  of  heat  transfer  coefficient  associated  with  the  cooling 
medium,  initial  cooling  temperature  and  the  thermal  contraction  mismatch  on 
the  surface  crack  size  was  statistically  significant.  However,  the  influence  of 
thermal  contraction  mismatch  stresses  on  the  induced  cracks  were  dominated 
by  the  thermal  tempering  stresses. 

3 .  Silicone  oil  viscosity  in  the  range  of  50  to  5,000  centistokes  had  no  major 
influence  on  the  induced  crack  size  and  thus  on  the  magnitude  of  surface 
compression. 

4.  Although,  the  measured  surface  crack  sizes  in  disks  which  were  quenched  in 
silicone  oil  were  smaller  than  the  tempered  specimen  by  air  jet,  quenching  in 
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silicone  oil  from  high  initial  temperatures  was  associated  with  a  greater  risk  of 

thermal  shock  failure. 

In  general,  under  the  same  quenching  conditions,  the  body  porcelain  disks 
(Aoc  =  0)  showed  a  better  resistance  to  thermal  shock  than  the  bilayered  disks. 
Quenching  in  silicone  oil  may  result  in  a  more  favorable  residual  stress  distribution 
because  of  more  uniform  cooling  effects;  however,  the  success  of  quenching  in 
silicone  oil  is  highly  dependent  on  the  proper  selection  of  initial  temperatures  for 
each  specimen  thickness. 

Biaxial  Flexure  Strength 

The  mean  values  of  flexural  strength  of  specimens  tempered  from  initial 
temperatures  of  650,  750,  and  850°C  were  significantly  greater  than  those  for 
slow-cooled  and  fast-cooled  specimens.    For  an  initial  cooling  temperature  of 
982°C,  tempering  of  the  thermally  compatible  system  (Ace  =  0.0),  resulted  in  a 
mean  strength  value  of  116.2  MPa,  while  the  corresponding  value  for  the 
specimens  cooled  by  free  convection  was  39.7  MPa  [65].  These  results  support 
the  hypothesis  that  tempering  improves  the  strength  of  porcelain  structures.  The 
flexural  strength  of  the  tempered  disks  increased  with  the  initial  cooling 
temperature.  However,  the  influence  of  tempering  temperature  was  much  greater 
between  650  and  750°C  (Fig.  4.51).   This  study  showed  that  there  was  only  a 
small  and  statistically  insignificant  difference  between  the  mean  biaxial  strength 
values  for  tempered  specimens  for  initial  specimen  temperatures  between  982°C 
and  850°C.  To  improve  the  effectiveness  of  tempering  and  to  minimize  structural 
distortion  because  of  viscous  flow  and  air  blasting,  tempering  should  be  achieved 
from  an  initial  temperature  of  850°C  or  lower.    The  mean  flexural  strength  of 
thermally  compatible  tempered  disks  (Aoc  =  0.0)  was  1.8  times  greater  than  the 
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corresponding  values  for  the  Act  =  +3.2  ppm/°C  group  [Anusavice  et  al.,  65]. 
Perhaps,  the  incompatibility  stresses  at  the  interior  of  the  disk  which  were  tensile  in 
nature  had  reduced  the  strength  of  these  specimens  despite  the  tempering  treatment. 
The  major  findings  of  the  biaxial  strength  are  summarized  as  follows: 

1 .  The  mean  flexure  strength  of  the  tempered  specimens  was  significantly 
greater  than  that  of  slow-cooled  and  fast-cooled  specimens.  Tempering 
treatment  improved  the  strength  of  porcelain  disks  up  to  3  times  compared 
with  conventionally  cooled  specimens. 

2.  Although  the  influence  of  the  incompatibility  stresses  on  the  size  of  induced 
surface  cracks  was  dominated  by  tempering,  these  stresses  had  a  major 
influence  on  the  flexural  strength  of  tempered  specimens.  This  was 
probably  because  the  effect  of  the  thermal  incompatibility  stresses  were 
more  pronounced  at  the  interior  than  the  surface  region. 

3 .  The  increase  of  initial  cooling  temperature  (TO  from  650  to  850°C  did  not 
have  a  major  effect  on  the  flexural  strength  of  the  conventionally  cooled 
specimens.  However,  for  specimens  tempered  in  air  the  flexural  strength 
increased  by  34%. 

4.  The  flexural  strength  of  the  tempered  disks  increased  with  the  initial  cooling 
temperature.  However,  the  influence  of  initial  tempering  temperature  was 
much  greater  between  650  and  850°C. 

5.  To  improve  the  effectiveness  of  tempering  and  to  minimize  structural 
distortion  because  of  viscous  flow  and  air  blasting,  tempering  should  be 
achieved  from  an  initial  temperature  of  850°  C  or  lower. 

It  should  be  emphasized  that  the  condition  required  to  justify  the  use  of  thermal 
tempering  is  that  the  interior  tensile  stresses  which  are  produced  because  of  applied 
forces  do  not  cause  any  significant  risk  of  failure.     To  strengthen  dental 
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restorations  or  any  other  ceramic  structure  in  which  failure  normally  initiates  at  the 
surface,  tempering  is  a  promising  technique. 

Comparison  of  Predicted  and  Experimentally  Determiner!  Stress 

The  compressive  surface  stresses  predicted  from  the  viscoelastic  model 
were  much  greater  than  those  calculated  from  the  biaxial  flexure  test  data  (Table 
4.9).  For  a  2-mm-thick  disk  tempered  from  an  initial  temperature  of  982°C  in 
air,  the  compressive  stress  determined  from  the  biaxial  test  was  69.4  MPa,  while 
the  corresponding  stress  predicted  from  the  viscoelastic  model  was  171.4  MPa. 
In  explaining  this  difference  one  should  consider  the  validity  of  the  following 
points: 

a)  The  data  obtained  from  the  biaxial  flexure  test  were  not  significantly 

influenced  by  the  flaws  and  microcracks  within  the  body  porcelain  disks. 

b)  The  tempered  specimens  used  in  flexure  tests  were  not  associated  with 
significant  warping. 

c)  The  compressive  stress  calculated  from  the  indentation  test  (eq.  1 1)  showed 
better  agreement  with  the  values  determined  from  flexure  tests  than  with  the 
viscoelastic  model. 

If  these  points  are  correct,  one  may  conclude  that  one  or  more  of  the  following 
factors  which  were  not  considered  in  the  viscoelastic  model  were  responsible  for 
these  differences: 

1 .  The  effect  of  the  cooling  rate  and  specimen  temperature  on  thermal 
contraction  coefficient,  glass  transition  temperature,  and  the  modulus  of 
elasticity  of  porcelain  were  significant. 

2.  The  relaxation  function  and  some  viscoelastic  properties  of  glass 
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inappropriately  assumed  to  be  equivalent  to  the  corresponding  properties  of 
dental  porcelain. 

3.  The  bottom  surface  of  the  porcelain  disk  was  not  perfectly  insulated. 

4.  For  the  viscoelastic  model,  specific  conditions  of  a  two-dimensional  stress 
state  in  which  the  horizontal  and  the  vertical  stress  components  were 
identical  were  assumed.  The  experimental  tempering  condition  may  be 
more  accurately  represented  by  a  true  two-dimensional  viscoelastic  model. 

Recommendations  for  Future  Studies 

This  study  utilized  various  thermal  stress  analysis  techniques  to  test  the 
hypothesis  that  tempering  by  forced  convective  cooling  in  air  is  the  most  effective 
strengthening  method  for  layered  feldspathic  porcelain  structures.  To  simplify  the 
theoretical  and  experimental  analyses,  porcelains  were  prepared  in  the  form  of 
bilayered  disks.  However,  the  major  long-term  goal  is  to  apply  the  thermal 
tempering  technique  to  improve  the  longevity  of  various  ceramic  based  structures 
such  as  dental  ceramic  restorations. 

From  the  foregoing  study  and  our  previous  studies,  it  is  evident  that  the 
strength  of  layered  dental  porcelains  with  simple  geometric  shapes  such  as  disks  or 
plates  can  be  improved  by  thermal  tempering  treatment.  The  major  requirement  for 
development  of  compressive  surface  stresses  by  thermal  tempering  treatment  is  the 
viscoelastic  behavior  of  the  material  at  temperatures  below  the  initial  tempering 
temperature.  Therefore,  these  beneficial  stresses  can  also  be  developed  at  the 
surfaces  of  structures  with  more  complex  geometric  shapes  such  as  crowns  and 
bridges.  However,  because  heat  transfer  coefficient  and  ceramic  thickness  are 
important  variables  which  determine  the  effectiveness  of  tempering,  it  is  necessary 
to  control  the  cooling  rate  at  various  locations  for  a  complex  geometric  structures 
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with  a  non-uniform  thickness.  Perhaps  a  uniform  cooling  can  be  established  by  an 
adjustable  array  of  multi-orifice  nozzles  which  are  positioned  parallel  to  the  surface 
of  the  structure. 

Based  on  the  present  findings,  the  following  recommendations  are  provided  for 
the  future  studies: 

1 .  Experiments  should  be  conducted  to  analyze  the  effectiveness  of  thermal 
tempering  treatment  for  more  complex  geometric  shapes.  To  reduce  the 
number  of  variables,  it  is  recommended  to  start  with  shapes  such  as  a 
semicircular  hoop,  with  a  uniform  thickness. 

2.  Appropriate  viscoelastic  models  should  be  developed  to  analyze  thermal 
stresses  developed  in  multilayered  specimens.  These  models  should  also  be 
capable  of  accounting  for  the  effect  of  cooling  rate  and  specimen 
temperature  on  the  thermal  contraction  coefficient  and  glass  transition 
temperature. 

3.  More  reliable  material  properties  such  as  glass  transition  temperature  and 
viscosity  at  high  temperatures  are  needed  for  porcelain. 

4.  For  tempering  of  structures  with  curved  boundaries,  the  effectiveness  of  a 
multi-orifice  nozzle  which  is  located  in  a  surface  parallel  to  the  surface  of 
structure  should  be  investigated  experimentally. 


APPENDIX  A 

CALCULATION  OF  THERMAL  STRESSES  BY  VISCOELASTICITY  AND 
STRUCTURAL  RELAXATION  [66] 

1-  The  actual  temperature  versus  time  profile  in  a  porcelain  plate  cooled  from 
an  initial  isothermal  state  is  calculated  from  the  solution  for  a  semi-infinite  plate 
subjected  to  convective  cooling: 

N  2  sin  (pn)  cos  (&-*■)  exp  (-    ^  ?  t  ) 

T(x,t)  =  Ta  +  (Tl-Ta)    £   { '   S  '       LL gg       j 

n=l  Pn  +  Sin  pn  COS  pn  ' 

(A-l) 
where  ah  is  the  plate  half-thickness,  x  is  the  distance  from  the  midplane,  k  is 
the  thermal  diffusivity  and  t  is  the  time.  J3n  are  the  positive  roots  of  the 
transcendental  equation  : 

Pntan(pn)-   ^=0  (A-2) 

where  h  is  the  convective  heat  transfer  coefficient  and  K  is  the  thermal 
conductivity  of  porcelain. 

2-  The  fictive  temperature  distribution  is  calculated  from  the  actual  cooling 
profile  by  using  numerical  integration.  The  fictive  temperature  is  defined  as  the 
temperature  from  which  glass  equilibrium  must  be  quenched  in  order  to 
produce  a  particular  structural  state: 

Tf(x,t)  =  T(x,t)-    fexp  K5(*,t)-5(x,u)3b  }  Bj^u}^  ^ 


where    £  (x,t)  =  J  {0  (T(x,u)  ,  Tf  (x,u))}  du  (A-4) 
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AH*     1  2 


and        Ln  ty  (T,Tf)]  =  -  ^r  [  ^  +  ^  -  —  ]  (A.5) 

S  (x,t)  is  referred  to  as  "reduced"  time,  based  on  the  assumption  that  porcelain 
behaves  as  a  thermorheologically  simple  viscoelastic  material.    This  is  a 
material  in  which  the  relaxation  functions  measured  at  two  different  states  are 
related  by  a  shift  in  time  scale  by  a  factor.   In  general,  the  reduced  time,   2; 
(x,t),  at  a  reference  temperature  Tr,  is  the  time  required  for  the  material  to  relax 
the  same  amount  as  in  time  t  at  temperature  T.    Ts  is  the  strain  point 
temperature  (temperature  at  which  the  viscosity  is  10  13-5  N  cm  /m2).  AH*  is 
the  activation  energy  (kJ/mol),  and    R*  is  the  universal  gas  constant  (8.317 
J/mol.  K).  Ty  is  the  change  of  volume  as  a  function  of  time  due  to  step  change 
in  temperature,  b  is  normally  determined  from  measured  relaxation  data,  but  in 
this  study  the  value  of  0.68  which  corresponds  to  glass  [27]  was  used.  %  (x,t) 
also  depends  on  the  fictive  temperature  and  viscosity  of  porcelain.  Once  the 
cooling  profile  and  the  viscosity  as  a  function  of  temperature  for  porcelain  are 
found,  the  fictive  temperature  can  be  determined.    The  integration  variable  u  in 
equation  (A-3)  corresponds  to  any  time  between  t  =  0  and  t  for  the  indicated 
variable. 

3-    The  transient  and  stress  distribution  is  determined  from  the  following 
equations: 


t 
G(x,t)  =  JRn  [£(x,t)  -  S(x,u)  ]      G^X'U)  du  (A-6) 

where: 

°g  (x,t)  =  (  — )  {  e  (t)  -  en,  (x,t) }  (A-7) 
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and 

*n  (t)  -    R(o)  (A-8) 

where  E  and  v  are  the  Young's  modulus  and  Poisson's  ratio  of  porcelain  at 
room  temperature.  R  (t)  is  ,the  axillary  modulus  Rn  (t)  is  a  normalized 
auxiliary  function  which  is  related  to  the  uniaxial  relaxation  function  of 
porcelain  and  is  determined  by  numerical  integration.  e(t)  is  the  strain  in  the 
plate  at  the  time  t.  £&  (x,t)  is  referred  to  as  free  thermal  contraction  from  the 
initial  temperature  state  which  is  determined  from  the  thermal  contraction  data 
and  the  actual  and  fieri ve  cooling  profiles. 


£th  (x,t)  =  a0  [  T(x,t)  -  Tf  (x,t)  ]  +  &1  [  T2  (x,t)  -  Tf2  (x,t)]  + 

a2  [  Tf  (x,t)  -  T,  ]  +   a3  [  Tf2  (x,t)  -  T{2  ]  (A-9) 

At  the  initial  stage  of  cooling  the  stress  and  strain  values  within  the  entire 
specimen  are  zero: 

<*(x.0)  =  0  ,  e(0)  =  0  (A-10) 

and  the  thermally  induced  stresses  at  each  time  must  be  in  internal  equilibrium 
such  that : 

i 
Ja  (x,t)dx  =  0    .  (A_n) 
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This  requires  that : 

l 

fP(x»tn)dx 

£(tn)  =    1 (A-12) 

f[l  +  *„(x,tn.1)]  dx 

0 

Using  the  calculated  values  of  e^  from  the  computer  programs  which  had  been 

developed  in  BASIC  language  for  IBM  microcomputers  [66],  a  solution  for 
PCx.tjj)  at  each  point  in  space  (x)  and  time  (t^)  is  obtained  by  an  iterative 
technique.  Strain,  e(t),  at  initial  condition  is  zero  and  at  a  time  L.  is  calculated 
from  (A-12)  via  numerical  integration.  Then  the  transient  stresses  within  the 
plate  are  determined  from  the  calculated  strain  values  by  using  equation  (A-7). 


APPENDIX  B 
FINITE  ELEMENT  METHOD 

In  the  finite  element  analysis,  a  structure  is  divided  into  discrete  elements  which 
are  connected  at  a  finite  number  of  nodal  points.  The  behavior  of  each  element 
under  internal  and  external  forces  or  thermal  loads  applied  to  the  structures  is 
expressed  by  appropriate  equations  based  on  theory  of  elasticity  and  laws  of 
thermodynamics.  In  the  displacement  formulation  of  finite  element  method, 
unknown  nodal  displacements  are  determined  by  solution  of  a  system  of  linear 
equations  through  the  application  of  the  equilibrium  conditions  between  internal 
and  external  forces  at  the  nodes.  The  nodal  forces  in  thermal  analysis  are  due  to 
temperature  and  in  static  analysis  are  due  to  applied  loads.  The  following  are 
basic  formulations  of  finite  element  method  for  each  problem: 

I.  Thermal  Analysis  ["1 Q2] 

The  general  equation  for  energy  balance  in  the  absence  of  any  internal  heat 
generation  is  : 

pCp¥=v'kvt  (B-l) 

For  constant  thermal  conductivity  and  axisymmetric  condition  this  equation 

becomes: 

dT  82T      i  3T       92T 

5    ~  K[  ~frr+r~dT  +  "ST  3  (B-2) 

If  a  shape  function  {N}  is  defined  such  that  the  temperature  in  any  arbitrary  point 
within  an  element  in  terms  of  the  element  nodal  temperature,  Te,  can  be  obtained 
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from 


T={N}1[Te] 


Then  the  following  matrix  of  temperature  gradient  can  be  written 
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(B-3) 


(S}  =  [B]   {Te} 


were 


(S}  = 


ax 

d  T 

d  T 


[B]  = 


The  heat  flow  rates  at  a  point  can  be  written  as : 


r  d    N  .T 
1    9y    J 

r  ^    N.T 

v  [  s  ]     y 


(B-4) 


(B-5) 


(Q}  = 


Qy 

vQzy 


(kx  0    0  > 

0  ky     0 

\0     0   kzJ 

ax 
a  t 

dy 
a  t 


=  [D]{S}  (B-6) 


Then, 


T 


{Q}  =  [B]1{Q} 


(B-7) 


The  heat  flow  rate  is  related  to  conductivity  matrix  as  follows 


[KJ  (Te}  =  {Qe} 
So  the  conductivity  matrix  becomes  : 


(B-8) 


(Ke)  =     J  [B]T[D]  [B]  d(voi) 

vol 


(B-9) 


Finally  the  matrix  form  of  energy  equation  can  be  written  in  the  following  form 
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dT 


ra  f^ri  +  fK]  {T}  ={qe) 


(B-10) 


where,  [C]  =  Specific  hear  matrix  =  p  CP     J"  {N}  (NjT  d(voi)  (B-ll) 


vol 


ve. 


on 


[K]  -  Effective  thermal  conductivity  matrix  which  can  include  conducti 

convective  and  radiative  surfaces. 

[QE]  =  Effective  nodal  heat  flows  which  can  include  heat  flow  by  convecti 

surface  and  internal  heat  generation. 

An  implicit  time  integration  scheme  can  be  used  to  integrate  the  above  equation 

with  respect  to  time.   Then  a  set  of  linear  simultaneous  equation  is  obtained 

which  is  solved  by  numerical  techniques. 


II.  Static  Analysis 

The  stress-strain  relations  for  isotropic  materials  and  under  a  plane  stress 
condition  are  [103]: 


=  [E/(1-V)2] 


(I  V  0 

1  0 

^Symmetric  (l-v)/2 


£y  -a  AT 

H         Yxy     J 


(B-12) 
General  displacement  expressions  as  a  function  of  nodal  displacements  are  used 
to  describe  the  behavior  of  the  element  within  the  structure  and  are  incorporated 
in  the  strain-displacement  relations: 
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dx 
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(B-13) 


The  corresponding  stress-strain  and  strain-displacement  relations  for  an 
axisymmetric  problem  are: 


faT~\ 


=  [E/(l+v)(l-2v)] 
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(B-15) 


where  u  (x,y),  v(x,y),  u(r,z)  and  w(r,z)  correspond  to  the  displacement  at  any 
point  within  a  finite  element.  For  static  analysis  by  the  finite  element  method,  a 
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set  of  simultaneous  equations  is  solved.  This  has  the  form  : 

IK]    {UN}  =  {F}  (B.16) 

or  (%}  =  [K]-1{F}  (B-17) 

where  [K]  is  the  total  stiffness  matrix  and  is  obtained  by  superposition  of  the 
individual  elements,  {uNJ  is  the  vector  of  unknown  nodal  displacements  and  {F} 
is  the  equivalent  vector  of  applied  loads  which  is  obtained  by  lumping  the  edge 
and  element  loads  at  the  nodal  points. 
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